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This  exploratory  development  provides  the  initial  concepts  for  a digital  imagery 
test  set  from  which  subsequent  detailed  specifications  would  result.  The  diffi- 
cult questions  of  image  evaluation  and  quality  were  explored  with  the  following 
categories  of  performance  suggested  for  the  test  set  criteria;  amplitude  re- 
sponse, frequency  response,  geometric  fidelity,  and  noise.  Both  subjective 
tests  with  an  experienced  photointerpreter  and  objective  tests  using  test  pat- 
tern concepts  were  made.  Computer  programs  were  generated  to  simulate  digital  ^ 
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imagery  degradations  for  the  subjective  tests,  and  these  programs  could  be  ex- 
tended, to  include  sensor  simulations . The  use  of  Moire'  pstterns  provided 
simple  objective  tests  as  a result  of  the  array  structure  of  2-dlmenslonal 
digital  imagery.  Further  development  is  recommended  in  exploring  the  trade- 
offs between  MTF  response  and  aliasing  distortion. 
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ABSTRACT 


The  work  described  in  this  report  provides  basic  information  which  will 
allow  partial  production  of  a standardized  image  quality  test  set.  Factors 
affecting  the  production  of  a standardized  imagery  test  set  concerned  digi- 
tal image  processing  procedures  and  image  reconstruction.  Imagery  consisting 
of  one  scene  from  a Texas  Instruments  FLIR  sensor  tape,  one  scene  from  a 
UPD-4  coverage  of  Griffiss  AFB,  NY  and  two  scenes  from  a Wild  RC-8  Mapping 
camera  were  qualitatively  evaluated  by  a photo  interpreter  from  Rome  Research 
Corp.  (Mr.  Petroski).  The  imagery  was  evaluated  both  in  hard  copy  (digitized 
scenes  with  perturbation  on  film)  and  soft  copy  (digitized  scenes  with  per- 
turbations using  the  RADC/Image  Processing  Facility).  A total  of  144  images 
was  derived  from  four  test  scenes  and  were  evaluated  against  an  original  set 
of  imagery  that  was  used  as  a base  for  comparison. 

The  information  gathered  included  relationships  between  perturbations 
and  image  quality  as  expressed  by  per  cent  utility  under  perturbation,  iso- 
preference contours  (pleasing  to  the  eye)  across  modes  of  perturbation, 
verbal  description  of  effects,  and  image  quality  as  a transitive  (comparative) 
characteristic. 

Results  indicate  that  the  combination  of  sampling  rate  inclusive  of 
source  (sensor)  imagery  characteristics,  structure  (characteristics  of  dis- 
play) and  optimization  (modulation  transfer  function)  of  the  entire  digital 
(processing)  system,  and  the  physiology  of  the  visual  system,  lead  to  bounds 
on  digital  image  quality. 
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SECTION  I 


INTRODUCTION 

Digital  processing  techniques  have  become  commonplace 
in  their  applications  to  imaging  systems  such  that  the  need 
is  long  overdue  for  a digital  imagery  test  set.  The  objective 
of  this  exploratory  investigation  is  to  determine  preliminary 
criteria  for  such  a test  set.  Necessary  prerequisites  which 
are  first  investigated  include  considerations  of  image  quality, 
detailed  analysis  of  the  elements  peculiar  to  digital  imagery 
and  subjective  tests  of  imagery  degraded  by  basic  digital 
pa  rameters . 

The  image  quality  considerations  are  introduced  in  Sec- 
tion II  wherein  an  imagery  source,  digital  processor,  and  ob- 
server are  each  evaluated  as  elements  in  the  imagery  chain. 

The  representation  of  imagery  by  digital  techniques  and  the 
processing  of  digital  imagery  is  placed  in  proper  perspective 
by  considering  the  source  characteristics  such  as  its  spectral 
properties,  and  the  observer  characteristics. 

The  sampling  and  reconstruction  of  two-dimensional  imagery 
signals  are  the  basic  functions  peculiar  to  digital  imagery  and 
these  are  analyzed  in  Section  III.  Aliasing  and  quantization 
noise  are  emphasized. 

Section  IV  provides  the  conditions  for  and  the  results  of 
subjective  tests.  A trained  photointerpreter  is  used  to 
evaluate  the  effects  of  basic  digital  degradations.  Different 
imagery  sources  are  used  and  include  visual,  FLIR  and  SAR. 
Digital  degradations  include  additive  noise,  line  jitter, 
quantization,  and  aliasing. 

Finally,  in  Section  V preliminary  criteria  for  the  digital 
imagery  test  set  are  given  and  several  examples  of  candidate 
tests  are  provided.  The  use  of  Moire'  patterns,  which  are 
Ideally  suited  for  analyzing  the  str  ture  of  digital  image 
patterns,  is  included  in  this  Section. 
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EVALUATION 


Surprisingly  little  is  found  in  the  open  literature  on  digital  image 
quality.  This  preliminary  investigation  including  both  subjective  and  ob- 
jective techniques  forms  a basis  for  the  definition  of  a digital  imagery 
test  set. 

The  recommended  test  set  criteria  should  include  measurements  in  four 
categories  of  system  performance:  non-linear  amplitude  transfer  character- 
istic, non-uniform  MTF  response  (for  frequencies  up  to  one-half  the  Nyquist 
frequency),  non-linear  geometric  response,  and  noise  effects  in  general. 

Digital  imagery,  based  on  its  periodic  or  array  structure,  has  inter- 
esting properties  which  provide  simple  tests  for  many,  if  not  all,  of  the 
categories  of  system  performance. 

A wide  variety  of  alternatives  were  considered  in  the  search  for  dom- 
inant factors  of  digital  image  quality.  These  included  physiology  of  the 
visual  system,  source  imagery  characteristics,  structure  and  optimizations 
of  digital  image  systems,  and  experimental  evaluation  of  a preliminary  sub- 
jective digital  image  test  set. 

A list  of  the  objective  digital  image  quality  factors  was  determined: 

a.  Image  extent  in  number  of  pixels. 

b.  Image  sampling  rate. 

c.  Image  sampling  structure. 

d.  System  MTF  (Modulation  Transfer  Function)  below  the  Nyquist  rate. 

e.  Prefilter  MTF  response  above  the  Nyquist  rate. 

f.  Postfilter  MTF  response  above  the  Nyquist  rate. 

g.  Number  of  quantization  levels. 

h.  Quantization  intervals. 

i.  Gray  level  reconstruction  (false  contouring). 

Preliminary  objective  test  images  suitable  for  non-mechanized  use  by  an 
observer  were  developed  to  determine  the  presence  or  absence  of  many  of  the 
image  quality  factors  listed  above.  The  most  novel  of  these  are  relatively 
unique  being  based  upon  Moire'  pattern  comparisons  against  reference  images. 
A quick  guide  to  system  MTF  up  to  the  Nyquist  limit  was  also  presented. 
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The  following  is  required  in  the  continuation  of  efforts  to  develop  a 
standardized  digital  image  test  set: 

a.  Subjective  evaluation  of  the  complete  set  of  objective  image  qual- 
ity factors  previously  listed. 

b.  Development  of  finished  forms  of  test  imagery  from  the  preliminary 
concepts  described. 

c.  A survey  of  current  state-of-the-art  digital  image  subsystems  in 
reference  to  the  list  of  objective  image  factors. 

d.  Digital  image  processing  procedures  and  algorithms  and  their  in- 
herent affect  upon  image  quality. 

e.  Display  technology  concerning  the  modulation  transfer  function  and 
equipment  characteristics  as  they  pertain  to  image  quality  assessment. 

f.  Image  reconstruction  equipment  as  they  pertain  to  the  production 
of  a standardized  imagery  test  set. 

This  initial  study  has  provided  the  framework  for  a digital  image  test 
set  with  continuing  efforts  recommended  for  detailed  specifications  and 
construction  of  a test  set.  One  development  which  was  started  in  this  ef- 
fort and  worth  exploring  further  is  the  relationship  between  MTF  response 
and  aliased  distortion.  In  addition,  this  study  has  generated  computer 
programs  for  the  subjective  evaluation  of  digital  degradations  which  would 
prove  useful  for  simulating  the  performance  of  the  sensor  portion  of  the 
total  digital  imaging  system. 
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SECTION  II 
THE  IMAGERY  CHAIN 

The  objective  of  this  Section  is  to  provide  introductory 
discussions  of  a digital  imagery  system  in  order  to  define  the 
criteria  for  a Test  Set.  The  essential  elements,  or  functions, 
of  an  imagery  system  are  first  discussed,  and  are  then  followed 
by  considerations  of  performance.  Thus,  the  notions  of  imagery 
evaluation  and  image  quality,  so  important  to  the  test  set 
definition,  are  first  introduced. 

The  linking  of  imagery  information  between  an  input  source 
and  the  output  observer  is  termed  an  imagery  chain.  The  sim- 
plest case  in  which  only  a source  and  observer  are  Dresent  is 
called  the  natural  imagery  chain  as  depicted  by  the  Figure  of 
2-1.  The  observer  is  a human  who  visually  observes  the  source 
representation  be  it  FLIR,  synthetic  aperature  radar,  or  a 
visual  photographic  sensor.  The  human,  being  the  end  user  of 
the  imagery  information,  is  the  final  judge  of  what  constitutes 
"good"  performance.  This  subjective  nature  of  this  performance 
and  the  limited  body  of  knowledge  on  the  psy cho- vi s ua 1 aspects 
of  the  human  observer  presents  a rather  primitive,  although  to- 
date  adequate,  status. 

An  image  processing  function  when  added  between  the  source 
and  observer  constitutes  the  general  imagery  chain  of  Figure  2-2. 
With  the  mushrooming  of  computer  technology  and  the  ability  to 
process,  or  transform,  imagery  information  almost  at  will,  the 
basic  questions  of  imagery  evaluation  and  performance  become 
more  paramount.  For  example,  the  image  quality,  presumed  for 
a given  observer's  task,  can  be  drastically  altered  by  con- 
trolling the  detailed  spectral  response  of  the  imagery.  The 
developments  of  digital  image  processing  have  advanced  to  the 
point  where  basic  questions  of  image  quality,  especially 
digital  image  quality,  need  to  be  answered. 

This  Section  provides  an  introduction  to  the  essential 


Figure  2-1.  The  Natural  Imagery  Chain 


Figure  2-2.  Inclusion  of  Imagery  Processing 
Into  the  Natural  Imagery  Chain 
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functions  of  the  digital  imagery  chain  in  the  following  para- 
graphs. The  source  and  observer  of  the  natural  imagery  chain  are 
discussed  in  Sections  1 1 - A and  1 1 - B respectively.  The  governing 
characteristics  of  the  image  processing  subsystem  under  consid- 
eration here  will  be  digital,  rather  than  photographic,  electric, 
or  acoustic  in  nature.  The  fundamental  characteristics  of  the 
digital  image  processing  subsystem  are  discussed  in  Sections 
1 1 - C , 1 1 - D and  1 1 - E . In  each  of  these  subsections  the  broader 
view  is  taken  in  lieu  of  more  specific  and  detailed  concentra- 
tion upon  the  narrow  present  state-of-the-art  considerations. 
These  narrower  considerations  are  addressed  in  Section  III. 

Once  a digital  processing  subsystem  has  been  inserted  into 
the  natural  imagery  chain,  the  impact  on  performance  must  be 
evaluated.  Imagery  evaluation  procedures  are  considered  in 
Section  1 1 - F . The  question  of  what  constitutes  image  quality 
will  also  be  found  relevant  at  this  point.  The  prime  indica- 
tion of  the  state  of  relevant  knowledge  is  that  no  consistent 
definition  or  measurable  indicator  of  image  quality  has  yet 
emerged.  This  problem  is  further  considered  in  Section  1 1 - G . 

A.  THE  NATURAL  IMAGERY  SOURCE 

Not  totally  as  conjecture,  the  natural  imagery  source  can 
be  pictured  as  shown  in  Figure  2-3.  With  some  abstraction, 
various  naturalistic  relationships  combine  with  laws  of  nature 
to  produce  objects;  i.e.,  material  bodies  of  specified  composi- 
tion occupying  specified  positions  in  space.  These  material 
bodies,  as  a function  of  their  material  properties,  produce 
potential  visual  features.  Usually  these  features  are  surface 
phenomena  (shape,  size,  roughness,  texture)  but  are  occasionally 
volumetric  phenomena  (glow,  trans 1 ucence , density).  Under  self- 
or  external  Illumination,  these  features  are  converted  to  a 
luminous  field,  and  then  are  s ubject  to  various  atmospheric  and 
optical  effects.  Upon  being  subjected  to  a final  perspective 
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transformation,  the  luminous  field  is  considered  an  aerial 
image  suitable  for  Input  to  man  or  machine. 

The  simpler  portions  of  the  natural  imagery  source  have 
been  reduced  to  theory.  Certainly  the  nature  of  projective 
transforms  has  been  understood  for  a few  centuries.  The  various 
a tmospheri c and  optical  effects  have  been  catalogued  and  to  a 
large  degree  manageably  reduced  to  formula.  Basic  laws  of 
interrelationship  of  reflectance,  transmittance  and  absorption 
as  a function  of  illumination  have  been  developed.  In  effect, 
all  of  the  image  modifiers  of  the  image  source  are  reasonably 
well  understood.  What  is  insufficiently  well  understood  is 
the  basic  structural  information  contained  in  the  image.  For 
instance,  it  is  not  presently  possible  to  reduce  to  formula 
answers  to  the  questions:  "How  are  the  image  reflectance  and 
transmittance  functions  of  the  image  features  developed?"  or 

"How  are  the  object  functions  of  the  environmental  relation- 
ships developed?" 

Present  efforts  toward  characterization  of  the  structural 
information  of  an  image  is  simplistic.  Nevertheless  the 
data  so  derived  is  useful  in  determining  the  impact  of  current 
state-of-the-art  image  processing  subsystems.  The  characteri- 
zations include  measurement  and/or  modelling  of  the  image  spec- 
trum, amplitude  distribution,  second-order  statistics,  etc. 

The  spatial  imagery  spectrum  is  currently  the  most  useful 
form  of  structural  data  to  be  extracted  from  the  image.  This 
spectrum  data  which  varies  with  classes  of  source  imagery  is 
particularly  useful  in  determining  desirable  digital  sampling 
frequencies,  filter  rolloffs,  and  bandwidth  reduction  parameters. 
The  approach  taken  here  is  to  intuitively  arrive  at  an  appro- 
priate model  for  image  source  spectra  and  then  demonstrate  that 
the  model  is  consistent  with  the  actual  spectra  of  several 
Image  sources.  One  dimensional  signals  are  assumed;  however, 
the  results  are  easily  extended  to  two  dimensions  as  required. 
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Consider  the  hypothetical  one-dimensional  image  gray  level 
profile  shown  in  Figure  2-4 . The  various  gray  levels  could 
correspond  to  different  objects  in  a scene,  each  having  a 
different  gray  shade  value.  The  transitions  then  are  the 
object  boundaries. 

In  order  to  develop  the  image  source  model,  it  is  necessary 
to  determine  the  probability  of  one  or  more  gray  level  boundary 
transitions  in  an  interval  of  length  t.  In  a scene,  these  gray 
level  transitions  usually  occur  at  the  boundaries  of  material 
objects.  Since  the  boundaries  of  objects  in  one  part  of  the 
scene  do  not  affect  the  boundaries  in  another  part  of  the  scene, 
and  these  boundaries  are  randomly  distributed  in  most  natural 
or  man-made  scenes,  the  boundaries  are  easily  seen  to  be 
Poisson  distributed.  We  shall  assume  that  for  a given  class  of 
imagery,  there  are  a gray  level  transition  boundaries  per  unit 
length;  of  course,  a can  vary  from  one  class  of  imagery  to  the 
next. 

With  these  assumptions  the  probability  of  k transitions  in 
the  interval  t is  given  by 

Pk ( t ) = ( ot ) k e aT/k I , 

in  which  case  it  can  be  shown  that  the  probability  of  no  tran- 
sitions in  r is  P0(x)=e"aT.  (It  should  be  noted  that  we  are 
considering  the  probability  density  function  for  the  number  of 
gray  level  transitions  in  an  interval  t ; this  is  not  to  be 
confused  with  the  probability  density  function  for  the  gray 
levels  themsel ves . ) 

The  autocovariance  function,  C(x),  is  given  by  the 
def 1 nl ti on 

C ( x ) « E{(f(x)-Ef(x) ) ( f ( x+t )-Ef ( x+t ) ) } 

1 • ~ 
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fix) 


on  Assumed) 


where  E{-}  denotes  the  expectation  operator.  It  will  be 
assumed  that  Ef (x)*Ef (x+t)=p.  Let  T denote  the  number  of 
transitions  in  gray  shade  in  (x,x+t).  Furthermore  let  A and 
B represent  the  gray  shade  value  at  x and  (x+t)  respectively. 
Then 


C ( t ) = E [( A-p ) ( B-u ) | T = o]  pr (T=0  ) 

+ E[(A-y)(B-u)|T*l)  pr(Til) 

It  is  easy  to  see  that 

E [(A-p) (B-w) | T=0  j = E(A-y)2  . 

It  is  also  true  that 

eJ{A-(j)(B-u)|T>i]  = 0 . 

This  can  be  shown  as  follows.  It  is  trivially  true  that 
E(A-w)=0.  Furthermore,  the  variates  (A-p)  and  (B-p)  are 
conditionally  independent  given  Til  since  if  a transition 
occurred  between  x and  (x+t)  then  (f(x)-y)  will  be  unrelated 
to  the  value  of  (f(x+x)-u)  for  imagery.  Hence  the  result 
E ( A - u ) (B-w)=0  follows.  It  should  be  noted  that  the  variates 
(f(x)-x)  and  (f(x+T)-p)  are,  in  general,  not  independent:  it 
Is  only  when  they  are  conditioned  by  the  event  Til  do  they 
become  independent. 

It  then  follows  from  the  above  that,  if 

o2  - E [f(x)-w]2 

then 

C(x)  = E[f(x)-p)2  e'a|T| 

so  that 
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S(a)J  - a — 2 ” h{a  e ) • 

u + a 

In  this  case  the  image  spectrum  rolls  offat  6 dB/octave  (20  dB/ 
decade)  and  is  characterized  by  a 3 dB  corner  frequency  of  u>c  = a 
radians/sec.  This  type  of  spectrum  is  referred  to  an  RC  spectrum 
since  it  would  also  result  from  passing  white  noise  through  an 
RC  lowpass  filter.  It  has  thus  been  shown  that  an  RC  spectrum 
is  a plausible  model  for  imagery  source  spectra. 

In  order  to  demonstrate  the  validity  of  the  RC  model,  it 
is  necessary  to  compare  actual  imagery  source  spectra  to  the  RC 
spectrum  model.  The  results  of  real-time  spectral  analysis  of 
analog  television  signals  M indicate  that  television  signal 
spectra  rolloff  at  6 dB/octave  until  the  upper  frequency  limit 
of  the  bandlimited  television  signal  is  reached.  Further  sup- 
port for  this  model  for  source  imagery  spectra  was  obtained 
during  the  study  by  computing  and  examining  the  Discrete  Fourier 
Transform  (DFT)  of  horizontal,  vertical,  and  diagonal  lines  of 
several  images  of  interest  from  Section  IV.  It  should  be  noted 
that  these  DFT's  are  not  the  source  spectra  but  are  based  upon  a 
representation  of  the  source  spectra  obtained  after  scanning, 
sampling,  and  quantizing  the  original  source  image.  This  scan- 
ning and  sampling  process  introduces  secondary  bandl imi tings  due 
to  the  lowpass  nature  of  optical  lens  and  the  finite  bandwidth 
of  photodetectors  and  the  associated  electrical  circuitry.  More 
importantly,  the  original  image  was  undoubtedly  purposely  band- 
limited  (prefiltered)  to  less  than  one-half  the  sampling  rate 
before  sampling  and  quantizing  to  minimize  aliasing  degradations. 
For  these  reasons,  the  spectrum  that  exists  just  prior  to  sam- 
pling and  quantizing  Is  more  properly  called  bandlimited  RC 
(BLRC).  O'Neal ^3  has  derived  an  analytical  expression  for 
the  autocorrelation  of  BLRC  spectra,  with  the  corner  frequency 
(normalized  by  total  bandwidth)  as  a parameter.  Autocorrelation 
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curves  corresponding  to  BLRC  spectra  with  various  normalized 
corner  frequencies  are  shown  in  Figure  2-5.  With  this  in  mind, 
the  DFT  results  obtained  in  the  Image  Processing  Laboratory, 
Figures  2-6  through  2-8,  were  compared  to  the  BLRC  spectrum. 

The  good  agreement  also  supports  the  use  of  an  RC  image  source 
spectrum  model,  if  it  is  assumed  that  secondary  bandlimiting  is 
due  to  the  presampling  filter  in  the  scanner. 

B.  THE  OBSERVER 

We  only  consider  the  human  observer.  Without  too  much 
conjecture,  the  human  observer  can  be  pictured  as  shown  in 
Figure  2-9.  An  image  is  perceived  by  one  of  the  two  sensor 
modes,  to  be  discussed  below.  The  output  of  the  sensor 
operation  is  then  subjected  to  various  feature  and  object 
extractors  which  operate  on  the  sub-conscious  level  of  the 
brain  and  are  thus  automatical ly  performed  by  the  human. 

The  resulting  features  and  objects  are  then  passed  on  to  the 
higher  levels  of  the  brain  where  scene  interpretation  occurs, 
also  on  the  sub-conscious  level.  Where  there  is  time  and 
opportunity,  the  conscious  level  of  the  brain  can  perform 
intellectual  processes  to  further  refine  the  visual  response. 

Little  is  understood  about  the  internal  processes  of 
the  human  observer.  What  is  known  is  something  of  the  physio- 
logy and  performance  of  the  sensor  precesses.  The  sensor  of  the 
human  visual  perception  system  is  the  eye/optic  nerve  combina- 
tion. It  is  well  known  that  the  first  stages  of  this 
combination  are  a lens,  followed  by  a field  of  photodetectors, 
fallowed  by  postprocessing.  The  light  distribution  formed  by 
the  iens  of  the  living  human  eye  on  the  retina  has  been 
measured  and  found  to  be  the. line  spread  function. 
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Autocorrelation  Curves  Corresponding  to  BLRC  Spectra 
lormalized  Corner  Frequency  w /o>  as  a Parameter 


n ( x ) » exp  ( - a | x | ) 


For  white  light  and  a normal  pupil  diameter  of  3 mm , a has  a 
value  of  about  0.7.  The  one-dimensional  MTF  of  the  lens  can 
thus  be  computed  to  be 


HU)  - — P-y 

a +w 


(0.7)2+q>2 


which  is  found  to  be  that  of  a lowpass  filter.  The  -3  dB 
point  occurs  at  approximately  6.6  cycles/degree. 

There  are  two  types  of  detectors  in  the  retina,  rods  and 
cones.  These  each  have  different  functions.  The  rods  function 
primarily  at  low  light  intensities  and  are  insensitive  to  color. 
The  cones  function  at  normal  light  level  and  are  sensitive  to 
color.  The  cones  are  distributed  at  maximum  density  at  the 
center  of  the  retina  and  are  responsible  for  normal  vision. 

They  are  power-law  detectors;  i.e.,  their  output  is  proportional 
to  the  input  light  energy.  This  is  experimentally  expressed 
as  Fechner's  logarithmic  law.  The  base  of  the  logarithmic 
function  has  experimentally  been  placed  at  .29  to  .33  depending 
on  whether  the  source  is  psychological  or  psychophys i ca 1 data. 
The  logarithmic  detector  action  can  be  used  to  explain  the 
phenomena  of  "brightness  consistency"  - the  fact  that  the  same 
object  observed  under  different  illumination  levels  appears  to 
maintain  a near  uniform  brightness.  Figure  2-10(a)  contains  the 
Intensity  profile  of  an  object,  one  under  normal  light  level  and 
the  other  under  five  times  as  intense  an  illumination.  The 
object  to  background  Intensity  difference  Is  1 under  normal 
light  and  Is  5 under  the  more  intense  light.  Figure  2-10(b), 
however,  shows  the  relative  intensities  after  detection  by  the 
eye.  The  object  to  background  and  background  perceived  inten- 
sities now  remain  reasonably  constant  in  both  cases  thereby 
illustrating  the  phenomenon  of  "brightness  constancy." 
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Obviously,  however,  the  eye  detection  system  is  not  an  exact 
logarithmic  function.  Figure  2-11  is  more  exacting.  Figure 
2- 1 1 ( a ) shows  the  subjective  brightness  as  a function  of  the 
image  luminance  for  viewing  of  an  image  in  a darkened  room. 
Figure  2 - 1 1 ( b ) shows  the  same  relationship  for  viewing  of  an 
image  in  a lighted  room.  In  both  cases,  it  must  be  noted  that 
the  eye  can  slowly  adapt  to  the  average  background- image 
luminance.  The  various  curves  of  Figure  2-11  represent  dif- 
ferent adaptation  points.  For  any  particular  adaptation  and 
background  light  level,  there  exists  a roughly  exponential 
curve  relating  image  luminance  and  subjective  brightness. 

From  the  above  it  can  be  shown  that  the  first  two  sensor 
stages  of  the  eye  cannot  account  for  the  observed  MTF  response 
of  the  eye.  Postprocessing  of  the  visual  image  is  required. 
Anatomically  this  can  be  modelled  by  the  backward  inhibitor 
network  of  Figure  2-12.  This  is  one  of  four  such  model  types 
proposed  for  neural  networks.  Assuming  no  self-inhibiting 
action  and  exponentially  decreasing  inhibitory  interaction 
otherwise,  this  postprocessing  stage  can  be  shown  to  be  a 
highpass  filter  with  a response  which  can  be  widely  adjusted. 
Under  these  conditions,  the  human  MTF  can  be  seen  to  be  a 
nonlinear  function  of  the  image  modulation  level  as  shown  in 
Figure  2-13.  The  utility  of  such  a response  is  not  fully 
understood  at  this  time  but  matches  experimental  data.  This 
response  Is  thought  to  reflect  one  extreme  of  the  self-adjust- 
ment behavior  of  the  eye  following  some  internal  criterion  of 
adjustment  for  image  quality.  Figure  2-14  shows  another  such 
Indication  of  the  self-adjusting  behavior  of  the  eye  to  image 
content  as  a function  of  visual  angle,  contrast,  and  background 
brightness.  When  any  combination  of  these  factors  falls  below 
the  curved  surface,  image  discrimination  will  fail,  while  a 
combination  falling  on  a point  above  the  surface  permits  seeing. 


28 


Figure  2-11. 
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Figure  2-13.  Overall  Spatial  Frequency 
Response  for  Several  Values  of 
Sinusoidal  Modulation  m 
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Figure  2-14.  Ability  to  See  as  a Function 
of  Visual  Parameters 
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C.  MAJOR  COMPONENTS  OF  THE  DIGITAL  IMAGERY  PROCESSING 

SUBSYSTEM 

A Digital  Imagery  Processing  Subsystem  is  composed  of  the 
elements  shown  in  Figure  2-15.  Generically,  it  consists  of  a 
conversion  process  which  converts  the  aerial  image  into  a digi- 
tal representation.  Once  the  digital  image  has  been  appro- 
priately processed,  it  is  converted  back  to  an  aerial  image 
for  viewing  by  the  human  observer. 

It  is  most  instructive  to  the  purposes  of  the  study  to  first 
view  the  process  on  an  abstract  level.  Just  as  one  can  plot 
individual  numbers  on  a graph  or  as  positions  in  three-dimenr 
sional  space,  so  can  an  entire  image  be  thought  of  as  occupying 
but  a single  point  in  an  infinite  dimension  space.  Thus  image 
a and  image  6 can  occupy  positions  in  this  space  as  shown  in 
Figure  2-16(a).  The  digital  input  conversion  process  partitions 
this  space  into  regions  such  that  if  an  input  image  lies  in  a 
particular  region  R,  then  the  output  of  the  digital  image  input 
conversion  process  is  the  code  symbol  R.  In  Figure  2 - 16(b) 
three  such  regions  A,  B,  and  C are  shown.  The  symbols  A,  B, 

C thus  a re  three  separate  digital  images.  Now  if  input  images 
a and  B are  input  to  the  digital  input  image  conversion  process, 
the  digital  images  A and  B would  be  respectively  output. 

The  digital  image  processor  most  generally  is  a many  to 
one  mapping  between  digital  images  as  shown  in  Figure  2-  16 ( c ) . 

In  the  example  shown,  the  digital  images  A and  B are  processed 
into  the  digital  images  3 and  2 respectively. 

The  digital  image,  output  conversion  process  is  then  a 
mapping  back  into  image  space  on  a 1 to  1 basis,  i.e.,  one 
digital  image  gives  one  point  (one  aerial  image)  in  image 
space.  This  action  is  shown  in  Figure  2 - 1 6 ( d ) . 

If  the  action  of  the  entire  digital  image  processing  sub- 
system is  sought,  then  the  composite  input  conversion,  digital 
processing  mapping,  and  output  conversion  action  must  be  deter- 
mined. This  is  shown  in  Figure  2-16(e),  the  composite  of  the 
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Figure  2-15.  The  Structure  of  the  Digital  Image  Processing  Subsystem 
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(c)  The  General  Digital  Image  Processing  Mapping 

Figure  2-16.  Abstract  Representation  of  the  Action  Taken  By 
A Digital  Image  Processing  Subsystem 
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(Corresponds 
to  Digital 
Image  2) 


(d)  The  Digital  Image  Output  Conversion  Process 

\ 


(e)  The  Composite  Image  Processing  Subsystem  Action 


Figure  2- 16 (Continued) . Abstract  Representation 
of  the  Action  Taken  By  A Digital  Image 
Processing  Subsystem 
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actions  of  Figures  2-16(b),  (c),  and  (d).  Here  it  is  seen 
that  any  input  aerial  image  lying  in  Regions  A,  B,  or  C will 
result  in  output  aerial  images  I,  II,  or  III  respectively. 

Thus  if  image  a is  input,  image  III  will  be  output,  or  if 
image  B is  input,  image  II  will  be  output. 

From  the  above  discussion  three  important  conclusions 
emerge.  First,  in  a digital  image  processing  subsystem,  many 
input  images  may  result  in  the  same  output  image,  or  conversely, 
any  two  specific  input  Images  may  result  in  the  same  output 
image.  Second,  the  output  image  need  not  bear  any  resemblance 
to  its  input  image,  or  conversely,  it  may  bear  any  resemblance 
to  its  input  image.  Finally,  the  response  of  the  digital  image 
subsystem  to  any  image  (such  as  a test  image)  in  no  way  i$ 
necessarily  indicative  of  its  response  to  any  other  Image. 

While  the  above  conclusions  are  true,  current  useful 
theoretical  approaches  to  the  digital  image  processing  sub- 
' system  are  more  restrictive.  These  processes  are  discussed  in 

Sections  1 1 - D and  1 1 - E . Current  state-of-the-art  hardware 
implementations  are  even  more  restrictive.  These  are  con- 
sidered separately  in  Section  III. 

D.  THE  DIGITAL  IMAGE  CONVERSION  PROCESS 

Since  there  are  no  digital  imagery  sources  in  nature  and 
since  the  perceptual  system  of  the  human  observer  is  ill- 
equipped  to  directly  deal  with  the  digital  image  once  formed, 
there  must  exist  conversion  processes  as  shown  in  Figure  2-15. 

Current  examples  of  hardware  realizations  of  the  input  conver- 
sion process  are  CCD-TV,  half-tone  lithography,  and  laser  beam 
scanners.  Current  examples  of  the  hardware  realizations  of  the 
output  conversion  process  are  digital  CRT  monitors,  half-tone 
printers,  and  laser  beam  reproducers.  Implementation  specific 
data  is  contained  in  Section  III.  The  utility  of  the  present 
discussion  is  that  it  is  more  specific  than  that  of  the  prior 
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section  but  is  general  enough  to  encompass  the  majority  of 
future  state-of-art  developments. 

At  the  usual  theoretical  level,  a digital  image  can  be 
taken  as  a finite  collection  of  finite  valued  integers.  Given 
this,  an  immediate  relevant  question  arises;  "Is  It  possible 
to  develop  digital  image  Input  and  output  conversion  processes 
which  will  in  all  cases  produce  an  output  image  identical  to 
the  input  image?"  The  response  is  NO,  since  there  are  infinite 
numbers  of  possible  input  images  but  only  a finite  number  of 
digital  images.  This  succinct  observation  brings  one  rapidly 
to  the  crux  of  the  matter  - "With  what  quality  does  a given 
output  image(s)  approximate  its  input  image(s)?"  Since  a 
general  definition  of  image  quality  definition  has  not  yet  been 
agreed  upon,  the  above  question  of  image  quality  will  be 
indirectly  pursued. 

In  Figure  2-16(a),  if  image  a is  output  when  image  a is 
input,  then  the  image  quality  of  the  intervening  processes  can 
be  said  to  be  perfect  for  that  image.  If,  however,  an  image  is 
output  which  is  sufficiently  close  to  image  a,  then  the  image 
quality  of  the  intervening  processes  should  again  be  good.  If, 
however,  an  image  is  output  at  a yet  greater  distance  from 
image  a in  image  space,  then  it  is  reasonable  to  assume  that 
the  image  quality  is  not  as  good  since  the  two  images  will  be 
more  dissimilar.  Thus  it  is  possible  to  loosely  associate  image 
quality  with  distance  in  image  space,  in  that  as  the  output 
Image  draws  closer  to  the  input  Image  it  is  reasonable  (but  not 
strictly  necessary)  that  the  image  quality  of  the  output  image 
should  Increase.  It  Is  possible  therefore  to  replace  the 
desire  for  digital  conversion  processes  of  high  Imagfc  quality 
by  the  desire  to  minimize  the  average  distance  between  an  Input 
Image  and  output  image  In  image  space.  This  average  distance 
is  widely  referred  to  as  the  average  mean  square  error  of  the 
output  Image  relative  to  the  Input  image.  It  Is  possible  to 
show  that  the  average  mean  square  error  can  be  made  as 
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E. 


DIGITAL  IMAGE  PROCESSOR 


The  digital  image  processor  (DIP)  is  a digitized  version 
of  the  "Imagery  Processing  Subsystem"  of  Figure  2-2.  Hence  it 
is  typically  a specialized  digital  computer  performing  one  or 
more  algorithms  on  the  digital  data. 

There  are  several  different  types  of  digital  image  pro- 
cessing in  use  today.  One  of  the  first  uses  was  that  of  data 
compression  where  the  DIP  reduces  the  amount  of  digitized  data 
needed  to  represent  an  image;  typically  the  DIP  would  then 
include  a TV  camera,  A/D  converter,  orthogonal  transform 
algorithm,  a communications  channel,  a receiver,  detrans former , 
D/A  converter  and  finally,  a TV  monitor. 

Another  common  type  of  DIP  is  one  used  for  image  restora- 
tion. Here,  the  digital  computer  attempts  to  correct  or  com- 
pensate for  the  optical  imperfections  of  the  imagery  processing 
subsystem.  For  example,  the  DIP  may  introduce  an  effective 
change  in  the  MTF  to  compensate  for  the  TV  camera's  MTF  so  that 
the  visual  high  frequency  content  of  the  scene  is  more  readily 
visible  to  a human  observer. 

Yet  another  type  of  DIP  is  that  of  image  enhancement.  This 
consists  of  performing  operations  on  the  image  to  make  it  appear 
more  pleasing  to  the  eye  or  increasing  its  utility  to  a human 
observer.  To  date,  image  enhancement  consists  mostly  of  ad  hoc 
techniques,  although  the  potential  to  incorporate  knowledge  of 
the  human  visual  system  appears  to  be  good.  One  arbitrary  exam- 
ple of  image  enhancement  is  given  in  Figure  2-17  for  which  the 
output  is  not  necessarily  a replica  of  the  input  but  has  certain 
desirable  features  enhanced. 
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Figure  2-17.  Results  of  an  Arbitrary  Digital 
Image  Processing  Action 
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F.  IMAGE  EVALUATION 

Evaluation  of  the  changes  effected  by  a digital  processing 
subsystem  can  be  accomplished  by  any  of  the  three  measurement 
methods  of  Figure  2-18.  In  contrasting  the  three  methods  it 
must  be  noted  that  the  three  methods  use  significantly  different 
evaluation  structures,  and  thus,  the  evaluation  outputs  can  be 
expected  to  be  quite  different  in  content.  The  choice  of  eval- 
uation technique  will,  therefore,  be  conditioned  on  the  end  use 
of  the  output  evaluation. 

The  mathematical  approach  is  well  suited  to  the  measurement 
of  system  parameters  against  mathematically  defined  specifica- 
tions such  as 

• Modulation  Transfer  Function 

• Si gna 1 - to- Noi se  Ratio 

• Gray  Level  Mapping 

• Resolution 

• Area  Equivalent  MTF 

• Dynamic  Range 

The  subjective  approach  is  well  suited  to  the  cataloging 
of  human  perceived  changes  engendered  by  the  digital  system  such 
as 

• Crispness 

• Brightness 

• Loss  of  Detail 

• Appearance  of  Spurious  Patterns 

t Sparkle 

9 Contrast 

The  objective  approach  is  well  suited  to  the  determination 
of  impairment  or  facilitation  of  the  observer's  reactions  such 
as 


(a)  Mathematical  Evaluation 


(b)  Subjective  Evaluation 


(c)  Objective  Evaluation 


Figure  2-18.  Evaluation  of  a Digital  Imagery  Subsystem 
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• Detection  Time 

• Accuracy  of  Recognition 

• Mensuration  Ability 

G.  IMAGE  QUALITY 

There  exists  no  commonly  accepted  definition  of  image 
quality  in  the  literature.  In  the  case  of  interest  here,  i.e., 
the  insertion  of  a digital  imagery  processing  system  into  the 
natural  imagery  chain,  image  quality  can  reasonably  be  taken 
as  the  percentage  measure  of  how  well  output  images  approxi- 
mates their  respective  input  images  for  the  task  at  hand.  The 
following  problems  arise,  however: 

• What  is  the  appropriate  measure  of  image  quality? 

• What  is  the  set  of  input  images  to  be  considered? 

• What  is  the  task  at  hand? 

An  appropriate  measure  of  image  quality  is  not  yet  recognized  in 
the  literature.  However,  it  is  generally  thought  that  the  par- 
ticular measure  of  average  mean  square  error  between  input  and 
output  image  is  not  appropriate.  For  the  purposes  of  this 
report,  and  particularly  Section  IV,  the  measure  is  indirectly 
defined  as  the  percentage  utility  of  an  output  reconstructed 
digital  image  to  the  input  Imaqe,  as  verbally  expressed  by  a 
trained  photointerpretator . 

The  set  of  input  images  to  be  considered  is  usually 
indirectly  defined  by  examples  or  modelling  such  as  performed 
in  Section  1 1 - A . Section  IV  considers  that  the  input  image 
ensemble  can  be  adequately  specified  by  a few  well  chosen 
examples . 

The  task  at  hand  is  a critical  question.  It  could  range 
from  artistic  evaluation  to  target  detection.  In  Section  IV  the 
task  at  hand  is  indirectly  defined  by  asking  the  trained 
photointerpretator  to  assess  the  images  based  upon  his  expe- 
rience in  photointerpretation. 
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Although  image  quality  has  not  and  can  not  be  definitized 
here,  it  is  interesting  to  note  what  some  of  the  desirable 
properties  of  a measure  of  image  quality  should  be  if  it  is  to 
be  mathematically  useful  and  logically  sound. 

First,  the  condition  of  equality  of  image  quality  should 
be  an  equivalence  relation.  That  is,  the  following  should  be 
obeyed : 

(1)  Symmetric  Relationship 

Any  given  output  image  should  possess  the 
same  image  quality  relative  to  the  input  when 
compared  with  itself.  This  is  patently  true. 

(2)  Reflexive  Relationship 

If  two  output  images  A and  B possess  equal 
image  quality,  it  should  not  matter  that  B is 
being  compared  to  A or  that  A is  being  compared 
to  B.  This  also  is  patently  true. 

(3)  Transitive  Relationship 

If  (1)  two  output  images  A and  B possess 
equal  image  quality  and  (2)  the  two  output  images 
B and  C possess  equal  image  quality,  then  (3)  A 
and  C possess  equal  image  quality.  Many  observers 
question  if  this  is  true.  For  instance,  (1)  if 
an  output  image  A containing  much  noise  seems  to 
possess  upon  comparison  with  image  B containing 
no  noise  but  low  dynamic  range  the  same  image 
quality,  and  (2)  if  image  B seems  to  possess 
upon  comparison  with  image  C containing  jitter 
the  same  image  quality,  then  (3)  it  is  not 
necessarl ly  apparent  that  image  A containing 
noise  will  appear  to  possess  the  same  image 
quality  as  image  C containing  jitter  when 
directly  compared.  This  problem  is  experimentally 
addressed  in  Section  IV. 


The  equivalence  relation  is  basic  to  the  formation  of  any 
mathematical  measure.  Should  the  equality  of  image  quality  not 
exist,  there  will  be  little  means  by  which  a body  of  image 
quality  theory  can  emerge. 

Second,  images  should  be  rankable  in  terms  of  image  quality. 
That  is,  if  the  quality  of  image  A appears  better  than  that  of 
B,  and  independently,  the  quality  of  B appears  better  than  C, 
then  when  directly  compared,  A should  also  appear  better  than  C. 
Obviously,  if  this  does  not  hold,  it  will  be  impossible  to 
obtain  a measure  of  image  quality.  This  condition  is  further 
considered  in  Section  IV. 

Third,  image  quality  should  ideally  be  a metric,  i.e., 
the  closer  the  output  image  is  to  the  input  image,  the  better 
should  be  the  quality.  One  common  example  of  a metric  is  the 
mean  squared  error  (MSE)  criterion.  In  some  applications,  this 
is  likely  not  to  be  true  since  for  instance,  one  can  make  an 
output  image  very  similar  to  the  input  image  except  for  a 
blackening  of  a very  small  area.  The  quality  will  drop  pre- 
cipitantly,  however,  if  that  blackened  area  happens  to  contain 
the  target  of  interest  but  the  MSE  will  be  little  affected  if 
the  blackened  area  originally  contained  the  target  of  Interest 
in  a dark  background. 
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SECTION  III 


SCANNERS  AND  REPRODUCERS 

In  this  section,  the  theoretical  and  practical  aspects  of 
the  scanning  and  reproducing  of  source  images  are  considered. 

The  image  degradations  introduced  by  spatial  sampling  and 
reconstruction  are  discussed  in  Section  1 1 1 - A followed  by  a 
discussion  of  the  image  degradation  due  to  quantizing  noise  in 
Section  III-B.  Finally,  some  of  the  practical  aspects  of  image 
degradations  introduced  in  the  scanning  and  reproducing  process 
are  examined  in  Section  III-C. 

A.  SPATIAL  SAMPLING  AND  RECONSTRUCTION 

Two  fundamental  operations  are  performed  by  the  scanner 
in  any  digital  imagery  system:  two-dimensional  spatial  sampling 
and  quantizing  each  of  the  resulting  analog  values  to  a finite 
number  of  bits;  i.e.,  a form  of  the  general  finite  orthonormal 
approximation  and  digitization  process.  Section  1 1 - E . The 
result  of  these  two  operations  is  a digital  representati on  of 
the  source  image  that  is  discrete  in  both  space  and  amplitude. 

The  digitized  image  may  be  thought  of  as  a two-dimensional 
matrix  of  sample  values,  in  which  each  sample  value  is  allowed 
to  have  one  of  a finite  set  of  values.  For  convenience,  the 
Image  is  considered  to  be  one-dimensional  for  the  remainder  of 
this  section  since  the  conclusions  are  readily  extended  to  the 
two-dimensional  case.  In  addition,  some  of  the  discussion  is  in 
1 n terms  of  time  signals  rather  than  spatial  images.  The  corre- 
spondence between  one-dimensional  spatial  images  and  time-signals 
is  s t ra i gh t forwa rd . 
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1.  The  Nyquist  Sampling  Theorem  and  Time-Limited  Real  Signals 

The  Nyquist  Sampling  Theorem  states  that  any  function  f(t) 
containing  frequencies  no  higher  than  W can  be  perfectly  recon- 
structed from  its  samples  taken  — seconds  apart.  Signals 
which  are  sampled  less  frequently  than  f$=?W  samples/second 
suffer  a type  of  degradation  called  "foldover"  or  aliasing,  in 
which  high  source  frequency  components  masquerade  as  lower 
frequency  components  in  the  reconstruction  process.  A subtlety 
involved  in  any  practical  application  of  the  Nyquist  Sampling 
Theorem  is  due  to  the  fact  that  any  real-world  signal  of  prac- 
tical interest  is  finite  in  extent  (time  or  space)  and  cannot 
therefore  be  truly  bandlimited  in  the  strict  sense.  This 
implies  that  the  sampled  signal  cannot  be  reconstructed  exactly 
for  any  finite  sampling  rate;  however,  the  energy  of  the  error 
signal  can  be  made  as  small  as  desired  by  increasing  the  sam- 
pling rate.  A result  due  to  Landau  and  PollakL3-!  quantifies 
the  tradeoff  between  sampling  rate  and  error  signal  energy 
as  follows: 

if  x ( t ) is  time-limited  to  T seconds  and  the 
signal  energy  that  falls  outside  the  bandwidth 
B Hz  is  less  than  or  equal  to  k times  the  total 
signal  energy,  then  the  energy  of  the  error 
signal  between  x(t)  and  an  approx imation  con- 
structed using  2TB  signal  samples  is  less  than 
12  kE. 

Using  the  results  of  reference  [3]  , Slepian^J  has  proposed 
meaningful  definitions  for  bandlimited  and  time-limited  signals 
using  the  notion  of  "d i s t i ngui sha b i 1 i ty  at  level  E." 

One  of  the  measures  of  image  quality  that  is  directly 
related  to  the  spatial  sampling  rate  is  spatial  resolution; 
another  measure  of  image  quality  related  to  spatial  sampling 
rate  is  the  degree  of  aliasing  degradation  present.  The 
relationships  and  tradeoffs  among  image  bandwidth,  image 
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resolution,  sampling  rate,  and  aliasing  effects  a re  classified 
in  the  following  paragraphs. 

2.  Image  Resolution  and  Aliasing 

The  original  image  to  be  transmitted  possesses  a continuous 
spectrum  of  spatial  frequencies  much  as  shown  in  Figure  3-1. 

When  this  image  is  (ideally)  sampled  at  the  scanner,  the  sam- 
pling process  drastically  modifies  the  original  image  spectrum 
as  shown.  The  sampled  spectrum  is  composed  of  an  infinite 
number  of  original  image  spectra  spaced  at  intervals  f$  (the 
sampling  rate)  apart.  When  the  transmitted  sampled  image  is 
received  at  the  recorder,  it  is  desired  to  display  not  the 
sampled  image  itself,  but  a reconstruction  of  the  original 
image.  This  is  accomplished  by  trying  to  retrieve  the  spectrum 
of  the  original  image  by  isolating  the  fundamental  spectrum 
(crosshatched  in  Figure  3-1  and  other  figures)  of  the  sampled 
image.  This  may  be  accomplished  by  placing  the  sampled  image 
through  the  lowpass  reconstruction  filter  shown  in  Figure  3-1 
which  "removes"  all  image  data  above  spatial  frequency  f s / 2 to 
provide  the  reconstructed  spectrum  for  display.  Note  that  if 
the  repeated  spectra  of  the  sampled  image  overlap  as  shown 
in  Figure  3-1,  then  the  isolation  action  of  the  reconstruction 
filter  is  imperfect  and  some  spurious  detail  called  aliasing, 
due  solely  to  the  sampling  process,  enters  the  displayed  image. 
Since  it  does  not  visually  relate  to  the  basic  image  data  this 
spurious  detail  is  normally  perceived  as  image  noise  degrading 
the  quality  of  the  display  image. 

A graphic  example  of  the  aliasing  effect  is  shown  in  Figure 
3-2,  in  which  a sinusoid  of  1 i nea rly- i ncreasi ng  (with  space) 
spatial  frequency  (a  "chirp"  waveform)  which  has  also  been 
sampled  at  a fixed  frequency  f$  is  displayed  on  a recorder 
hardcopy  output.  In  this  figure,  the  abscissa  is  spatial 
location  and  since  the  source  frequency  is  increasing  with 
time,  a spatial  frequency  can  be  identified  with  each  spatial 
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Figure  3-1.  Sampled  Imagery  System  with  Aliasing 
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location  on  the  abscissa.  Because  the  MTF  of  the  recorder 
attenuates  all  spatial  frequencies  greater  than  f $ / 2 , only 
the  aliased  frequency  is  visible  on  the  recorder  output.  In 
Figure  3-2  the  spatial  frequency  at  spatial  location  xq  is 
zero.  For  source  frequencies  higher  than  f$/2  (at  Xj  in  the 
figure),  the  MTF  of  the  recorder  attenuates  the  fundamental 
(principal  alias)  and  passes  the  aliased  frequency.  At  X£, 
the  source  frequency  is  equal  to  fg  and  the  aliased  frequency 
is  zero  while  at  x4  the  aliased  component  associated  with  2f$ 
has  zero  frequency;  the  zero  frequency  display  is  evidenced 
as  the  three  widest  black  stripes  in  the  figure  at  Xq  , x? , and 
x 4 . Although  it  is  incidental  to  the  discussion  of  aliasing, 
note  the  change  in  contrast  along  the  vertical  axis  of  Figure 
3-2  which  indicates  that  aliasing  is  not  a function  of  contrast. 

In  general,  there  are  only  two  ways  theoretically  to 
remove  any  spurious  aliasing.  The  first  is  to  move  the  spatial 
sampling  frequency  up  to  an  f$  which  is  twice  the  spatial 
frequency  content  of  the  original  image  and  to  move  the  recon- 
struction filter  lowpass  cutoff  to  one-half  that  amount.  The 
results  are  seen  in  Figure  3-3  for  this  case.  Since  the 
repeated  copies  of  the  original  image  spectrum  within  the  sam- 
pled image  do  not  now  overlap,  they  can  be  correctly  separated 
by  the  lowpass  reconstruction  filter  without  incurring  any 
image  degradation  due  to  aliasing.  The  only  problem  with  this 
approach  is  that  is  is  counterproductive  in  terms  of  conserving 
image  transmission  bandwidth  as  it  requires  the  transmission  of 
more  data  due  to  the  higher  sampling  rate.  The  second  approach 
places  a prefilter  prior  to  the  spatial  sampler  as  shown  in 
Figure  3-4  in  order  to  "remov.e"  all  original  image  spatial 
frequencies  higher  than  one-half  the  sampling  frequency  f . 

This  prevents  overlapping  of  the  repeated  copies  of  the  pre- 
filtered original  image  spectrum.  Since  these  repeated'copies 
do  not  overlap,  they  can  now  be  correctly  separated  by  the 
lowpass  reconstruction  filter  without  incurring  any  image 
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degradation  due  to  aliasing.  The  net  result  has  been  to  reduce 
the  system  resolution  and  eliminate  aliasing  degradation.  As 
pictured  in  Figure  3-4,  the  final  reconstructed  image  possesses 
none  of  the  degradations  of  aliasing,  although  the  resolution 
is  poorer  when  compared  to  that  of  the  original  image. 

3.  Aliasing  Considerations  in  Image  Resolution  Reduction 

In  many  practical  image  scanning  and  recording  systems 
there  exists  a requirement  for  more  than  one  available  resolu- 
tion. This  requirement  can  be  met  by  always  sampling  at  the 
highest  rate  required  and  reducing  the  number  of  pixels  in  the 
digitized  image  whenever  lower  resolution  is  required.  There 
are  several  ways  in  which  the  number  of  digitized  image  pixels 
could  be  reduced;  the  simplest  of  these,  involves  merely 
retaining  only  the  nth  pixel  of  every  nth  line.  Theoretically 
this  action  is  equivalent  to  inserting  a supplementary  spatial 
sampler  into  the  system  of  Figure  3-4  at  the  location  shown  in 
Figure  3-5.  The  result  of  this  reduced  sampling  rate  is  a 
repetition  of  the  imagery  spectrum  at  frequency  intervals  fg/n 
which  causes  a heavy  overlapping  of  the  fundamental  image 
spectrum  (crosshatched  in  Figure  3-5).  Since  there  is  no  re- 
construction  filter  which  can  isolate  the  fundamental  image, 
there  will  be  considerable  spurious  detail  within  the  displayed 
image  as  shown  in  the  frequency  spectrum  of  the  reconstructed 
image  in  Figure  3-5.  This  aliasing,  of  course,  visually 
appears  as  noise-like  detail  degrading  the  image  quality. 

What  modifications  to  this  redundancy  removal  technique 
are  appropriate  in  order  to  reduce  the  excessive  aliasing 
effects  Incurred?  From  the  sampling  theory  discussion  of 
Figure  3-4,  it  is  clear  that  an  appropriate  prefilter  must  be 
provided  with  a cutoff  of  0.5  f$/n  and  that  the  reconstruction 
filter  cutoff  must  also  be  decreased  to  0.5  f$/n. 

If  one  now  abandons  the  one-dimensional  signal  assumption 
and  generalizes  the  results  to  a two-dimensional  spatial  image. 
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Figure  3-4.  Sampled  Imagery  System,  Prefiltered  to  Prevent 

Aliasing 


it  is  clear  that  the  "redundancy  removal"  filter  function  of 
Figure  3-5  must  be  a two-dimensional  filter  function.  Imple- 
menting such  a two-dimensional  filter  function  requires  enough 
memory  to  "sequentially  assemble"  the  pixel  values  in  a suffi- 
ciently large  neighborhood  of  each  image  pixel  as  the  image  is 
point  scanned.  These  neighboring  pixel  values  are  then  pro- 
perly weighted  and  summed  to  provide  a single  output  value 
at  each  sampling  instant. 

The  weighting  function  corresponding  to  the  ideal  lowpass 
filters  shown  in  Figure  3-6  is  a two-dimensional  — ^ -x-  function, 
which  is  infinite  in  extent.  Truncated  s--y ---  weighting  func- 
tions could  be  used  to  approximate  the  lowpass  filter  charac- 
teristic. As  an  example  of  the  aliasing  degradation  that  can 
be  encountered  in  practical  resolution  reduction  operations, 
assume  that  the  weighting  function  (one-dimensional  case)  is 
uniform  and  encompasses  a number  of  samples  corresponding 
to  twice  the  one-dimensional  resolution  reduction  ratio,  i.e., 
for  a 4:1  resolution  reduction,  eight  samples  would  be  uniformly 
wei ghted . 

The  model  for  such  a resolution  reduction  approach  and 
corresponding  spectra  at  key  points  are  shown  in  Figure  3-7.  In 
Figure  3-7,  hj(t)  is  the  digital  decimating  filter  and  h^Ct)  is 
the  interpolating  filter  which  includes  the  model  for  the  analog 
Interpolation  process  (D/A  converter,  s amp  1 e- and- ho  1 d post- 
filter). Note  the  presence  of  aliasing  between  0 and  f$/8  and 
the  unwanted  spectral  components  between  fs/8  anf  f$/2. 

An  alternative  decimation  approach  that  requires  only  one- 
fourth'(for  the  two-dimensional  case)  of  the  computational 
complexity  of  that  shown  in  Figure  3-7  involves  using  only 
four  input  samples  rather  thaneight  (see  hj(t)  of  Figure  3-7) 
to  obtain  each  "representative"  pixel.  The  penalty  for  this 
reduced  computation  requirement  is  that  the  aliasing  becomes 
more  severe.  The  principal  (desired)  alias  of  S^(f)  in 
Figure  3-7  is  contaminated  only  by  side  lobes  of  the  remaining 
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Figure  3-7.  Non  ideal  Two-Dimensional  Low -Pass  Filtering 
Using  Uniformly  Weighted  Input  Values 


aliases  for  eight-sample  decimation,  whereas  four-sample 
decimation  causes  the  principal  alias  of  SE(f)  to  be  contamin- 
ated by  the  main  lobe  of  the  nearest  alias,  in  addition  to  the 
side  lobes  of  all  the  remaining  aliases. 

B.  QUANTIZATION  OF  IMAGE  SOURCES 

In  order  to  transmit  a representation  (approximation)  of 
an  aerial  image  to  another  location  using  a digital  channel, 
the  source  must  first  be  sampled  and  then  quantized.  The 
image  degradations  introduced,  viz.,  aliasing  and  resolution 
reduction,  in  the  sampling  process  were  discussed  in  Paragraph 
1 1 1 - A . The  degradations  incurred  as  a result  of  the  scanner 
quantizing  Figure  3-8(a)  and  reproducer  gray  level  mapping, 
Figure  3-8 ( b ) , processes  taken  together.  Figure  3-9,  will  be 
treated  here.  Three  approaches  of  increasing  complexity  and 
performance  are  described.  This  is  followed  by  a discussion 
of  the  quantizing  performance  criterion  and  the  determination 
of  optimal  quantizing  characteristics. 

1.  Three  Approaches  to  Quantizing 

Assume  that  the  sequence  of  numbers  that  results  from  the 
sampling  of  the  image  source  produces  Gaussian  samples  that  are 
independent.  The  simplest  way  to  convert  this  sequence  of  image 
sample  values  into  a sequence  of  bits  is  to  utilize  a quantizer 
having  M uniformly  spaced  output  levels  to  map  each  input  sample 
value  into  one  of  the  M output  levels. 

The  quantizing  distortion  that  results  is  a function  only 
of  the  number  of  levels  M and  the  spacing  between  levels.  This 
simple  approach  of  digitizing  the  source  sequence  on  a sample- 
by-sample  basis  requires  r bits/sample  where 

r i log2  M . 
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(a)  Scanner  Gray  Level  Quantizatii 
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(b)  Reproducer  Gray  Level  Mapping 


Figure  3-8.  Scanner  and  Reproducer  Gray  Level  Mappings 


I 


z 

o 


3 >» 
O <T3 
L. 
O 


62 


Figure  3-9.  Development  of  the  Composite  Gray  Level  Mapping 


A more  complex  quantizing  approach  treats  the  quantizer 
output  as  a M-ary  source  and  assigns  binary  codewords  to 
sequences  or  blocks  of  quantizer  outputs.  As  the  number  of 
quantizer  output  symbols  per  block  becomes  large,  the  minimum 
number  of  bits/sample  can  approach  log^M,  where  M need  not  be 
an  integer  power  of  the  base  two.  Note  that  the  two  approaches 
above  yield  an  identical  number  of  bits/sample  in  the  special 
case  M=2r,  integer  = r. 

The  third  and  most  complex  of  the  three  quantizing  ap- 
proaches takes  advantage  of  the  fact  that  a di screte- ti me 
Gaussian  source  will  not  have  equiprobable  outputs.  For  any 
source  with  a nonuniform  amplitude  probability  density  function, 
the  M-level  quantizer  output  will  exhibit  an  entropy  that  is 
less  than  log^M  bits/sample.  The  output  of  such  a source 
quantizer  can  be  entropy  coded  with  a number  of  bits/sample 
that  approaches  the  informational  entropy  of  the  quantizer 
output. 

2.  The  Quantizing  Performance  Metric 

The  performance  metric  for  quantizers  is  usually  taken 
as  the  mean-squared  error  (MSE)  between  the  quantizer  input 
(source)  and  quantizer  oi  iput.  Although  it  is  recognized  that 
mean-squared  error  is  an  inadequate  metric  of  subjective  image 
quality,  it  is  a useful  measure  of  image  quality  in  the  sense 
that,  if  the  MSE  is  made  sufficiently  small,  the  subjective 
image  quality  will  be  perceived  as  being  "good"  rather  than 
"fair"  or  "poor".  Whichever  of  the  three  quantizing  approaches 
discussed  is  chosen  for  quantizing  imagery,  the  quantizing 
characteristics  can  be  specified  so  that  the  quantizing  dis- 
tortion (MSE)  is  its  minimum  possible  value.  All  that  is 
required  is  a knowledge  of  the  amplitude  probability  density 
function  p(x)  of  the  source,  which  can  be  used  to  specify  the 
quantizer  threshold  values  xi  and  representation  values  y ^ 
by  solving  (numerically)  the  simultaneous  equations 
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/*x  i + 1 

y , = / xp(x)dx  i = 1 , 2 . . . 

xi 

and  using 

xi  = (yi+yi-i>/2  i = 2,3  ... 


It  is  clear  from  the  first  of  Ihese  equations  that  the  optimum 
representation  level  yi  is  the  centroid  of  the  area  of  p(x) 
between  x^  and  x.j  + ^. 


3.  Source  Coding  Subject  to  a Fidelity  Criterion: 
Shannon's  Ra te- Di s tort i on  Theory 


For  all  useful  quantizing  characteristics  the  magnitude 
of  the  MSE  decreases  as  the  number  of  bits/sample  increases, 
i.e.,  the  channel  bit  rate  can  be  traded  for  quantizing 
distortion  power  e . 

The  lower  bound  for  this  tradeoff  is  given  by  Shannon's 
rate-distortion  theory  result 


R(bi ts/sample ) = 1/2  log^  (o2/e^) 


whe  re 

2 

a = input  signal  variance 

and 

2 2 

a /e  = s i gna 1 - to-quanti zi ng  noise  ratio  . 

It  is  of  interest  to  compare  the  performance  of  the  three 
quantizing  approaches  described  previously  to  the  rate-distortion 
bound  above.  Figure  3-10  (taken  from  [6])  indicates  that  entropy 
coding  of  each  uniformly-quantized  image-source  sample  (see  the 
HQ ( v ) curves)  performs  within  1/4  bit  of  the  rate-distorti on 
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2 2 

bound.  The  curves  in  Figure  3-10  designated  M(e  ) and  B(e  ) 
refer  to  block  coding  and  s amp  1 e- to- s amp  1 e quantizing,  respec- 
tively, and  show  that  even  these  relatively  simple  techniques 
exhibit  rate -distort ion  performance  that  is  reasonably  close 
to  the  Shannon  bound  for  moderate  values  of  MSE. 

As  might  be  expected,  nonlinear  (nonuniform)  M-level 
quantizers  yield  smaller  MSE  than  linear  (uniform  ) M-level 
quantizers  for  any  nonuniform  probability  density  function, 
although  the  difference  is  relatively  small  for  Gaussian 
si gnal s (see  Max  t7^ ) . 

Curiously  enough,  if  the  quantizer  output  is  to  be  entropy 
coded,  uniform  (i.e.,  linear)  quantizing  gives  lower  MSE  than 
nonuniform  quantizing,  as  proved  first  by  Gish  and  Piercet8] 
and  then  by  Woodt9].  This  result  has  been  shown  to  hold  for 
any  reasonably  smooth  probability  density  function,  not  just 
Gaussian  distributed  signals. 

Because  of  their  relative  simplicity  of  implementation , 
linear  quantizers  without  entropy  coding  are  often  used  to 
digitize  image  sources  on  a sample-by-sample  basis.  It  can 
be  easily  shown  that  the  quantizing  noise  power  is  given  by 

p 

A / 12  (A=quantizer  stepsize),  regardless  of  the  probability 
density  function,  if  the  number  of  quantizer  levels  is  not  too 
small.  The  s i gnal - to-quan ti zi ng  noise  ratio  (SNR)  can  then  be 
obtained  if  the  stepsize  can  be  expressed  in  terms  of  signal 
power.  One  commonly  used  quantizer  is  the  4o  quantizer  in 
which  the  maximum  quantizer  output  level  is  designed  to  be 
four  times  the  rms  signal  level;  i.e.,  Vpea|<  = 4<7-  Using  the 
relationships 


A = total  quantizer  range/number  of  levels  -1 
2V. 


8o 


( 2n" 1 ) 


bits/sample 


and 


Nq  = A / 12 
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the  resulting  expression  for  SNR  is 


! 


? 12  ( 2n  ) 2 
SNR  = a / Nq  = ^ 


or , in  dB  , 


SNR  (dB)  = 10  log  SNR  = (6n-7.27)  dB 

Of  course,  this  expression  is  only  accurate  for  signals  whose 
rms  value  a is  one-fourth  the  maximum  quantizer  output  level. 

The  corresponding  derivation  for  full  load  sinusoidal 
signals  yields 

SNR(dB)  = ( 6n  + l . 76  ) dB  . 

A problem  encountered  in  practice  is  that  real  world 
signals  are  seldom  truly  stationary  with  respect  to  either 
signal  power  or  amplitude  probability  density  function.  This 
fact  results  in  suboptimum  quantization.  This  effect  can  be 
ameliorated  by  designing  a suboptimum  logarithmic  p-law  quan- 
tizer that  performs  well  with  a wide  variety  of  probability 
density  functions  over  a wide  dynamic  range  of  input  power 
levels.  The  performance  of  p-law  quantizers  is  given  (io]  by 

SNR  (dB)  = 20  log  / — \ dB 

\ 1 og ( 1 + p ) / 

and  for  the  commonly  used  p = 255  system,  this  reduces  to 

SNR(dB)  = ( 6 n- 2 . 86 ) dB  . 

The  rate  distortion  bound  on  SNR  performance  can  be 
obtained  from 


\ 

l 
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m 

" sample  * R (e2>  ' * '°92  <°2/y  ■ >>  '092(SNR) 
which  implies 


SNR  (dB)  = 10  log1Q  (22n)  = 6n(dB)  . 

Since  one  bit  is  equivalent  to  6 dB,  the  4 a quantizer  and  y-law 
( u = 2 5 5 ) quantizer  perform  within  1.2  bits  and  0.48  bit  of  the 
rate  distortion  bound,  respectively. 

C.  PRACTICAL  DIGITAL  IMAGE  DEGRADATION  MECHANISMS 

It  has  been  shown  above  and  in  Section  1 1 - C that  any  digi- 
tal image  conversion  process  must  introduce  image  degradations. 
It  was  also  shown  that  theoretically,  these  degradations  could 
l be  made  as  insignificant  as  desired. 

In  addition  to  any  theoretical  degradation  incurred,  there 
normally  exists  a large  repertory  of  degradations  encountered 
due  to  practical  implementation  considerations.  For  instance, 
any  structural  or  parameter  change  from  optimal  will  result  in 
degradation.  These  changes  can  be  due  to  noise,  spurious  sig- 
nals, device  limitations,  design  choice,  the  environment  and 
component  failures.  In  the  most  general  state-of-the-art  case, 
general  orthonormal  operation,  the  effects  of  these  changes  from 
optimality  are  blurred  across  the  image.  Thus  the  degradation 
will  subjectively  appear  as  a general  worsening  of  image 
qual i ty-general ly  with  little  specific  to  point  to  other  than 
loss  of  resolution,  contrast,  and  mudding  of  image  detail. 

In  the  case  of  spatially  sampled  orthonormal  conversion, 
the  effects  of  image  processing  system  changes  from  optimality 
are  much  more  distinct  and  differentiated  than  in  the  case  of 
general  orthonormal  conversion.  Since  spatially  sampled  ortho- 
normal conversion  Is  used  in  the  vast  bulk  of  current  state-of- 
the-art  imagery  scanners  and  recorders,  it  becomes  Important  to 
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look  more  closely  at  their  practical  degradation  effects. 

First,  note  that  the  orthonormal  basis  used  is  the  impulse 
function.  As  a result,  it  will  typically  be  found  that  any 
system  change  from  optimality  will  affect  only  a very  small 
region  of  the  image.  This  is  largely  due  to  the  fact  that  in 
this  system  any  unit  of  data  passing  through  the  process  con- 
tains information  only  about  a corresponding  local  region  of 
the  image.  If  this  data  is  perturbed  by  a system  change,  it 
is  normally  seen  as  a localized  phenomenon  (which  may  apply 
equally  across  the  image).  The  eye  of  the  observer,  however , 
considers  local  Imagery  changes  as  imagery  detail.  Thus,  those 
system  changes  which  result  in  nonhomogeneous  imagery  degrada- 
tion, i.e.,  localized  changes,  will  be  especially  apparent. 

In  the  case  of  digital  scanners  and  reproducers,  the 
following  is  a list  of  commonly  encountered  practical  departures 
from  theory,  possible  sources  of  the  departure  In  terms  of  com- 
monly encountered  scanners  and  recorders,  and  the  most  likely 
perceived  form  of  the  resulting  degradation. 

Scanner 

The  Prefilter  MTF  Cutoff  Frequency  is  Below  Optimal 

Possible  Sources: 

• System  design 

• The  scanner  optical  path 

• The  scanner  electronic  circuits  and 
f i 1 ters 

Perceived  Degradations: 

• Loss  of  fine  detail,  blurring  (higher 
spatial  frequencies) 

The  Prefilter  MTF  Cutoff  Frequency  is  Above  Optimal 

Possible  Sources: 

t System  design 

• The  scanner  optical  path 
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0 The  scanner  electronic  circuits  and  filters 
Perceived  Degradations: 

0 False  low  spatial  frequency  patterns 
instead  of  and  at  an  angle  to  strong 
high  spatial  frequency  image  detail 
(such  as  waves  at  sea,  plowed  fields, 
shingled  roofs  , etc.  ) 

0 Periodic  dropouts  of  thin  lines  at  a 
slight  angle  to  the  scanning  raster 
(such  as  seeing  a dashed  line  where 
a power  line  is  known  to  be  present). 

0 Regular  structure  or  detail  on  a 

boundary  in  the  image  known  to  be 
short  and  smooth 

0 General  mottling  of  the  image 

The  Prefilter  MTF  Curve  Shape  Departs  from  Optimality 
Possible  Sources: 

0 System  design 

0 Optics 

0 Detector 

0 Preamplifier 

Perceived  Degradations: 

0 Decreased  resolution  (blurring) 

0 Ringing  (appearance  like  ripples 

In  water) 

0 Increased  raster  visibility 

0 Image  too  dark  or  too  light 
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The  Sampling  Rate  is  Reduced  Appreciably 
Possible  Sources: 


t System  design 

• Electrical  component  failure  in  clocking 

circuit 

• Operator  error 


Perceived  Degradations: 


• Loss  of  fine  detail  (above  1/2  the  sampling 
rate ) 

• Image  degenerates  ultimately  to  a matrix  of 
dots  of  varying  intensity 

• The  effects  (above)  caused  by  the  MTF 
going  above  optimal 


The  Sampling  Rate  is  Systematically  or  Randomly 
Perturbed 


Possible  Sources- 


• System  design 

• Electrical  ccmponent  failure  in  clocking 
circuit 

• Mechanical  scanning  imperfections 
Perceived  Degradations: 

• General  mottling  (for  high  spatial 

frequency  perturbations) 

• Image  waviness  or  patterns  in  the  image 

structure  (for  moderate  spatial  frequency 
perturbati ons  ) 

§ Geome tr i c-  d i s tort i on  accompanied  by 
effects  (above)  caused  by  MTF  cutoff 
going  above  and  below  optimal  (for 
low  spatial  frequency  perturbations). 
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T he  Spa  t i al  Sampling  Structure  Departs  from  Optimality 


Possible  Causes: 


• System  design 

• Electrical  component  failure  in  sync 
circuit,  line  start  circuit 


Perceived  Degradations: 


a Effects  (above)  caused  by  MTF  cutoff 

going  above  optimal-general ly  along  one 
direction  in  the  image,  and  simultaneously 
• The  effects  (above)  caused  by  MTF  cutoff 
going  below  opti ma 1 -a  1 ong  some  other 
direction  in  the  image 


The  Digitizer  Levels  are  All  Chosen  Too  Low 


Possible 

Causes : 

• 

System  design 

• 

Preamplifier  gain  increased 

• 

A/D  malfunctioning 

• 

Input  image  too  light 

Pe  rce i ved 

Degradati ons  : 

• 

Total  loss  of  detail  in  highlights 
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The  Digitizer  Levels  are  All  Chosen  Too  High 


Possible  Causes: 

• System  design 

• Preamplifier  gain  decreased 

• A/D  malfunctioning 

• Input  image  too  dark 

Perceived  Degradation: 

• Total  loss  of  detail  in  shadows 


One  or  More  Digitization  Intervals  is  Enlarged 
Possible  Sources : 


• System  design 

• De tec  tor/p  reamp  1 i fi e r gain  curve  departs 

from  optimality 

• A/D  malfunctioning 


Perceived  Degradations: 


• False  contouring  seen  in  detail-less  areas 

• Portions  of  picture  appear  painted 

• Loss  of  low  contrast  detail 

• Apparent  increase  in  image  contrast 


One  or  More  of  the  Pi qitization  Intervals  is  Randomly 
or  Systematically  Perturbed 


Possible  Sources: 


• Mechanical  vibrations 

• Detector  noise 

• Preamplifier  noise 

• A/D  malfunctioning 


Perceived  Degradations: 


• Spot  noise  in  the  image 

• General  mottling 

• Granularity 
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Reproducer 


The  Gray  Level  Mapping  Levels  are  all  Chosen  Too  Low 

Possible  Causes: 

0 System  design 

0 Insufficient  light  output 

0 Insensitive  output  medium 

0 Amplifier  gain  too  low 

Perceived  Degradations 

0 Image  is  too  dark 

0 Cannot  see  shadow  detail 


The  Gray  Level  Mapping  Levels  are  all  Chosen  Too  High 

Possible  Causes: 

0 System  design 

0 Too  intense  a light  output 

0 Too  sensitive  an  output  medium 

0 Amplifier  gain  too  high 

Perceived  Degradations: 

0 Image  is  "washed  out,"  flat 

0 Low  contrast 

One  or  More  of  the  Gray  Level  Mapping  Levels  is 
Mischosen 

Possible  Causes: 


0 System  design 

0 Amplifier/light  modulator/output  medium 

gain  curve  departs  from  optimality 
0 D/A  malfunctioning 


Perceived  Degradations: 


0 False  containing 

0 Masking  of  image  detail 

0 Loss  of  contrast 
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One  or  More  of  the  Gray  Level  Intervals  is  Randomly 


or  Systematically  Perturbed 

Poss i bl e 

Sources : 

• 

Mechanical  vibration 

• 

Light  modulator  noise 

• 

Amplifier  noise 

• 

Channel  noise 

• 

A/D  malfunctioning 

Perce i ved 

Degradati ons  : 

• 

Spot  noise  in  the  image 

• 

General  mottling 

• 

Scan  pattern  becomes  visible 

• 

Granul ari ty 

The  Postfilter  MTF  Cutoff  Frequency  is  Below 

Opti ma 1 

Possible 

Causes : 

• 

System  design 

• 

The  reproducer  optical  path 

• 

The  reproducer  electronic  filters 

Perceived 

Degradati ons  : 

• 

Loss  of  fine  detail,  blurring 

The  Postfilter  MTF  Cutoff  Frequency  is  Above 

Optimal 

Possible 

Sources : 

• 

System  design 

• 

The  reproducer  optical  path 

• 

The  reproducer  electronic  circuits 

Percei ved 

Degradations: 

• 

General  mottling  of  the  image 

• 

The  image  raster  will  become  more 
pronounced 
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• Object  boundaries  known  to  be  smooth 
will  become  ragged 

The  Postfilter  M TF  Curve  Departs  from  Optimality 
Possible  Sources: 

• System  design 

• Optics 

0 Electronics 

• Light  modulator 

Perceived  Degradations: 

• Decreased  resolution  (blurring) 

0 Ringing  (appearance  like  ripples  in  water) 

• Increased  raster  visibility 

• Image  too  dark  or  light 

The  Reproduced  Pixel  Rate  is  Reduced  Appreciably 
Possible  Sources: 

• System  design 

• Electrical  failure  in  clocking  circuit 

• Operator  error 

Perceived  Degradation: 

• Image  will  degenerate  ultimately 
to  a matrix  of  dots  of  varying 
intensity 

• The  effects  (above)  caused  by  the  post- 

filter MTF  going  above  optimal 

The  Reproduced  Pixel  Rate  is  Systematically  or 
Randomly  Perturbed 

Possible  Sources: 

0 System  design 

t Electrical  failure  in  clocking  circuit 

0 Mechanical  scanning  imperfections 
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Perceived  Degradations: 

• General  mottling  (for  high  spatial 
frequency  perturbations) 

• Image  waviness  or  patterns  across  the 
image  structure  (for  moderate  spatial 
frequency  perturbations) 

• Geometric  distortions  caused  by  incorrect 
MTF  of  the  postfilter 

The  Reproduced  Pixel  Structure  Departs  from  Optimality 

Possible  Causes: 

• System  design 

• Electrical  failure  in  sync  circuit  or 

line  start  circuit 

Perceived  Degradations: 

• Systematic  jitter  of  image  detail 

• Incorrect  postfilter  MTF  effects  (above) 
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SECTION  IV 


SUBJECTIVE  EFFECTS  OF  SELECTED  DIGITAL  DEGRADATIONS 

Sections  II  and  III  consider  the  structural  aspects  of  the 
imagery  chain,  methodology  ot  evaluation  and  the  definition  and 
desirable  features  of  a measure  of  image  quality.  These  sec- 
tions do  not,  however,  directly  consider  the  question  of  how 
sensitive  and  consistent  the  human  observer's  reaction  is  to 
various  digital  degradations  likely  to  be  encountered  in 
practice.  This  is  experimentally  considered  here. 

Section  IV-A  discusses  the  test  methodology.  Section  IV-B 
presents  preliminary  image  comparative  test  sets  for  a variety 
of  digital  degradations  and  source  imagery,  and  Section  IV-C 
provides  indications  of  the  sensitivity  of  a trained  photo- 
interpreter to  various  common  digital  degradation  as  well  as 
the  consistence  of  his  observation. 

The  comparative  test  sets  and  corresponding  indications 
of  the  photointerpreter's  sensitivity  to  the  various  digital 
degradations  should  prove  useful  to  a variety  of  hardware 
programs  in  the  future. 

A.  TESTING 

The  testing  procedure  used  was  generally  that  of  Figure 
2- 18 ( b ) , i.e.,  subjective  evaluations  were  sought  of  all  per- 
ceived digital  effects. 

In  order  to  choose  representative  examples  of  the  above 
imagery,  a trained  photointerpreter  intimately  familiar  with 
the  RADC  Imagery  Data  Bank,  Richard  Petroski  of  Rome  Research 
Corporation,  was  asked  to  make  an  appropriate  selection.  The 
following  predigitized  images  were  chosen  for  evaluation. 
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• B/W  Vertical  Photography  (RC-8  camera)  over 

2 

Griffiss  AFB,  NY  - approximately  (1%  ground  ft.) 
per  pixel. 

• B/W  Vertical  Photography  (RC-8  camera)  over  Griffiss 
AFB,  NY  - approximately  (8  ground  ft.)  per  pixel. 

• Infrared  image  (sensor  not  identified)  over  Parked 
Aircraft. 

• Synthetic  Aperature  Radar  (UPD-4)  over  Griffiss  AFB, 

NY. 

Representative  256x512  pixel  samples  of  each  image  were 
chosen  and  used  as  the  test  reference  (Figure  4-1).  Each  image 
was  then  subjected  to  controlled  intensities  of  digital  degra- 
dation in  accordance  with  Table  1.  Appendix  A contains  the 
documentation  of  the  methodology  and  computer  programming 
used.  This  methodology  and  programming  has  the  unique  feature 
of  being  easily  extendable  to  include  other  digital  effects 
including  those  encountered  in  a digital  scanner. 

The  degraded  images  shown  in  Section  IV-B  were  obtained 
in  hardcopy  and  softcopy  form.  The  hardcopy  medium  was  3-M's 
dry  silver  paper.  The  softcopy  display  was  the  Sony  Color 
Digital  TV  at  RADC. 

A trained  photointerpreter,  again  Richard  Petroski  of 
Rome  Research  Corporation,  was  then  given  the  directions  con- 
tained in  Table  2.  The  photoi nterpreter  was  permitted  to 
select  his  working  environment  and  use  whatever  aids  he  desired. 
The  results  of  the  above  testing  were  collected  and  compiled. 
Results  are  given  in  Section  IV-C  and  Appendix  B. 

B.  PRELIMINARY  IMAGE  COMPARATIVE  TEST  SETS 

The  figures  of  this  section  contain  high  quality  photo- 
graphic reproducti ons  of  the  series  of  hardcopy  degraded 
imagery  used  in  the  experimental  testing  of  this  section. 
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£ B/W  vertical  photography  (RC-8 
camera)  over  Griffiss  AFB,  NY  ■ 
approxi mately  (1-1/4  ground 
ft.)2  per  pixel. 


3/W  vertical  photography  (RC-8 
camera)  over  Griffiss  AFB,  NY  - 
approximately  (8-3/4  ground 
f t . ) 2 per  pixel. 


Infrared  image  (sensor  not 
identified)  over  parked 
aircraft 


[^T  ?/ 
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Synthetic  Aperture  Radar 
(UPD-4)  over  Griffis  AFB,  NY 


X . 


Figure  4-1.  Master  Digital  Images 
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TABLE  1.  LIST  OF  DIGITAL  DEGRADATION  & VALUES  CONSIDERED 


Additive  Gaussian  Spot  Noise 

Mean  = 0 

Standard  Deviation  = .3,  .6,  1.22,  2.44, 
4.88,  9.77,  19.52,  39.04%  of  the  gray 
level  dynamic  range 

Additive  Gaussian  Line  Noise 

Mean  = 0 

Standard  Deviation  = .3,  .6,  1.22,  2.44, 
4.88,  9.77,  19.52,  39.04%  of  the  gray 
i level  dynamic  range 

Gaussian  Line  Jitter 

Mean  = 0 

Standard  Deviation  = .05,  .1,  .2,  .39,  .78, 
1.56,  3.13,  6.25  pixels 

Quanti zation 

No.  of  levels  = 32,  16,  8,  4 gray  levels 

Sampling  Reduction,  High  Aliasing 

Reduction  Size  = 1/2,  1/3,  1/4,  1/5  of 
original 

Sampling  Reduction,  Reduced  Aliasing 

Reduction  Size  = 1/2,  1/3,  1/4,  1/5  of 
original 
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TABLE  2.  PHOTOINTERPRETATOR  TEST 
INSTRUCTIONS 


( 1 ) Subjective  Verbal  Assessment  of  Image  Quality 

A.  Look  at  each  digitally  perturbed  image 
in  reference  to  the  original  and  define 
its  percent  utility  to  photoi nterpreti ve 
tasks . 

100%  utility  - same  as  original 
0%  utility  - of  no  use  whatsoever 

B.  Note  in  words  the  nature  and  extent  of 
interference  to  image  quality  relating 
to  your  percentile  utility  rating 
assi gned . 


( 2 ) Determination  of  Isopreference  Contours  of  Image 

Quality 

Using  selected  images  from  one  digitally 
perturbed  series,  match  them  to  images  of  other 
series  such  that  they  visually  seem  to  possess 
the  same  image  quality. 


( 3 ) Crosscheck  of  the  Transitiveness  of  Image  Quality 

Repeat  Number  (2)  above. 


( 4 ) Specific  Comparison  of  Aliased  Imagery 

For  each  image  pair  of  digitally  reduced 
spatially  sampled  images  (high  aliasing,  low 
aliasing)  determine  which  possess  the  higher 
image  quality. 
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They  are  included  as  a key  component  of  this  report  since  via 
comparative  techniques,  they  may  be  used  to  evaluate  actual 
digital  systems  for  similar  degradations. 

Figures  4-2  through  4-5  illustrate  the  effects  of  pro- 
gressive digital  degradation  on  a high  resolution  B/W  vertical 
photograph,  the  top  image  of  Figure  4-1. 

Figures  4-6  through  4-11  illustrate  the  effects  of  pro- 
gressive digital  degradations  on  a low  resolution  B/W  vertical 
photograph,  the  second  image  of  Figure  4-1. 

Figures  4-12  through  4-17  illustrate  the  effects  of  pro- 
gressive digital  degradations  on  an  infrared  image  of  parked 
aircraft,  the  third  image  of  Figure  4-1. 

Finally  Figures  4-18  through  4-23  illustrate  the  effects 
of  progressive  digital  degradations  on  a synthetic  aperature 
radar  image,  the  bottom  image  of  Figure  4-1. 

C.  RESULTS 

The  results  of  subjective  photoi nterpretati ve  evaluation 
of  the  preliminary  test  set  of  Section  IV-Baregiven  here. 

The  preliminary  subjective  image  test  set  was  evaluated 
according  to  the  procedure  of  Section  IV-A  and  Table  3.  The 
results  of  the  subjective  hardcopy  utility  assessment  of  photo- 
interpretative  image  content  relative  to  the  origina.l  digitized 
images  are  shown  in  Figures  4-25  through  4-29.  The  various 
absolute  sensitivities  recorded  should  be  useful  in  comparative 
evaluations  and  the  design  of  other  imagery  systems. 

The  figures,  4-24  through  4-29  show  that  the  subjective 
image  quality  measure  is  a function  of  the  form  and  degree 
of  the  digital  degradation.  The  subjective  image  quality  is 
also  seen  to  be  dependent  upon  the  characteristics  of  the  source 
image.  The  Interdependency  of  the  form  and  degree  of  digital 
degradation  with  the  characteristics  of  the  source  image  are 
found  to  be  most  complex  and  deserving  of  further  study.  No 
immediate  correlation  between  the  image  spectrum  of  the  original 
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Standard  deviation 
of  gray  level  range 


Standard  deviation 
of  gray  level  range 


Standard  deviation 
of  gray  level  range 


Standard  deviation 
of  gray  level  range 


figure  4-2.  Additive  Gaussian  Spot  Noise 
Applied  to  a High  Resolution  B/W  Vertical 
Digitized  Photograph 


V * 


Standard  deviation  = 4.88% 
of  gray  level  range 


Standard  deviation  = 9.77% 
of  gray  level  range 
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Standard  deviation  = 19.55% 
of  gray  level  range 


Standard  deviation  = 39.04% 
of  gray  level  range 


Figure  4-2  (Continued).  Additive  Gaussian  Spot  Noise 
Applied  to  a High  Resolution  B/W  Vertical 
Digitized  Photograph 


Standard  deviation  = .3 
of  gray  level  range 


tandard  deviation  = .6 
f gray  level  range 


tandard  deviation  = 1 
f gray  level  range 


Standard  deviation  = 2 
of  gray  ldvel  range 


Figure  4-3.  Additive  Gaussian  Line  Noise  Applied 
to  a High  Resolution  B/W  Vertical  Digitized 
Photograph 


Standa  rd  deviation  = 
of  gray  level  range 


Standard  deviation  = 
of  gray  level  range 


Standard  deviation  = 
of  gray  level  range 


Standard  deviation  = 
of  gray  level  range 


4.88% 


9.77% 


19.52% 


39.04% 
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gure  4-3 
Appl  i 


(Continued).  Additive  Gaussian  Lire 
ed  to  a High  Resolution  B/W  Vertical 
Digitized  Photograph 


Noise 
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Sampling  Reduction  = 2 


Sampling  Reduction  = 3 


Sampling  Reduction  * 4 


Sampling  Reduction  = 5 


Figure  4-5.  Varying  the  Sampling  of  a High  Resolution  B/W 
Vertical  Digitized  Photograph  Under  Reduced  Aliasing 
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Standard  deviation  = .3% 
of  gray  level  range 


Standard  deviation  = .6% 
of  gray  level  range 


Standard  deviation  = 1.2% 
of  gray  level  range 


Standard  deviation  = 2.44% 
of  gray  level  range 


Figure  4-6.  Additive  Gaussian  Spot  Noise  Applied 
to  a Low  Resolution  B/W  Vertical  Digitized  Photograph 
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Standard  deviation  = 
of  gray  level  range 


Standard  deviation  = 
of  gray  level  range 


Standard  deviation  = 
of  gray  level  range 


Standard  deviation 
of  gray  level  range 


Figure  4-6  (Continued).  Additive 
Applied  to  a Low  Resolution  R 4 . • 

Photograph 
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Standard  deviation  = .3% 
of  gray  level  range 


Standard  deviation  - .6% 
of  gray  level  range 


Standard  deviation  * 1.22* 
of  gray  level  range 


Standard  deviation  ■ 2.44% 
of  gray  level  range 


Figure  4-7.  Additive  Gaussian  Line  Noise  Applied 
to  a Low  Resolution  B/W  Vertical  Digitized  Photograph 
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Standard  deviation  = 
of  gray  level  range 


Standard  deviation  = 
of  gray  level  range 


Standard  deviation  ■ 
of  gray  level  range 


Standard  deviation  * 
of  gray  level  range 


4.88* 


9.77% 


19.52% 


39.04% 


Figure  4-7  (Continued).  Additive  Gaussian  Line  Noise  Applied 
to  a Low  Resolution  B/W  Vertical  Digitized  Photograph 


Standard  deviation 
pixels 


Standard  deviation 
pixels 


Standard  deviation 
pixels 


Standard  deviation 
pixels 


= .05 


= .2 


= .39 


Figure  4-8.  Gaussian  Line  Jitter  Applied  to  a 
Low  Resolution  B/W  Vertical  Digitized 
Photograph 
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Standard  deviation  = 6.25 
pixels 


Figure  4-8  (Continued).-  Gaussian  Line  Jitter  Applied 
to  a Low  Resolution  B/W  Vertical  Digitized 
Photograph 
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Number  of  gray  levels  ■ 32 


Number  of  gray  levels  = 16 


Number  of  gray  levels  * 8 


Number  of  gray  levels  * 4 


Figure  4-9.  Varying  the  Quantization  of  a Low 
Resolution  B/W  Vertical  Digital  Photograph 


Sampling  Reduction 


Resol utlon 


Figure  4-10. 
B/W  Vertical 


Varying  the  Sampling  of  a Low 
Digitized  Photograph  Under  High  Aliasing 
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Sampling  Reduction  ■ 2 
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Sampling  Reduction  » 3 
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Sampling  Reduction  - 4 


Sampling  Reduction  • 5 


Figure  4-11.  Varying  the  Sampling  of  a Low  Resolution 
B/W  Vertical  Digitized  Photograph  Under  Reduced  Aliasing 
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Standard  deviation  = .6 
of  gray  level  range 


Standard  deviation  = 1. 
of  gray  level  range 


Standard  deviation  = 2. 
of  gray  level  range 


Figure  4-12.  Additive  Gaussian  Spot  Noise  Applied 
to  an  Infrared  Image  of  Parked  Aircraft 


Standard  deviation  = .3 
of  gray  level  range 
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Standard  deviation  - 
of  gray  level  range 


Standard  deviation  ■ 
of  gray  level  range 


Standard  deviation  * 
of  gray  level  range 


Standard  deviation  ■ 
of  gray  level  range 
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9.77X 


19 . 52% 


39.04X 


Standard  deviation  = .3% 
of  gray  level  range 


Standard  deviation  = .6% 
of  gray  level  range 


Standard  deviation  = 1.22% 
of  gray  level  range 


Standard  deviation  = 2.44% 
of  gray  level  range 


Figure  4-13.  Additive  Gaussian  Line  Noise  Applied  to 
an  Infrared  Image  of  Parked  Aircraft 
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Standard  deviation 
pixels 


Standard  deviation 
pixels 


Standard  deviation 
pixels 


Standard  deviation 
pixels 


Figure  4-14.  Gaussian  Line  Jitter  applied  to  an  Infrared 
Image  of  Parked  Aircraft 
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Standard  deviation 
pixels 


Standard  deviation  = 1.56 
pixels 


Standard  deviation  = 3.13 
pixels 


Standard  deviation  = 6.25 
pixels 


Figure  4-14  (Continued).  Gaussian  Line  Jitter  Applied 
to  an  Infrared  Image  of  Parked  Aircraft 
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Figure  4-15.  Varying  the  Quantization  of  an  Infrared 
Image  of  Parked  Aircraft 
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Sampling  Reduction  = 2 


Sampling  Reduction  = 3 


Sampling  Reduction  = 4 


Sampling  Reduction  * 5 


Figure  4-16.  Varying  the  Sampling  of  an  Infrared 
Image  of  Parked  Aircraft  Under  High  Aliasing 
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Sampling  Reduction  = 2 


Sampling  Reduction  = 3 


Sampling  Reduction  * 4 


Sampling  Reduction 


Figure  4-17.  Varying  the  Sampling  of  an  Infrared 
Image  of  Parked  Aircraft  Under  Reduced  Aliasing 


Standard  devtatlon  = .3% 
of  gray  level  range 


Standard  deviation  = .6% 
of  gray  level  range 


Standard  deviation  = 1.22% 
of  gray  level  range 


Standard  deviation  * 2.44% 
of  gray  level  range 


Figure  4-19.  Additive  Gaussian  Line  Noise  Applied 
to  a Synthetic  Aperture  Radar  fmage 


Standard  deviation  * 4.88% 
of  gray  level  range 


Standard  deviation  ■ 9.77% 
of  gray  level  range 


Standard  deviation  = 19.52% 
of  gray  level  range 


Standard  deviation  * 39.04% 
of  gray  level  range 


Figure  4-19  (Continued).  Additive  Gaussian  Line  Noise 
Applied  to  a Synthetic  Aperture  Radar  Image 
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Standard  deviation  * .78 
pixels 


Standard  deviation  = 1.56 
pixels 


Standard  deviation  = 3.13 
pixels 


Standard  deviation  * 6.25 
pixels 


Figure 


4-20  (Continued).  Gaussian  Line  Jitter  Applied 
to  a Synthetic  Aperture  Radar  Image 
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Figure  4-21.  Varying  the  Quantization  of  a Synthetic 
Aperture  Radar  Image 
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Sampling  Reduction  = 2 


Sampling  Reduction  = 3 


Sampling  Reduction  = 4 


Sampling  Reduction  = 5 


Figure  4-22.  Varying  the  Sampling  of  a Synthetic 
Aperture  Radar  Image  Under  High  Aliasing 


116 


Sampling  Reduction  * 2 


Sampling  Reduction  = 3 


Sampling  Reduction  * 4 


Sampling  Reduction  * 5 


Figure  4-23.  Varying  the  Sampling  of  a Synthetic 
Aperture  Radar  Image  Under  Reduced  Aliasing 
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X UTILITY  % UTILITY 


Figure  4-24.  Effect  of  Additive  Gaussian  Spot  Noise 


Figure  4-25.  Effect  of  Additive  Gaussian  Line  Noise 

NOTE:  the  units  of  standard  deviations  correspond  to  the  actual 

values  chosen  for  the  tests  and  provide  an  approximate 
logrithmic  scale. 
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STANDARD  DEVIATION  OF  GAUSSIAN  JITTER  IN  PIXEL  WIDTHS 


Figure  4-26.  Effects  of  Jitter 


NUMBER  OF  GRAY  LEVELS 


Figure  4-27-  Effects  of  Number  Gray  Levels 
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Figure  4-28-  Effect  of  Sampling  Rate,  High  Aliasing 


Vert.  Photograph 
) Resolution 
Vert.  Photograph 


NOTE: 


Figure  4-29.  Effect  of  Sampling  Rate,  Low  Aliasing 

The  Photointerpreter  was  allowed  use  of  a magnifying 
glass  to  compensate  for  the  change  in  scale  in  these 
sampling  rate  tests. 
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TABLE  3. 


IMAGE  EVALUATION  FACTORS 


1.  Shape  - Shape  relates  to  the  general  configuration 
of  an  object  depicted  on  the  image  scene.  Although 
the  shape  of  an  object  as  seen  from  above  is  some- 
times difficult  to  interpret,  it  often  provides 

the  most  important  clue  to  identification.  Shape 
is  often  distorted  on  radar  and  infrared  imagery 
because  of  such  factors  as  spot  size,  pulse  length, 
beamwidth  and  fold-over  (radar),  and  temperature, 
spot  size,  and  field  of  view  (infrared).  The  air- 
craft attitude  may  introduce  errors  by  pitch,  roll, 
yaw,  and  improper  altitude  and  velocity.  These 
errors  are  common  on  all  three  general  imagery 
types  (photo,  infrared,  and  radar).  Distortions 
on  the  image  scene  may  also  be  caused  by  digitizing 
the  original  data.  When  evaluating  an  image  scene, 
the  interpreter  must  consider  shape  distortion. 

2.  Size  - Size  relates  to  the  dimensions,  surface,  and 
vol ume  of  an  object  depicted  on  the  image  scene. 

The  size  of  one  object  can  often  be  determined  from 
its  relative  size  to  other  objects.  Where  accurate 
dimensions  are  required,  the  scale  of  the  digitized 
image  must  be  calculated.  The  closeness  of  the  ob- 
ject size  on  the  digitized  image  to  the  actual 
object  size  will  be  the  second  factor  considered 
when  ranking  an  image  scene  by  quality. 

3.  Surroundi ngs  - Surroundings  refer  to  the  spatial 
arrangement  of  objects  and  their  relationship  to 
natural  objects  (background).  For  example,  a ther- 
mal electric  power  plant  will  have  various  target 
components  (objects)  surrounding  it  (such  as  a 
transformer  yard,  fuel  supply,  power  line  cuts) 
which  will  enable  the  interpreter  to  correctly 
identify  the  target  by  its  surroundings.  The  loca- 
tion of  the  target  also  falls  Into  the  category  of 
surroundings.  When  an  image  is  degraded,  certain 
of  the  target  surroundings  may  not  be  identifiable. 
This  will  be  the  third  criterion  used  when  ranking 
an  image  scene.  This  category  can  also  be  referred 
to  as  the  "minimum  image  concept"  or  "target  skele- 
ton concept."  Here  the  interpreter  must  decide  the 
mi n i mum  number  of  objects  (target  and  surroundings) 
that  must  be  identifiable  to  correctly  determine  the 
type  of  target  present  on  the  image  scene. 
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TABLE  3.  IMAGE  EVALUATION  FACTORS 
(Continued ) 


4.  Shade  - Shade  or  tone  refers  to  the  brilliance  with 
whi ch  light  is  reflected  by  an  object.  Without  tone 
differences  the  shapes  of  objects  could  not  be  dis- 
cerned on  imagery.  The  overall  image  shade/tone 
quality  will  be  determined  during  this  fourth  step 
of  the  evaluation. 

5.  Shadow  - Shadow  describes  the  condition  wherein  an 
intervening  object  prevents  the  direct  sun  rays  from 
striking  certain  areas  on  the  imagery.  Shadows  are 
important  because  their  shapes  are  indicative  of 

the  profile  views  of  objects,  and  can  facilitate 
recognition,  as  in  the  case  of  a church  steeple 
with  the  cross  on  the  top,  or  a tall  tower. 


images.  Figures  2-6  through  2-8,  could  be  found  with  the  curves 
of  performance  against  degradation.  This  is  thought  to  be  due 
to  the  fact  that  the  photointerpreter  tends  to  concentrate  on 
targets  while  the  spectrum  data  is  composed  largely  of  back- 
ground data. 

Figures  4-30  through  4-33  show  comparisons  of  subjective 
image  utility  as  a function  of  display  type.  The  results 
indicate,  as  expected,  that  hardcopy  and  softcopy  utility  is 
roughly  comparable  when  the  same  image  detail  degradation  and 
extent  are  shown  on  both  outputs.  There  appears  to  be  a slight 
bias  toward  preference  of  softcopy  over  hardcopy  products.  This 
bias  may  be  a function  of  the  testing  procedure  although  the 
softcopy  evaluations  were  held  more  than  30  days  after  the 
hardcopy  evaluations  due  to  the  nonavailability  of  a softcopy 
display  at  Harris  ESD  during  this  time  frame.  The  apparent 
preference  of  softcopy  over  hardcopy  is  probably  due  to  the 
fact  that  both  the  softcopy  and  hardcopy  curves  for  no  noise 
start  at  100%  utility.  This  is  probably  not  true:  the  basic 
utility  of  the  noiseless  hardcopy  is  probably  better  than  the 
basic  utility  of  the  noiseless  softcopy.  Furthermore,  as  noise 
is  introduced,  the  softcopy  would  be  affected  less  than  the 
hardcopy  because  the  softcopy  originally  had  more  noise  in  the 
form  of  scan  lines,  phosphor  noise  on  the  screen,  etc. 

Other  tests  were  held  to  indicate  the  repeatability  and 
transitiveness  (see  Section  1 1 - G ) of  the  photointerpreter's 
utility  assessment.  The  repeatability  was  on  the  order  of  80- 
90%.  The  utility  assessment  did  seem  to  be  transitive,  thus  the 
subjective  utility  assessment  is  an  equivalence  relation.  As  a 
result,  it  appears  reasonable  to  talk  of  equal  quality  of  two 
or  more  images  even  though  they  have  been  degraded  by  quite 
different  mechanisms. 
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V 


Figure  4-30  . Effect  of  Additive  Gaussian  Spot  Noise 
(Hardcopy  vs.  Softcopy) 
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Figure  4-31.  Effect  of  Additive  Gaussian  Line  Noise 
(Hardcopy  vs.  Softcopy) 


124 


% UTILITY  X UTILITY 


Figure  4-32.  Effects  of  Jitter 
(Hardcopy  vs.  Softcopy) 
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Figure  4-33.  Effects  of  Number  Gray  Levels 
(Hardcopy  vs.  Softcopy) 
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SECTION  V 

POSSIBLE  DIGITAL  IMAGE  SYSTEM  OBJECTIVE 
TEST  PATTERNS 

A set  of  test  patterns  will  now  be  given  which  can  be  used 
to  measure  the  various  types  of  system  degradations.  Four  basic 
types  of  image  quality  degradations  of  the  imagery  processing 
subsystem  are  proposed  as  criteria  for  the  Test  Set:  (a) 
amplitude  distortion,  i.e.,  incorrectly  representing  the  gray 
shade  at  a particular  point  in  the  image,  (b)  the  system  induced 
modulation  transfer  function  (MTF),  (c)  geometric  distortion, 
i.e.,  the  "warping"  of  the  output  scene  due  to  system  geometric 
non-linearities,  and  (d)  noise,  which  includes  channel  noise 
(bit  errors),  quantization  noise,  aliasing  noise,  etc. 

This  section  provides  simple  visual  methods  of  measuring 
the  first  three  types  of  image  degradations  mentioned  above. 

A.  MEASUREMENT  OF  SAMPLING  RATE 

The  sampling  rates  of  a digital  imagery  system  can  easily 
be  measured  through  use  of  Moire'  patterns.  The  following 
basis  for  a test  procedure  to  determine  the  reproduction  sam- 
pling rate  is  illustrative: 

(1)  Scan  or  reproduce  a test  image.  Figure  5-1, 
consisting  of  vertical  lines,  with  a density 
comparable  to  the  scan  density  to  be  measured, 
and  obtain  a reproduced  output. 

(2)  Lay  the  comparative  test  image  consisting  of 
radial  lines  emerging  from  the  center,  Figure 
5-2,  over  the  reproduced  test  image  obtained 
from  the  above  step. 
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(3)  The  radial  distance  to  the  most  prominent 
Moire'  pattern  is  inversely  proportional 
to  the  spatial  sampling  frequency  of  the 
digital  system. 

Example  results  are  shown  in  Figure  5-3  where  the  distance 
c is  inversely  proportional  to  the  sampling  frequency  of  the 
system  under  test. 

The  reasoning  behind  the  test  procedure  can  be  deduced 
from  the  details  provided  in  Appendix  C. 

The  following  questions  and  considerations  are  among  those 
which  should  be  answered  to  reduce  the  above  procedure  to  actual 
practice. 

(a)  Two  particularly  interesting  and  important 
things  to  note  are  that  the  figure  of  Step  1 
is  input  into  the  digital  system  and  not  the 
figure  of  Step  2,  and  that  the  figure  of 
Step  1 has  a pattern  frequency  of  one-half 
the  digital  spatial  sampling  frequency.  See 
Appendix  C for  the  necessary  expanation. 

(b)  The  test  pattern  of  the  figure  of  Step  1 is 
specified  to  have  a basic  spatial  frequency 
one-half  the  sampling  rate  of  the  system  to 
be  tested.  If  the  digital  image  reproducer 
only  is  to  be  tested,  how  can  this  be  accom- 
plished? If  a complete  digital  system 
Including  a digital  image  scanner  is  to  be 
tested,  it  is  generally  impossible  to  meet 
the  above  specification.  How  should  the 
test  procedure  be  altered? 











III  Mil  II  III  II I Hill  II  lllll  II  HOI  li  tiiitmmti  

SSSSK-ss^SW 



; 



> 


'!i" 


\llu'S 

I1 

7 


THE  MOST 
PROMINENT 
MOIRE  PATTERN 





•<  /'i  Mi’i 1 


1 1 1 1 

MM 

f pr 


m 

i ■ 


in 


■H 


3 


(c)  How  can  the  sampling  frequency  be  indicated  on 
the  comparison  image  of  the  figure  of  Step  2? 

A partial  answer  is  to  use  concentric  circles 
indicating  the  spatial  frequency  corresponding 
to  the  most  significant  Moire'  pattern. 

(d)  Since  the  sampling  can  occur  along  more  than 
one  spatial  axis  s i mu  1 ta neous ly  , how  does 
the  technique  given  above  apply? 

(e)  How  can  the  above  procedure  be  altered  to 
determine  the  scanner  sampling  rate? 

Another,  very  precise,  determination  of  the  spatial  sampling  rate 
in  lieu  of  other  sampling  irregularities  follows: 

(1)  Input  the  test  pattern  of  Figure  5-1  into 
the  system. 

(2)  Overlay  the  reproduced  test  pattern  with 
the  original  test  pattern  canted  at  a 
slight  angle. 

(3)  The  spacing  between  the  Moire'  patterns 
seen  is  proportional  to  the  spatial 
reproduction  rate  of  the  system  under 
test. 

An  illustrative  example  of  the  Moire'  patterns  to  be  ob- 
tained is  shown  in  Figures  5-4  and  5-5.  The  slight  bending  of 
the  patterns  is  a result  of  slight  variations  in  the  reproduc- 
tion rate  aero  s the  image. 

B.  MEASUREMENT  OF  MTF 

MTF  is  a powerful  and  increasingly  used  guide  to  image 
quality  for  the  following  reasons: 

• The  MTF  of  multiply  cascaded  imagery  system 
components  can  be  easily  combined  into  a 
system  MTF. 
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Figure  5-4.  Determination  of  Spatial 
Sampling  Reproduction  Rate 
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Figure  5-5.  Precise  Determi na t i on  of  Spatial 
Sampling  Reproduction  Ratf? 
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• The  MTF  can  be  converted  into  the  "noise- 
equivalent  passband"  of  the  system,  a 
parameter  of  particular  interest  in  analysis 
of  communications  systems. 

t The  MTF  is  a guide  to  the  capability  of  the 

system  to  image  various  sized  imagery  detail. 

Considered  here  is  a technique  by  which  the  system  MTF  can  be 
measured.  The  technique  is  based  upon  a chirped  spatial 
frequency  procedure: 

(1)  Input  a test  image  containing  the  luminance 
function  given  by 

I ( x ,y  ) = 1 + y sinx2  (1) 

where  the  coordinate  system  is  so  chosen  that 
x corresponds  with  the  horizontal  frequency 
axis  and  y the  vertical  contrast  axis. 

(2)  Take  the  resulting  reproduced  test  image  (an 
example  is  shown  in  Figure  5-6)  and  lay  it 
down  beside  the  original  test  image. 

(3)  Back  up  from  both  test  images  until  it 
becomes  possible  to  draw  a curve  in  each 
separating  the  regions  containing  percep- 
tual sine  waves  from  those  not  perceptual. 

Call  the  curves  a(x)  and  b(x)  as  shown  in 
Figure  5-7. 

(4)  Derive  the  system  MTF  by  plotting 

MTF<*>  ’ -rtef- 

with  results  as  shown  in  Figure  5-8. 

The  reasoning  behind  the  procedure  follows  in  Appendix  E. 
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Figure  5-6 


A Reproduced  Chirped  Test  Image 
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The  following  questions  should  be  answered  prior  to 
reducing  the  above  procedure  to  practice: 

(1)  Can  the  system  truly  be  considered  linear? 

If  not,  how  does  the  concept  of  MTF  apply? 

If  so,  what  are  error  sources?  See  also 
Question  2. 

(2)  A chirped  sine  wave  is  not  the  same  as  a 
sine  wave  (as  required  by  the  definition 
of  MTF).  How  does  this  affect  accuracy 
and  how  can  any  errors  be  reduced?  Look 
at  variations  of  chirp  rate! 

(3)  Why  does  the  procedure  attempt  to  measure 
MTF  only  up  to  one  half  the  sampling  fre- 
quency of  the  system  under  test. 

(4)  Use  of  the  perception  threshold  characteris- 
tics of  the  eye  could  be  faulted  as  relying  on 
a marginal  and  highly  variable  phenomena . 

How  does  the  procedure  described  greatly  reduce 
any  possible  deleterious  effects? 

(5)  What  is  a great  advantage  of  the  use  of  the 
perception  threshold  in  the  procedure?  Con- 
sider that  most  nonlinearities  approach 
linearity  in  the  small  signal  case!  Does  the 
procedure  at  all  times  minimize  the  amplitude  of 
the  signal?  What  digital  non  1 i nea ri ti es  may  not 
approach  linearity  in  the  small  signal  case? 

What  Is  the  effect  of  digitization  upon  the 
procedure  and  indeed  on  any  procedure  designed 
to  measure  resolution  of  MTF? 
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(6)  How  could  the  division  indicated  in  Step  5 

be  simplified  to  a manual  rather  than  mental 
exercise?  Consider  a test  signal  in  Equation 
(1) 

2 

I ( x ,y ) = 1 + log(y)  sin(x) 

C.  MEASUREMENT  OF  SAMPLING  INHOMOGENEITIES 

A common  problem  in  digital  sampling  is  the  fact  that 
the  sampling  rate  often  varies  across  the  image  with  attendant 
effects.  These  inhomogeneities  can  easily  be  measured,  in  lieu 
of  other  sampling  degradations,  by  the  following  procedure: 

(1)  Scan  or  reproduce  the  test  Image,  Figure  5-9, 
consisting  of  horizontal  lines  and  obtain  a 
reproduced  test  image. 

(2)  Lay  a slightly  enlarged  transparent  copy  of 
the  test  image,  of  Step  1,  over  the  reproduced 
test  image  obtained  above,  canted  at  a slight 
angle. 

(3)  The  Moire'  pattern  observed  is  a good  approxi- 
mation to  the  curve  of  the  inhomogeneity  of  the 
scan  rate. 

As  an  example  of  the  above  procedure  let  Figure  5-10  be 
the  reproduced  image.  Note  that  it  has  a scan  nonhomogeneity 
which  while  present  Is  difficult  to  fully  quantify.  Laying  a 
slightly  enlarged  canted  copy  of  the  figure  from  Step  1 over 
Figure  5-10  gives  the  Figure  5-11. 

Since  the  theory  behind  the  appearance  of  this  particular 
Moire'  pattern  is  especially  Intriguing,  it  is  further  con- 
sidered in  Appendix  0.  It  becomes  apparent  there,  for  instance, 
that  theoretically  (and  in  practice  to  a large  degree)  any 
degree  of  magnification  of  scanning  Inhomogeneities  is 
possible. 
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Figure  5-11.  Example  of  Moire'  Scan  Nonlinearity 
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The  following  questions  are  among  those  which  should  be 
answered  to  reduce  the  above  procedure  to  actual  practice. 

(1)  How  do  other  forms  of  digital  sampling 
degradations  affect  the  accuracy  of  the 
test? 

(2)  What  are  the  second  order  effects  of  the 
approxi mati on  ? 

(3)  What  does  "slightly  enlarged"  mean  in  the 
2nd  step  of  the  procedure  (see  Section  V-E)? 

(4)  What  is  the  effect  of  changing  the  angle  a in 
Figure  5-11?  (The  first  order  effect  is  simple 
to  obtain.) 

D.  MEASUREMENT  OF  REPRODUCER  SAMPLING  STRUCTURE 

Some  indication  of  the  spot  shape  and  placement  structure 
of  the  reproducer  can  be  obtained  through  the  following  pro- 
cedure in  lieu  of  a microscope: 

(1)  Reproduce  the  mathematical  test  pattern  of 
Figure  5-12. 

(2)  Overlay  the  reproduced  test  pattern  with  a 
transparent  copy  of  the  test  pattern  at 
slight  cant. 

(3)  The  faint  checkered  pattern  is  a greatly 
amplified  copy  of  the  output  image. 

Figure  5-13  is  a representative  example  of  results  possible 
with  this  method.  In  this  case  the  Moire'  pattern  structure 
shows  the  microstructure  of  a 50%  half-tone  image. 
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Figure  5-12.  Mathematical  Spot  Test  Image 
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Figure  5-13.  Example  Results  of  Structural  De t jrmi na t i on 
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E.  OBSERVATION  OF  SYSTEMATIC  AND  RANDOM  IMAGE  PERTURBATIONS 

A very  common  form  of  digital  degradation  is  the  systematic 
or  random  perturbation  of  the  image  resulting  in  various  super- 
imposed image  textures  and  patterns.  While  it  is  conceivable 
that  these  perturbations  are  functions  of  image  content,  the 
much  more  usual  case  is  that  the  perturbation  will  be  independent 
of  the  image  detail  and  brightness.  The  most  sensitive  test  for 
this  case  will  be  to  input  a flat  (dc,  no  detail,  gray)  image 
of  appreciable  extent.  Any  observable  detail  within  the  output 
image  will  be  due  to  degradation. 

F.  DETERMINATION  OF  THE  IMAGERY  SYSTEM  GRAY-LEVEL  MAPPING 

AND  NUMBER  OF  LEVELS  OF  QUANTIZATION 

A test  image  consisting  of  a continuous  gray-level  wedge 
is  input  to  the  scanner.  The  gray-levels  of  the  output  wedge 
in  the  reproduced  image  can  be  visually  compared  to  the  ori- 
ginal image  gray-level  wedge  in  order  to  draw  a curve  of  output 
to  input  gray  levels.  Ideally  for  maximum  gray-level  realism, 
this  curve  should  be  a straight  45°  positive  slope  line. 

If  the  gray-level  of  the  reproduced  image  wedge  is  com- 
posed of  easily  discernible  gray-level  steps,  then  the  number 
of  observable  steps  is  equal  to  the  number  of  observable  quan- 
tization levels  and  of  the  potential  for  false  contouring 
within  an  image. 

A quick  test  for  gray-level  linearity  of  an  image  system 
Is  as  foil ows  : 

(1)  Input  a gray-level  wedge  into  the  scanner 

(2)  Overlay  the  output  gray-levei  wedge  from 
the  reproducer  with  a transparent  copy  of 
the  original  gray-level  wedge.  Turn  this 
copy  over  as  required  to  get  an  inverse  gray- 
level  wedge,  i.e.,  decreasing  density  rather 
than  increasing  density. 
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The  resulting  image  area  of  the  reproduced 
image  overlaid  with  an  inverted  transparent 
copy  of  the  original  should  be  uniformly 
gray  under  strong  light.  Any  deviation  from 
uniformity  will  be  an  indication  of  system 
gray-level  non  1 i nea ri ty . 
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SECTION  VI 

RESULTS  AND  CONCLUSIONS 

This  study  has  taken  an  in-depth  look  at  the  problem  of 
defining  and  measuring  the  image  ouality  of  a digital  image 
processing  system.  We  have  taken  two  approaches  in  an  attempt 
to  measure  image  quality.  The  subjective  method  whereby  a 
human  observer  judged  the  quality  of  degraded  imagery  in  a 
subjective  manner,  and  the  objective  method  whereby  mathe- 
matical or  physical  attributes  of  the  imagery  system  were 
measured  in  an  objective  manner. 

The  subjective  testing  of  imagery  involved  adding  various 
types  of  commonly  observed  noise  or  degradations  to  four  dif- 
ferent types  of  images.  The  subjective  "utility"  of  the 
degraded  imagery  was  observed  and  tabulated.  These  results 
will  enable  a system  designer  to  predict  the  performance  of 
an  imagery  system  for  a number  of  common  tasks  when  the  system 
has  been  corrupted  with  various  types  of  digital  noise  or 
degradations. 

The  objective  method  of  measuring  image  degradations 
involved  the  conceptual  measurement  of  the  following  type 
degradations : 

• amplitude  distortion, 

• modulation  transfer  function, 

• geometric  distortion 

The  fourth  type  of  degradation,  classified  as  noise,  was  dis- 
cussed In  this  report  although  no  simplified  objective  test 
was  included.  The  aliasing  noise,  discussed  in  Section  1 1 1 - A 
was  exemplified  by  the  chirp  MTF  test  of  Figure  5-6  although 
the  quantitative  measure  was  proposed  for  the  MTF  not  aliasing 
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distortion.  One  of  the  recommendations  for  further  work  is 
the  definition  of  the  trade-offs  in  performance  between  the 
conflicting  parameters  of  higher  resolution  and  aliasing  error 
i.e.,  if  the  resolution  of  a digital  imaging  system  is  made 
to  increase,  then  the  aliasing  error  will  also  increase. 

A model  of  imagery  was  developed  using  the  Poisson 
distribution  for  the  location  of  edges  in  imagery.  This  model 
generated  a spectrum  model  which  is  in  close  agreement  with 
empirically  derived  spectra  for  imagery. 

A.  SUBJECTIVE  TESTS 

The  following  four  types  of  images  were  evaluated:  (a) 
high  resolution  photograph,  (b)  low  resolution  photograph, 

(c)  infrared  image,  (d)  Synthetic  Aperture  Radar  image.  Each 
image  was  then  subjected  to  varying  amounts  of  spot  noise, 
line  noise,  line  jitter,  quantization  levels,  as  well  as 
reduced  sampling  rates. 

The  results  obtained  suggest  that 

• The  source  image  characteristics  form  strong 
contributing  factors  to  the  image  Quality  of 
a digital  image  processing  system. 

• Subjective  image  quality  is  a more  dominant 
image  quality  factor  than  objective  image 
quality. 

• The  system  parameters  of  the  state-of-art 
scann'-'-s  and  reproducers  can  serve  as  dominant 
image  quality  factors. 
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B. 

Several  simple  non-mechan i zed  tests  were  developed, 
including  the  use  of  Moire'  patterns,  to  measure  the  following 
objective  parameters  of  a digital  imagery  system: 

• Determination  of  Spatial  Sampling  Reproduction 
Ra  te 

• Modulation  Transfer  Function  Measurement 

• Scan  Nonhomogeneity 

• Measurement  of  reproducer  sampling  structure 

A simple  method  to  measure  the  imagery  system  gray-level  mapping 
and  the  number  of  levels  of  quantization  were  also  developed. 

All  of  these  objective  tests  are  easy  to  implement  and  require 
no  special  measurement  equipment. 

C.  CONCLUSIONS 

The  evaluation  of  digital  image  quality  has  lagged  behind 
developments  in  the  application  of  digital  techniques  to  imaging 
systems.  Surprisingly  little  is  found  in  the  open  literature 
on  digital  imagery  quality  even  the  important  considerations  of 
MTF  roll  off  versus  aliasing  errors.  This  preliminary  investi- 
gation including  both  subjective  and  objective  techniques  forms 
a basis  for  the  definition  of  a digital  imagery  test  set. 

The  recommended  test  set  criteria  should  include  measure- 
ments in  four  categories  of  system  performance:  non-linear 
amplitude  transfer  characteristic,  non-uniform  MTF  response 
(for  frequencies  up  to  one-half  the  Nyquist  frequency),  non- 
linear geometric  response,  and  noise  effects  in  general. 

These  criteria  have  their  roots  in  linear  system  theory  so 
they  should  be  good  descriptors  of  digital  image  quality 
especially  for  the  important  class  of  applications  for  which 
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the  desired  output  is  a close  replica  or  facsimile  of  the 
input  source. 

Digital  imagery,  based  on  its  periodic  or  array  struc- 
ture, has  interesting  properties  which  provide  simple  tests 
for  many,  if  not  all,  of  the  categories  of  system  performance. 

In  particular.  Moire'  patterns  can  provide  powerful  measure- 
ment techniques  which  exploit  the  array  structure  of  the 
digital  picture  elements. 

This  initial  study  has  provided  the  framework  for  a dig- 
ital image  test  set  with  continuing  efforts  recommended  for 
detailed  specifications  and  construction  of  a test  set.  One 
development  which  was  started  in  this  effort  and  worth  exploring 
further  is  the  relationship  between  MTF  response  and  aliased  dis- 
tortion. In  addition,  this  study  has  generated  computer 
programs  for  the  subjective  evaluation  of  digital  degradations 
which  would  prove  useful  for  simulating  the  performance  of  the 
sensor  portion  of  the  total  digital  imaging  system 
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DIGITAL  IMAGERY  TEST  SET 
PROGRAM  DOCUMENTATION 

Software  for  evaluation  of  the  Digital  Imagery  Test  Set  was 
developed  in  the  Certified  Prooessing  Laboratory  (CPL)  at  Harris 
ESD . 

The  hardware  configuration  of  the  CPL  is  shown  in  Figure 
A- 1 . Main  hardware  eguipments  in  this  configuration  are: 

• Datacraft  6024/5  CPU 

1.0  microsecond  cycle  time 

8-bit  and  24-bit  I/O  channels,  and  ABC  channels 

• Core  memory  ( 6 4 K , 24-bit  words) 

• Cartridge  Disc  System,  5.4  byte 

• Digital  Tape  Units 

1 each  (seven-track,  556/800  bits/inch) 

1 each  (nine-track,  800  bits/inch) 

• Punched  card  reader  (300  cards/minute) 

t Line  Printer  (200  lines/minute) 

• Anal og-to-Di gi tal  and  Di g i ta 1 -to- Ana  1 og 

Converters  (Harris  ESD  Design) 

• Teletypewriter,  ASR-33  Console  Teletype 

• Tektronix  Model  4014/15  Graphics  Terminal 

• Tektronix  Hard  Copy  Unit,  Model  4631 

• Harris  ESD  Imagery  Terminal 

• C0NRAC  CRT  Display 

An  overall  view  of  the  DITS  software  processing  system  is 
shown  in  Figure  A-2.  The  Digital  Processing  Subsystem  and  the 
Digital  Display  Subsystem  are  shown  in  Figure  A-3  and  A-4  re- 
specti vely . 

A digital  tape  with  6 image  files  was  received  from  Pattern 
Analysis  and  Recognition  Corp.  Contents  of  the  tape  are  listed 
in  Table  A-l.  The  data  was  digitized  at  8 bits  gray  level  per 
pixel  . 
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Detailed  descriptions,  documentation  and  program  listings 
of  the  software  used  to  evaluate  these  images  are  contained  in 
this  Appendix.  Section  A contains  program  descriptions.  Section 
B provides  program  listings,  and  Section  C contains  the  program 
flowcharts. 
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The  Digital  Display  Subsystem 


TABLE  A- 1 


I 

i 


1.  Synthetic  Aoerature  Radar  (UPD-4) 

over  Griffiss  AFB , NY  -(512x512)  pixels 

2.  B/W  vertical  photography  (RC-8  camera) 
over  Griffiss  AFB,  NY  -(512x512)  Dixels 

- approximately  (8-3/4  ground  ft.)1  per  pixel 

3.  B/W  vertical  photography  (RC-8  camera) 
over  Griffiss  AFB,  NY  -(1023x1024)  pixels 

- approximately  (1-1/4  ground  ft.)2  per  pixel 

4.  B/W  vertical  photography  (RC-8  camera) 

over  Floyd,  NY  resolution  target  -(1024x1024)  pixels 

- approximately  ( 1-1/4  ground  ft.)2  per  pixel 

F.  B/W  vertical  photography  (KC-1B  camera) 

over  Wright  Patterson  AFB,  Ohio  -(1024x1024)  pixels 

6.  Infrared  image  (sensor  not  identified  at  this  time) 
over  parked  aircraft  -(1024x1024)  pixels. 
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SECTION  A 


PROGRAM  DESCRIPTIONS 


PROGRAM  NAME: 

PARCON 

PURPOSE  : 

Convert  8-bit  packed  data  (32-bit  word) 
to  Harris  6024/5  8-bit  data. 

DESCRIPTION  AND 
METHODS: 

The  diqitized  images  received  from  PAR 
were  8-bit  data  packed  4 pixels/word. 

This  data  was  converted  for  use  on  the 
Harris  6024/5  computer. 

An  illustration  of  the  difference  in  data 
formats,  and  the  converted  format  used 
for  actual  processing  is  shown  in 

Figure  A-5. 

INPUT  VARIABLES: 

I T I 

Input  tape  device  number 

ITO 

Output  tape  device  number 

NWR 

Number  of  Words/Record 

NRF 

Number  of  records/file 

ITF1  - 

Tape  file  number  on 

Input  Tape 
(99=End  Program) 

ITF2  - 

Tape  file  number  on  output 
tape 

PROGRAMMING  LANGUAGE: 

FORTRAN 

PROGRAMMER  : 

G.R.  CUTHBERT 

HARRIS  ESD 
(305)  727-5284 

COMPUTER 

HARRIS 

6024/5 

OPERATING  SYSTEM: 

SERIES  6000  DISC  OPERATING  SYSTEM 

61516-01  EXTENDED  NON-SAU  FORTRAN 

COMPILER  REVISION  LEVEL  24-083076 
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Figure  A-5.  Data  Formats 


PROGRAM: 


DITS  MAIN 


PURPOSE: 


DESCRIPTION 

METHOD: 


This  is  a driver  program  used  to  execute 
one  or  a series  of  image  processing  sub- 
routines. The  subroutine  processed  is 
dependent  upon  the  value  of  the  variable 
ISUB  as  1 i sted  below: 


ISUB 

SUBROUTINE  CALLED 

1 

TODWIN  - Write  a desired  window 
to  tape  or  disk 

2 

SPAXFM  - Spatial  frequency  trans 
form 

3 

GRALEV  - Gray  Level  Histogram  to 
TEKTRONIX 

4 

DUMMY  - Not  used  at  this  time 

5 

DUMPSC  - Dump  a picture  to  the 
scanner 

6 

DUMPTV  - Dump  a picture  to 

CONRAC 

7 

GNCOPY  - Copy  a picture  from  one 
device  to  another 

8 

GRAMAP  - Gray  Level  Mapping 

9 

DUMMY  - Not  used  at  present  time 

10 

JITTER  - Introduce  jitter  into 
an  image 

11 

GNOISE  - Degrade  with  random 
Gaussian  noise 

12 

DIGDEG  - Digitization  Degraded 

13 

HGHALS  - Sampling  degraded,  high 
aliasing 

14 

REDALS  - Sampling  degraded,,  re- 
duced aliasing 

99 

END  PROGRAM 

AND  The  program  is  modular  in  construction, 

and  future  routines  can  be  incorporated 
with  ease. 

Format  for  input  to  all  programs  from 
cards  is  8110. 


INPUT  VARIABLES:  ISUB  - Subroutine  to  be  executed 

PROGRAMMING  LANGUAGE:  FORTRAN 


PROGRAMMERS:  DR.  R.H.  COFER,  G.R.  CUTHBERT 

HARRIS  "ESD 
(305)  727-5284 


COMPUTER:  Harris  6024/5 

OPERATING  SYSTEM:  Series  6000  DISC  OPERATING  SYSTEM 

61516-01  EXTENDED  NON-SAU  FORTRAN 
COMPILER,  REVISION  LEVEL  24.083076 


TODWIN 


PROGRAM  NAME: 
PURPOSE: 

DESCRIPTION  AND 
METHODS : 


INPUT  VARIABLES: 


PROGRAMMING  LANGUAGE: 
PROGRAMMER  : 

COMPUTER: 

OPERATING  SYSTEM 


Write  a desired  window  from  an  image  to 
a tape  or  disk  file 


Any  subset  of  an  image  can  be  written  to 
a new  output  file  by  specifying  the  line 
and  pixel  parameters  of  the  desired  win- 
dow 

M,N  - Main  picture  size 
(M  = no  of  lines) 

(N  = no  of  pixels/line) 

MM  - Location  of  first  line  in  window 

NN  - Location  of  first  spot  In  window 

II  - No.  of  lines  to  be  transferred 

JJ  - No,  of  spots/line  to  be  transferred 

I D I - Inout  file  number 

ITFI  File  no  of  input  tape/disc 

ID2-  Output  file  no 

ITF2  File  no  of  output  tape/disc 

FORTRAN 

G.R.  CUTHBERT 
HARRIS  ESD 
(305)  727-5284 

Harris  6024/5 

Series  6000  DISC  OPERATING  SYSTEM 
61516-01  EXTENDED  NON-SAU  FORTRAN 
COMPILER,  REVISION  LEVEL 
24.083076 
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PROGRAM  NAME: 

GRALEV 

PURPOSE: 

Plot  a histogram  of  image  grey  values  on 
TFKTRON IX 

DESCRIPTION  AND 
METHODS: 

Data  is  sequentially  read  from  the  input 
tape,  the  occurrences  of  grey  values  are 
counted,  and  the  mean  is  computed. 

These  results  are  then  displayed  on  the 
TEKTRONIX  CRT. 

INPUT  VARIABLES: 

M,N  Pictur,  size 

(M  = No  of  1 i nes ) 

(N  = No  of  pixels/line) 

PROGRAMMING  LANGUAGE: 

FORTRAN 

PROGRAMMER: 

G.R.  CUTHBFRT 

HARRIS  ESD 
(305)  727-5284 

COMPUTER: 

HARRIS  6024/5 

OPERATING  SYSTEM: 

Series  6000  DISC  OPERATING  SYSTEM 
61516-01  EXTFNDED  NON-SAU  FORTRAN 
COMPILER,  REVISION  LEVEL 

24.083076 
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PROGRAM  NAME: 

PURPOSE : 

DESCRIPTION  AND 
METHODS: 

INPUT  VARIABLES: 


PROGRAMMING  LANGUAGE: 
PROGRAMMER  : 

COMPUTER  : 

OPERATING  SYSTEM: 


SPAXFM 

Output  selected  spatial  transform  values 

A spatial  transform  is  executed  on  select- 
ed pixels  of  an  image  using  mathematics 
equivalent  to  a Fourier  transform.  Re~ 
suits  are  plotted  on  a TEKTRONIX  CRT. 

IU  - Row  number  of  pixel 

IV  - Column  number  of  pixel 

AMINR  - Minimum  value  for  X-axis  plot 

AMAXR  - Maximum  value  for  X-axis  plot 

AMINA  - Minimum  value  for  Y-axis  plot 

AMAXA  - Maximum  value  for  Y-axis  plot 

FORTRAN 

DR.  R.H.  COFER,  G.R.  CUTHBERT 
HARRIS  ESD 
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PROGRAM  NAME: 

PURPOSE  : 

DESCRIPTION  AND 
METHODS: 

INPUT  VARIABLES: 
PROGRAMMING  LANGUAGE: 
PROGRAMMER: 

COMPUTER: 

OPERATING  SYSTEM: 


DUMMY 

This  is  a dummy  subroutine  to  allow  for 
further  expansion  of  the  main  program. 

Establish  a temporary  subroutine  by  a 
call  statement  and  an  executable  state- 
ment 

None 

FORTRAN 

DR.  R.H.  COFER,  G.R.  CUTHBERT 
HARRIS  ESD 
(305)  727-5284 

HARRIS  6024/5 

Series  6000  DISC  OPERATING  SYSTEM 
61516-01  EXTENDED  NON-SAU  FORTRAN 
COMPILER,  REVISION  LEVEL 
24.083076 
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PROGRAM  NAME: 
PURPOSE  : 

DESCRIPTION  AND 
METHODS  : 

INPUT  VARIABLES: 


PROGRAMMING  LANGUAGE: 
PROGRAMMER  : 

COMPUTER: 

OPERATING  SYSTEM: 


DUMPSC 

Output  an  unpacked  6-bit  picture  to  the 
scanner 


8 bit  data  is  converted  to 
6 bit  data  and  read  out  to  the  scanner 
for  a hardcoDy 

ID1  - Input  device  no. 

ITF1  - File  no.  on  input  taoe 
II  - No.  of  lines  to  be  transferred 

JJ  - No.  of  pixels/line  to  be  transfer- 

red 

NFDPT-  No.  of  first  data  point  in  a line 
to  be  transferred 

NBITS-  No.  of  bits/pixel  on  input  record 
(6  bits  or  8 bits) 

FORTRAN 

DR.  R.H.  COFER,  G.R.  CUTHRERT 
HARRIS  ESD 
(305)  727-5284 

HARRIS  6024/5 

Series  6000  DISC  OPERATING  SYSTEM 
61516-01  EXTENDED  NON-SAU  FORTRAN 
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t 

PROGRAM  NAME:  DUMPTV 

PURPOSE:  Dump  a picture  to  the  CONRAC  CRT 

DESCRIPTION  AND  This  subroutine  has  not  yet  been  imple.4 

METHODS:  mented 


PROGRAM  NAME: 

PURPOSE: 

DESCRIPTION  AND 
METHODS: 

INPUT  VARIABLES: 


PROGRAMMING  LANGUAGE: 
PROGRAMMER: 

COMPUTER  : 

OPERATING  SYSTEM: 


GNCOPY 

Dump  a picture  from  one  device  to  another 

A utility  subroutine  which  gives  the 
capability  of  transfer  from  disc/tape, 
tape/disc  or  tape/tape.  Uses  BUFFER  IN, 
BUFFER  OUT  techniques 

ID1  - Input  device  no. 

ITF1  - File  no.  on  input  device 
ID2  - Output  device  no. 

ITF2  - File  no.  on  output  device 
II  - No.  of  records  to  be  transferred 

JJ  - No.  of  pixels/line  to  be  trans- 

ferred 

FORTRAN 

DR.  R.H.  COFER,  G.R.  CUTHBERT 
HARRIS  ESD 
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Series  6000  DISC  OPERATING  SYSTEM 
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PROGRAM  NAME: 

PURPOSE  : 

DESCRIPTION  AND 
METHODS: 

INPUT  VARIABLES 


PROGRAMMING  LANGUAGE: 
PROGRAMMER: 

COMPUTER: 

OPERATING  SYSTEM: 


GRAMAP 

Change  grey  level  values  as  desired 

Grey  values  are  changed  according  to 
breakpoints  and  the  new  desired  grey 
values  input  to  the  subroutine  on  data 
cards 

ID1  - Device  no.  for  input  pi c tore 
ITF1  - File  no.  on  input  device 
ID2  - Device  no.  for  outout  picture 
ITF2  - File  no/  on  output  device 
IMAX  - Total  number  of  grey  levels 
II  - No/  of  lines  in  picture 
JJ  - No/  of  pixels/line 
11,12  - Grey  level  mapping  values 
13  - Flag  to  stop  reading  11,12 

FORTRAN 

DR.  R.H.  COFER,  G.R.  CUTHBERT 
HARRIS  ESD 
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PROGRAM  NAME: 
PURPOSE: 

DESCRIPTION  AND 
METHODS  : 

INPUT  VARIABLES: 


PROGRAMMING  LANGUAGE: 
PROGRAMMER: 

COMPUTER: 

OPERATING  SYSTEM: 


JITTER 

Introduce  random  jitter  to  degrade  a pic- 
ture 

Jitter  is  introduced  at  random  points  in 
an  image  through  use  of  random  number 
generator 

ID1  - Device  no.  of  the  inout  picture 
ITF1  - File  no.  of  input  if  tape 
ID2  - Device  number  of  output  picture 

ITF2  - File  no.  of  output  if  tape 
II  - No.  of  lines  to  be  transferred 

JJ  - No.  of  spots/line  to  be  trans- 

ferred 

ISD  - Standard  deviation 

IMN  - Mean 

ISEED  - Seed  for  GRN  generator 

FORTRAN 

G. R.CUTHBERT 
HARRIS  ESD 
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PROGRAM  NAME: 
PURPOSE: 

DESCRIPTION  AND 
METHODS: 

INPUT  VARIABLES: 


PROGRAMMING  LANGUAGE: 
PROGRAMMER: 

COMPUTER: 

OPERATING  SYSTEM: 


GNOISE 

Introduce  Gaussian  random  noise  into  a 

picture 

Gaussian  noise  is  introduced  into  a pic 
ture  through  use  a random  number  genera 
tor.  The  noise  can  be  added  to  each  pi 
el  or  by  line. 

ID1  - Device  no.  of  input  picture 
ITF1  - File  no.  of  input  if  tape 
ID2  - Device  no.  of  output  picture 
ITF2  - File  no.  of  output  if  tape 
II  - No.  of  lines  to  be  transferred 

JJ  - No.  of  spot/line  to  be  trans- 

ferred 

IMAX  - Maximum  picture  value 
ISD  - Standard  Deviation 
IMN  - Specified  mean 
ISEED  - Seed  for  grn  generator 
ISL  = 0 - Add  a different  RV  to  each 
spot  in  picture 

= 1 - Add  a Different  RV  to  each 
line  in  picture 

FORTRAN 

G.R.  CUTHBERT 
HARRIS  ESD 
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1 


PROGRAM  NAME: 

PURPOSE: 

DESCRIPTION  AND 
METHOD: 

INPUT  VARIABLES: 


PROGRAMMING  LANGUAGE: 
PROGRAMMER: 

COMPUTER: 

OPERATING  SYSTEM: 

\ 


DIGDEG 

Picture  is  degraded  digitally 

The  number  of  bits/Dixel  is  changed  by 
shifting  bits  out  according  to  values  on 
input  cards. 

ID1  - Device  no.  of  input  picture 
ITF1  - File  no.  of  input  if  tape 
ID2  - Device  no.  of  output  picture 
ITF2  - File  no.  of  output  if  tape 
II  - No.  of  lines  In  output  picture 

JJ  - No.  of  spots/line  to  be  trans*- 

f erred 

NBIT  - No.  of  bits  to  be  shifted  (right) 

FORTRAN 

G.R.  CUTHBERT 
HARRIS  ESD 
(305)  727-5284 

HARRIS  6024/5 

Series  6000  DISC  OPERATING  SYSTEM 
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COMPILER,  REVISION  LEVEL 
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1- 

■ 

PROGRAM  NAME: 
PURPOSE: 

DESCRIPTION  AND 
METHODS: 

INPUT  VARIABLES: 


PROGRAMMING  LANGUAGE: 
PROGRAMMER  : 

COMPUTER: 

OPERATING  SYSTEM: 


REDALS 

Reduce  a picture  in  size  by  low 
aliasing 


Reduce  the  size 
on  it  in  Dowers 


of  an  image 
of  two. 


by  operating 


ID  1 

- Device  no. 

ITF 1 

- File  no . o 

ID2 

- Device  no. 

ITF2 

- File  no.  o 

II 

- No . of  1 i n 

JJ 

- No.  of  spo 
sferred 

ISIZE 

- Power  of  2 

IMAX 

- Maximum  pi 

of 

i npu 

t p 

i c 

ture 

i 

nput 

if 

ta 

pe 

of 

outo 

ut 

Pi 

cture 

0 

utput 

if 

tape 

s 

to  be 

tr 

an 

sferred 

s/ 

1 i ne 

to 

be 

tran- 

re 

ducti 

ons 

1 

n size 

tu 

re  va 

1 ue 
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PROGRAM  NAME: 

PURPOSE  : 

DESCRIPTION  AND 
METHODS  : 

INPUT  VARIABLES: 


PROGRAMMING  LANGUAGE: 
PROGRAMMER  : 

COMPUTER: 

OPERATING  SYSTEM: 


HGHALS 

Reduce  imaqe  usinq  hiqh  aliasinq 

Reduce  the  size  of  an  image  by  operating 
in  powers  of  two 

ID1  - Device  no.  of  inout  Dicture 
IFT1  - File  no.  of  input  if  tape 
ID2  - Device  no.  of  output  picture 
ITF2  - File  no.  of  output  if  tape 
II  - No.  of  lines  to  be  transferred 

JJ  - No.  of  SDOts/llne  to  be  trans- 

ferred 

ISIZE  - Power  of  2 reductions  in  size 
IMAX  - Maximum  Dicture  value 

FORTRAN 

DR.R.H.  COFER,  G.R.  CUTHBERT 
HARRIS  ESD 
(305)  727-5284 

HARRIS  6024/5 

Series  6000  DISC  OPERATING  SYSTEM 
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SECTION  B 


PROGRAM  LISTINGS 


NAME.  PAkCUN 

')  I PENSION  I HUP  1(344),  IHUFef(  103?) 
uATA  ICk,ILP, ITTY/7,6, 1/ 

C 

C CONVEkT  8-HIT  PACKED  DATA  (ip-pi!  «UkD)  Tu 

C 10  HAkWIS  60?a/b  H-hIT  DATA  ( (l  PIXEL/kUKO 

C wOPUS  PlOHT  JUSTIFIED 

c 

(.  PkOUP  AMtK  : OFkl  Cul  HbEPT 

L PkUoPfAF,  DATE: 

C 

C V Ak I AbLES 

C I T I 

C ITU 

C Nrt  k 

C NkF 

C I TP  1 

C 1TF? 

C 

kEAlH  ICH,8000) I T I , ITO 
'10  CUNT  iNUt 

KEAU(  lLR,S0l)0)NkF,N/«k,  ITF1,I1F? 

c advance  input  tape  FiLt  if  nelessapiy 

IF ( I TF 1 .E 0.49) GO  TO  VV9 
IFUlFl.Eu.DbO  TU  PO 
f)U  15  L = ?,  1 TF  1 

* 1 5’ Call  ADF ( 1 1 I ) 

c advance  output  tape  file  if  nelessaky 

Pu  IF  ( 1 lFP.tU.  1 )U0  TO  31) 

DU  <>5  L=?,1TF? 

<?5  CALL  ADF (ITO) 

30  CONTINUE 

C SET  DEFAULT  PIC1UPE  SIZt  TU  lo?4  ( 3^4  DC  Y.'UkDS  PACKED) 

C anU  CnAimOE  IF  NECESSAnY 

1EN0=34« 

*kEC=342 

IF  ( iV •» k . N E . 1 0 ? 4 ) 1 E N D = 1 7 i 
IF  (NKk.NE  . 10?4)nREC=  1 7 1 
i)U  500  1 1 = 1, NkF 
( CLEAk  OUT  8UTM  eUFFEHS 

f-0  35  L = 1 , 3 u 4 
3b  I bUF 1 (L ) =0 


used: 

- input  tape  uEvict 

- OUTPUT  1 APE  DEVICE 

- NU  OF  wUKOS/ktCUki) 

- NUMPEk  OF  WECOkuS/FlLE 

- TAPE  FILE  NO  ON  INPUT  TAPE  (94  = PkuUnA*  t 

- Tape  file  NU  ON  OUTPUT  1 APt 


Ni; ) 


BEST  AVAILABLE  COPY 
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IMl  i*  I . S J , 1 (i  32 
3»  1 SUF  2 IL  ) =0 

C 

C.  t,h  I 1 rj  A t*llF  l-t-K 

rit)F  F fcH  iMlTI,l  HUF  1 , -t , iwt  L , 1 S ) A 1 ) 

JO  CALI.  ST  A I US  ( I T 1 ) 

i ,»•  III  (<JO#txi*SO),]STAf 
C 

C rKuCfcc<S  1 's  I^C^t^tNTS  OF  J .‘.OPI'S 

Sil  M<  = -2 

uu  ioo  i=i,itNU»a 

00  20  0 J=1*J 

M\  = K K ♦ 3 
JJ=1 +J- 1 

J A = 1 HUF  1 (JJ)  .AM>.  ' 377 
insI<?jF  i (JJ)  .anu.  ’ 177400 

1 C = 1 PUF  1 (JJ)  .AND.  *7  7600000 
ImUF2(kk)=2SS-(  I C. ShI FT.- 16) 

1 t'UF  2 ( K K ♦ 1 ) =2S5-(  IH.sHlF  T ,-P) 
lrtUF2lK**£)=2SS-IA 
20  0 CuuT  1 1, lit 
300  CONTINUE 

c: 

r MljFFtP  r fiU  CUNTaIWS  OUlPUf  L»  A l a - Put  ll  U'M  lAPt 

^UF  FtP  ull  I ( i TO  , IHI'F  2 / *1 1 urtH  * ISlAf  } 

230  CALL  STATUSUTU) 

UO  TO  (23O.SOO,SO0) , 1SI AT 
SOO  CUNriNUt 

c 

C F J L t FlNlSHtO  akITE  EOF 

t »U  F 1 Lt  l TO 

C •**  IF  YOU  aunT  T U PNOLtSS  1 Ht  NtXl  F1LF.  ftllhOUT  Kt  «•  1 U>.)  i M- 
C *«*  1st  TAPES  AM)  A[)  v A YC  1 r.(i  THEN  AOAlN,  NaKE  11M  K IiF2  = l 
l,i)  I')  10 
>iVO  PtrtINu  1 T 0 
LmLL  EXIT 
MIl'O  FukM  AT  C * I 1 0 ) 

C NO 


* 
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BEST  AVAILABLE  COPY 


CCMM0N/IDATA/IBUF1 (4096) , IBUF2 (4096) 
DATA  ICR, ILP , ITTY/7 ,6,1/ 


C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

r 


PM0GP4F'  PUPPO.st  : THIS  16  4 IPIVFk  Pk  U,KM«  IF'Fu  TO 

c * C r I J T t < I N L ( • K A S t * I F.  b I • K 1 1 u 3 r 
PPOCt SSlub  SrriHiiiT)  MtS.  THF 
SI'HtfULiflNt  H-'Ol.F.SStO  la  liF  hF  T 
UPijN  T HF  VALUF  HF  I Tt  vAKjAbLt  I bUH 
AS  L IblFu  HELUtr  : 

ISob  SO 9POUT  I «F  C A L L t U 


1 TOOv- 1 0”nT<  1 1 1 4 OF  S 1 k1  fc  u >vll"UO»', 

I A^h  Til  O J Six 

?.  6 P 4 X F r 1 - S ► 4 ( I m L F k P (j  (J  t ft  C Y f H A N S F ij*«  "X 

3 bKALfV  - Okay  LtvFL  rl  1 6 T U 

1 U T t K | K I ; N 1 Y 

0 U6i«»  i — iv  U T 1 1 S f u m f P K F SF.  iv  T I 1 ifc. 

6 01P-.P6L  - HUSH  PlCTLWL  10 

6 C M a . i t F l 6-H IT) 

h b Uv  P I w - DbKP  P T C T 0 t TO  Tv 

7 bftCOPY  - copy  t PICTOaf.  From  OMt 

uf.  i/ 1 C t TO  m «0 T f-F P 

b GWA^AP  - Ui-«Y  LfcYtl  CAPPING 

9 oOkiX'Y  - Mil  USED  AT  PPcSEftT  (Ik 

10  J]|TF;k-JJT!£« 

11  bNUlSfc  - lit  UK  A u£  o I T H KAivOUft 

b A U S S I A fy  ft  U 1 S F 

12  UJbOtU  - niliiTi  7ATi  ;JN  i)Eti«AOtU 

13  HGhALS  - S A 4PL 1 A1G  UfUPaUci)  - 
N 1 b 0 A L 1 A S 1 0 b 

10  kFoALS  - SAPPLIUG  OFbPAOFO  - 

kFOULFu  A L 1 A 6 1 i\  b 
99  - TM)  P K 0 1-  K A M 


PkORKaM  PATt:  2/19/77 

PkORK  A kS  : b t k I ClITHFFKT/  Uk  k.H.  I’.OFF.P 


OOTf.:  ALL  SUBKUOflbEa  v.Il_L  Kt  KF.SP(  ns  IPLfc  T Oh  K t Al)  1 ilb 

u * I A KFiiOIkFP,  Kfc*lftOlNb  uF VICES,  « i'f>  O'k  I f 1 1 j 0 E uF  ' S . 

F OK'F  A T FOk  INPUT  fhtl'r-i  f 4 K u 6 1 6 Mil;) 
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( ' ’ r I r K '•*  I 'V.  fe  liMlth  SUt'KIiU  • i Nh.  lu  bXtCOlfc 

l ll  I UU  f I Ull*- 

AO  ( J C * , 6000  ) 1 SUn 
lr  ( ISUb.F U. V4)U0  lu 

UU  I U ( 1 (it  t 2 0 0. 30u  , 400  , *500 , 600  , /U0  * 80  0 , «O0  , 1 OuO  , 
i 1 1 oil,  1P00,  1 *0  0,  1 400)  , 1SU8 

C 

C rXtlOlF  SUFlPUOllNt  I 1 1 0 1~  l >M  - NPIIF.  AliMOUA  St  Ul/fc.  N t 

1 in.'  LuMI  1 -lit. 

C 


L 

SUbrtllUT  I l\ t 

vAuIASLtS  wfAL'  i nj  IUUvIu  - 

8110  FUKMA  T 

L 

» F W 

- 

KA  I N HJCll'Kt  S I / fc  ( P = N 0 IF 

L 1 UtS, 

r 

N = M)  OF  SPUTS/LINF) 

c 

nm 

- 

LuLATIOiv  UF  F 1 PST  LIMt  In 

W 1 uDu»» 

r 

•\N 

- 

LOLA i ION  »JF  »-lKSI  OcSiWtU 

SPOT  IN  N'M 

r 

1 1 

- 

NO  OF'  LlFiF.S  lU  8 h f w Ai'lSF  tKHt  0 

c 

JJ 

- 

MU  UF'  SPOT /LI  Alt  TO  bP  T n « 

USF  F.KPpD 

c 

101 

- 

I NPlJ  T FILL  uUMbt  P 

c 

ITF  1 

- 

FILt  NO  LU  INplJI  If  1 APt 

L 

ID? 

• 

OUTPUT  FJLt  No 

i 

ITF<> 

F l Lt  lull  OF  UU  1 HlJT  IF  1 APt 

c. 

C w t AD ( 1 C R , 80  0 0 ) M , N , MP , oN , 1 I , J J 

C Rt»D(  ICP,8P00)  I U t , l IF  1 , Iu<?,  1 1F<? 

C 

call  iou^im 
ou  to  io 
L 

?u0  CUoTlNlit 
L 

(.  FXtCUIt  SUhwUU  f lut  SPaxFi*'  - CUPHUTk  1 Hfc  .SPA  I I A L T h A .m  S r iJkm 

L 


f. 

IF  1 

- INF-ui 

P ILt  UUMBtH 

l. 

10<> 

- UISa 

Flit  UF  NlODllrttD  A 

ut  A 

c 

IPS 

- S T Aw  r 

Pun.  OF  fht  a IN  DO  a 

c 

N 

- SIZt 

OF  sr.  InoUn  (MUST  OF 

SuUARfc  ) 

CALL  1 01>A  I n 

wFA|)  j i\  lilt  NtCtSSMPY  VuwIAbLfcS 


w b A |)  ( 1C^»HOo0)  I L' «?  # I-j 
LmLL  SPAXFM  IU«r,N) 
i.n  to  lu 
c 
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f copy 


C 


c 

c 

c 


c 

c 


c 

f 

l. 

r. 

c 

c 

c 

c 

r. 

c 

c 

L 


l 

c 


tXfCUTE  b k L E - GkAy  LtvEL  hi  I f-  T 0 10  IkM  KuhUs 

J0O  f.GNTINuF 

EXECUTE  bPALtV-  PLUl  n l S 1 huh  <i.v  uF  UE  A Y VALUES 
(.ALL  (>K  a l e V 
bO  TO  10 

E.xeCU'lfc  SIJhkUU  T i Mk  - SF'ECTkAl  Fki  J xFE<p 

0 00  CALL  DUW'Y 
1,0  Tu  10 
S00  CONTlNUt 

kXtCUTE.  SUHfOUT  i NL  uUi-'PsC  - Uu^P  PiLTbkt  lu  Sc  A wA  c * ( l f ) 

SUhPuUTlUt  VAKlAHLtS  r,  E A 0 1m  OU'VPSC 
I D 1 - 0 1 Sf  ElLk  NummEk 

1TF  1 -FILE  NO  I'F  INPUT  ]F  TAPE 
11  - NO  UF  L i *0 F S III  HE  I h A iM  SF  t P K E 

JJ  - Nil  OF  SP0T/L1')E  T u f<E  1 KA'MSFknPtO 

nFuPT  - 1SI  l)A  I a PlJlAif  m HE  TkAMSFc.  YkEo 
ImtITS  - i\  U M S k k OF  H I Tb/H  1 /r  L ( F-Hll  uA|A  ;jk  o-r.ll  i ■ A f A ) 

WtAU(IChi,rtOOO)  101,  J IF  1 , I I , J J , i«FDPT  , >ME*  1 f S 

Coll  lhmpsl 
bO  iu  lo 

FxtCUTk  SUHkuUT  liYE  UUwPTV  - I'UHF  PICTURE  TO  Tv 


c 

r. 

c 

c 

c 

r. 

<; 

c 

c 

c 


#>00  LOLL  DUMMY 
UiJ  TU  1 0 

7 on  Cum i i \u k 

EXECUTE  SUUKOUTlWk  bOLOPY  - Copy  A PICTURE  FKUi'i 

U')L  UkViCt  TO  Arid  THY  f 
SUPkUUl  INF  VAFldHLko  fEAU  I M GNCUE’Y 
I D 1 - Input  UE  VICk  iMUf-'hE  h 

I T E 1 - FILE  NO  OF  IM^UT  IE  I APt 

102  - UUTPUT  UtvlLt  NUN^kF 

I T E 2 FILE  <vU  UE  OUTPUT  IF  TAPE 

II  - NO  UF  LlWtS  IU  HE  T H AimSE  kkKt  0 

JJ  - MU  OF  SE'U  I /L  1 imE  IU  Ok  Tf  AimSE  EkKkO 


M1  AP!  £ 

| il  Y -S  kU 


190 


r 

* R W ( ILK,  t)  000)  lull  » I 7 F 1 , 1 D<? » II  F J , I 1 , J J 

r 

• ALL  UN c DRY 

fiu  ro  i a 

i: 

MUD  L l)  "!  T 1 N U t 

I.  EXECUTE  SUBROUTINE  GK  «M  AP  - GHAY  LtVtl  MAPPING 

L VA-.  TABLES  TO  riE  t A l » 1 N BY  (,hfii<iAP 


c 

1 01 

- DEVICE  N 

0 OP  INPUT  PIlTURE 

c 

I tf  1 

- E I L E ivO 

Ur  INPUT  IF  TAPt 

c 

I D<? 

- DEV  ICE  N 

0 UP  UUIPUl  PICTURE 

c 

I TFP 

RILE  NU  OE 

UUrPUt  IE  1 APE 

c 

I max 

- MAXIMUM 

PICTURE  VaLUF 

c 

1 1 

NO  OF  LINE 

S Tu  mE  TPANSFEKRtU 

c 

1 

JJ 

- NO  UF  SPUTS/LINE  TO  BE  TRANSFERRED 

L 

c 

Re  Au ( I ( K , 

RO00HD1,]  IF],  ID<?,1TF2,IMAX,II,JJ 

C 

C «L  L GRAMAP 
fill  TU  10 


C EXECUTE  SUdROUT  i ME  LINInT  - LINE  INTENSITY  DEuRADAT  ION 

900  CALL  DUMMY 
00  TU  ID 
C 

l Execute  subroutine  jitter  - jitter 

1 iMirt  Cun  f i nue 

C VARIABLES  to  me  REAdIN  IN  BY  JITTER 

c ioi  - device  no  of  the  input  piciure 

C ITEl  - FILE  NO  UF  INPUT  IF  TAPt 

l-  102  - DEVICE  number  UF  OUTPUT  PICIURE 

C IIF2  FILE  NO  UP  UUTPUI  IF  TAPE 

C 11  NO  OF  LINES  TO  BE  TRANSFERRED 

( JJ  - NO  UF  SPUTS/LInE  TO  BE  TRANSFERRED 

C ISO  - STANDARD  DEVIATION 

C 1MN  - M E A i* 

C iSEfcD  - SEED  FOR  URN  uEneratgr 

c 

C RE  A0( ICR, BOOO) ID  I , I7F 1 , £02, I TFP, I I , JJ 

C REATM  ILR,oOOO) ISD, IMN, I SEED 

c 


DEGRADE  rtlTh  K A NOUN  IjAubSlAN 


CALL  JITTtK 
UU  TO  111 
1100  CONTINUE 

C EXECUTE  SUBKUUT  1 NE  GNuISE  - 

C MU  I St 

L SUBkOUTIM-  VAkIAHLES  READ  IN  UMOlSE 

C 101  - DEVICE  NU  OF  INPUT  PlCTUKt 

C I TF  I - ULt  NU  OF  INPUT  IF  TAPt 

C ID?  - DEVICE  NO  OF  OUTPUT  PlCTUKt 

C I T F ? F j LE  NO  OF  UUTPUI  IF  TAPE 

C II  - NO  OF  LINES  TO  HE  TkAnSFEKRU) 

L JJ  NO  OF  SPOT /LINE  TO  tiE  T h AmSFERPtD 

C " I n A X - M a X 1 Mom  PICTURE  VALUE 

C.  ISO  - STANDARD  DEVIATION 

C ' ~ IrtN  - SPECIFIED  Mt An 

C 1 SEED  - StEU  FUK  GRN  GENERATOR 

L " 1SL  = 0 - AuD  A DIFFERENT  KV  TO  EACH  SPOT  IN  PlLTUPc 

C = 1 - AUD  A DIFFERENT  KV  IU  EACH  LINE  In  PICTURE 

C 

C wt  Al>(  ILK, *000)101 , I TF  l , ID?,  I TF?,  1 1 , JJ 

C RE ADI  ICk, *000 ) ImAx , I SO, I BN, IStfc  V , 1SL 

C 

CALL  GNU  I St 
GO  TU  10 
IP 00  CONTlNUt 

C EXECUTE  SUFiKOUT  lNt  OIUDEG  - UIGI  T IZAT10N  DEGRADED 

C 

f VAKIAHLtS  KtAU  IN  DIUJEG 

r.  I D 1 - DEVICE  NO  OF  Input  PlCTuRt 

C I TF 1 - FILE  NO  OF  INPUT  IF  TAPt 

C 10?  - DEV  ICt  NU  OF  OUTPUT  PICTUkE 

C I TF?  FILE  NO  OF  OUTPUT  IF  (APE 

C II  NO  OF  LINES  Tn  uE  TrAmSF EKRtO 

C JJ  - NO  uF  SPOTS/LINE  TU  HE  IkanSeERPF.D 

C N0IT  - NO  OF  HITS  TO  ShIF  (CANNOT  ttt  0) 

C 

C Kt'AI)(ICR,P00  0)  101,  ITF  1 , ID?,  I TF?,  II,  JJ,NHi  T 

L 

call  oiadfg 

NO  TO  10 
1 3oO  CONTINUE 

C EXECUTE  SUHrOoT INE  hGhALS  - SAMPLING  OF  GRADE U Hluh  ALAIS» 

ING 


f 


I 


r. 

i; 

c 

c 

c 

c 

c 

c 

c 

c 

c 


VANI  Abl.ES 
11)1 
I TK  1 
I 02 
1 TK  2 
I I 
JJ 

1 S 1 2 1 


HEAD  ]H  HOMALS 

- DEVICE  NO  UK  1 r-!PiJ  I MCTOPE 
-KILE  NU  OK  INPIM  IK  TAPt 

- DEVICE  Ml  UK  uuTPul  I C I UK  t 
KILE  NO  UK  UU I HU  T IK  TAPE 

NO  UF  LINKS  Tl)  b£  TwAMSKEPPED 

- NU  Or  SPOTS/LINE  lU  bfc  TPANSKEPkEU 

- PUwER  Ul  2 REDUCTIONS  I IM  SlZt 


PEA|)(  ICp,UOOO)  Iul , 1 TH  , l IV,  1 TK2,  1 1,  JJ,  KSI2F 


CALL  HOmALS 
ou  10  10 
IaoO  CUM  I NOE 

C EXECUTE  SUBROUTINE  PEDALS  - SAPPLlMi  OEuPaDEU 

C PEDUCEO  ALIAS 1 NU 

CALL  PEDALS 
UU  TU  10 
C 

C END  PROGRAM 

RRRR  CONTINUE 
CALL  EXIT 
8000  K OPM AT  C«I  1 0) 
t NO 
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' SUhROUTINt  I UOw  I N 
CUMMON/IDAT A/IMUFI (4046) , IRUF2U046) 

OAT  A ICR,  ILP,  lTTY/7,6,  1/ 

SUHHOUT I Nfc  PUHPuSt:  WRITE  A OESIkEO  WINDOW  SEUMENT  FkOm 

TAPE  TO  UlSK 

INITIALIZE  lHUFl  ANO  1HUF2 
DTTlO  K= 1,4046 
IBUFI (K J=0 
"R  TflUF2(K)=0 


C 


READ  IN  NECESSARY  VARIABLES 

M,N  - MAIN  FICTukE  SIZE(F=NU  UF  LINES, 

NsMJ  OF  SPOTS/LINE) 

MM  - LOCATION  OF  FIRST  LINE  IN  wlNOuw 

NN  - LOCATION  OF  FIRST  UESIPtO  SPOT  IN  MM 


IT  - NO  OF  LINtS  TO  HE  T«  Af>SF  E'KHEU 

JJ  - NO  OF  SPOC/LINE  TO  HE  TRANSFERRED 

~ “1D1  - Input  file  vuMbER 

ITFl  - FILL  NO  OF  INPUI  IF  TAPE 

TD2  - OUTPUT  FILE  NO 

ITF2  FILE  NO  OF  OUTPUT  IF  IAPE 
READ ( ICR, H00  0 )M, N, MM, wN, 1 1 , J J 
Rt  AD( ICR , bOOO ) 101, ITFl, 102, ITF2 
"REMIND  101 
REWIND  102 

SET  THE  VALUES  OF  IWORU  AND  IIInE 
R0NDsINN  + JJ)-1 


IlInE  =(MM+II)-1 


c 

C ■ ADVANCE  TO  CUKPECT  FILE  NO  IF  NECESSARY 
IFCITFl.LE .1 )U0  TO  3 
UU  1 w = 2 » I I F 1 
1 CALL  A0K101) 
i CONI INUt 

IF  ( 1TF2.LE  .1)00  Tu  S 
00  4 Ks2,ITF2 

4 CALL  AUP(I02) 

5 CONTINUE 

c advance:  to  first  line  in  window 


BESI  AVAILABLE 
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IF (MM.Eu. | ) GO  TO  70 
j AUV=*M- 1 
00  bS  K=1,IADV 

riOFFEK  1N( Iul , IbUF 1,B,N, I ST  AT) 

50  CALL  STATUS (101) 

60  TO  (SO,  65, bO)  , 1STAT 
Dll  akITE(1ITY,9010) 

ST  OF* 

6S  COnTI Hi UE 
/O  CONTINUE 

L MEAL)  IK  t»  I NOOtt  A WO  WRITE  TO  DISw 

00  100  I=HM,ILINE 
HUFFER  INUD1  , IbUF  1 ,P,N,  ISTAT) 
bO  LALL  STATUS (101) 

60  TU  (60,85,90)  , 1ST  AT 
*0  •»HITE(1TTY,9020) 

SI  OP 

bS  CONTINUE 

c 

C CHECK  TU  SEE  IF  THE  WHOLE  LINE  SHOULD  60  TO  THE_DISK 

IF  ( IwOkl).EO.N)Gu  TO  300 
C 

C TPANSFEh  THE  NECESSARY  nORUS  TO  IBUF 2 

KK*0 

00  200  KsNN, I wOKO 

KlMKKt  1 

71 0 1HUF?(KM  = IHUF1  (K) 


WRITE  THIS  BUFFER  TO  THE  OISK 
BUFFER  OUT ( 102, IbUF 2,rt, JJ, ISTAT) 

Hi  0 CALL  STATUS (102) 

60  TO  (210, 220, 220), ISTAT 
'22o  GO  TO  100 

C WRITE  THE  WHOLE  LINE  TO  ThE  OISK 

300  hiiFFER  UUT(ID2,IBUFl,fi,N, ISTAT) 

2«0  CALL  ST  A T US ( 1 02 ) 

60  T0(290, 100,100), ISTAT 
tOO  CONTINUE 

C 

c DATA  TRANSFER  SHOULD  BE  COMPLETE 
END  FILE.  102 
REWIND  10) 

REWIND  102 
RETURN 

6000  FORMAT (61 10) 

90  10  FDKF-AT  ( IX, ’EOF  ENCOUNTERED  IN  ADVANCING  TO  wInDOW’) 

9020  FORMAT ( IX, 'FUF  ENCOUNTERED  DURING  WINDOW  PROCESSING’) 
tNl> 


SonROuT I Nfc  SP  A XF  M ( I O? , N ) 

C UUTPUI  SAMPLE  SPATIAL  TRANSFORM  VALUES 

CCiMP.(WI0AT4/IbUEl  (40R6),  Ib'JF2(  40Ro) 

DIMENSION  IU(S12),IV<512),A(S12),P(bl2),M(Sl?) 

DATA  IYES/3HYES/ 

C f 

C 1U  RUrt  NUMBER 

t IV  - COLUMN  NUMBER 

C ' R - SORT (IU**2+IV**2)--(VALUES  F (JR  X-AXIS  PLU1) 

C A AMPLITUDE  VALUES 

c p phase  values 

C I DP  - DISK  FILE  OF  PICTURE 

C 

PI=i. 1415S26S 

t 

C PE AO  IN  X MIN  AND  MAX 

E READ  IN  Y MIN  AND  MAX  -IN  OB) 

RtA0(7»S0)  AMI,\»,  AMAXR,  AMINA,  AMAXA 

SO  F ORmaT ( uF 1 0 . 3 ) 

nRUNsO 

— r o' o- co  nt  Ivor 

NS  = 0 

VhUNsNR'JN+1 

C READ  IN  SPATIAL  SAMPLINU  FREQUENCY  LOCAIIONS 

1 READ ( 7 » 2 ) 1U0,IVD, IFLAG 

2 F UWMA  T ( 2 I b » 1 1 ) 

MSsNStl 

A(NS)=0. 

P ( NS ) = 0 . 

IU ( NS ) *1UD 
1 V l NS) =1 VD 

R(RS)=S'jRT  (FLOAT  ( J 00 ) **2*FL0AT  (IVD)**2) 

IF ( If L Ao.EQ.  1 ) GO  TO  3 
GO  TO  1 

C 


C START  THE  TRANSFORM  PROCESS 

"*"J  C0N$T=-t>.2831ttS3/FLUAT(N) 

C DO  THE  ROwS  FIRST  FUR  MINIMAL  DISK  ACCESS 

REwIno  in? 
do  a ix 1=1 ,h 
Ik* IX l-i 

BUFFER  I N ( 1 02  # I RUE  1 » B » N, ISTAT ) 


n 


2'J  C«LL  ST  A 7 u S ( I U 2 ) 

GO  TO  ( 20 , 25 » 25 ) » J S f A T 
25  continue 

C l»0  THE  COLUMNS  next  FU«  MINIMAL  DISK  ACCESS 

OU  a IY1=1,n 
I Y=I Y 1-1 

C GET  IHt  SHOTS  VALUE 

Y4l  = FLI’aT(  IiUlf  1 ( IY1  ) ) 

V«L=1  . 1**VAL 

L ADO  In  EFKFCT  A l fcACfi  SAMPLE  SPATIAL  FKEJUF.nCV 

O'J  a k = l , n S 

AKG=CONS  I AFLOAT  ( IX*1U(K)  tl  YMlV(K)  ) 

O TFMHOPAhLY  STOKE  HEAL  PAKT  IN  H,  IMAG  PANT  lN  A 

AIF)=A(N)-VAL*S1NIAKG) 

P(H)=P(K)+VAL*LUS( AhG) 

« CONTINUE 

DO  FINAL  FMOCF.SS  I Mi 
i)0  5 K=  1 , NS 
A(a)=A(K)/F'LUAT  IN  J **2 

p(k)=f-(k)/F'loa r (nj**2 
Ad  SOLU  =1U.*ALUG10(AIK)**2+H(K)**2) 

PUW=ATAN<r(A(iO,P(K)) 

A ( K )=AHSULU 
5 CONTINUE 
VKlTElbrfe) 

h F uk MAT ( ' SAMPLE*, U, V, AMPLI luOE,  PhASE') 

00  7 1=1, NS 
ANG=P(1)*57. 2957b 

/ WIIE(6,M  I,IU  (1),1V  ( 1 ) , a ( i J » ANG 
b FUKMAT (IX, 3110, 2 (3X,F15.8)) 

CALL  PLOT  IT  (NKlJN,NS,K  ( 1 ) , m ( 1 ) , 0. 0, 0.0,0. 0, 0.0, 1 0.0,  10.0) 
IF  (NMjiV.LT.  3)  GO  Tu  10  0 
CALL  HOCtiPV 
«E lUKN 

e nu 
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SUHROUTINt  GPALEV 

COMMON  / 1 1*  A I A/IhUP  I ( Al  0 S»  o ) , I t>UP  ? ( A 0*to  ) 

DIMENSION  lGPA(2Sfc) , Ax(2S*i) , At  (25o) 

DATA  JCP,1LP,IFk  /7,6,2l/ 

C. 

SUHKUUl  lNt  PUnPuStt  PLOT  A PlClUWfc  hISTUUKAN  UN  TEKTkUNIX 
VAP1 aulES: 

M,N  - PICIUkE  SIZE  (Msfili  UP  LINES,  N = NU  OF  SmuTs) 
iui  - Input  hu 

PL  AIX  ICh, 8000) 101 ,w,N 
o 000  PUP  MAT  (ft]  10) 

C 

c initialize  igka  iuffe* 

DU  10  l=l,2b6 
10  IGWA(1)=0 
ASU*!  =0.0 

ADlVsFLOAT («)*FLO»r(N) 

c 

(.  GET  In  UATA,  COUNI  GNAT  VALUtS 

C ACCUMuLATt  SUM  F UP  MEAN  CALCULAllON 

00  SO  1=1, M 

nUPFtk  INI IUI , IBUF 1 ,b,M, 1ST  AT  ) 

20  CALL  STATUS(IOl) 

GO  TO  ( £0 , 30 , 30 ) , 1ST  A T 
30  UU  «0  J = 1 , N 

vaL  =IhUFl  (J)-IZ-H 
A SUMS A SUM  + V AL 
1 V = I BUP 1 ( J ) 

all  IGKA(IV)  = lGwA(IV)  + l 
SO  CONTINUE 
C 

i ' compute  f-.tAN 

ASUMs(A3UN/AUIV)+12d. 

C 

l St t up  the  tekTkumx  uisro 

OU  oO  1 = 1,  2St> 

A X ( I ) =F LU A T (1-1) 
f>0  AY(  IJsFLOAT  ( IGKAII  ) ) 

nhITE (h,9W00) UX( 1), AY U ) ,1=1 ,2S6) 

•\Pl  rE(fe,«)POn)ASUM 
VS»00  FUWKAT  CS(2X,F10.3)  ) 
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call  inith^u) 

CALL  H1NITT 

CALL  TfcKM  (I,10?«) 

CALL  NPTS(<?56) 

C 

r i>ti  i*e  min  k max  fu«  thf  y-axis 

YMIN  c 10000.0 
YMAXS  -10000.0 
CALL  PiMMX(AY,  YM1N,YMAX) 

CALL  DLJMX(0.0»2b5.0) 

CALL  ULIMY(YMIN,YMAX) 

CALL  VbAKSl (6,b,3) 

CALL  CHECK ( AX , A Y ) 

CALL  USPLaY(AX, AY) 

c 

c rtKirt'  out  the  mean 

CALL  F.N.T  AKS(60,  30) 

CALL  ANMCJut 
>.Kllt(llK,yO?0)ASUM 
vo<?o  Fi.inmak  ( lHt,  'MEAN  = • , r 1 0 . o ) 

C 

L f.El  A HAKbCllPY 

CALL  flNPoTCIC) 

CALL  K I M I T T ( 0 » 0 ) 

kETUPn 

ENO 


SUHHOUflNt  OUMMr 
C OUMMY  SUd 

Xs  1 

KETURN 

cnO 
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SUHNUUTINK  DUMPSC 

COMMON/ I UA  U/IbUF  1 U096)  , IbUF2(*409b) 

DATA  ICR, ILK, IT TY/7,6, 1/ 

SUBROUTINE  PURPOSE  0U1PUT  AM  UPACKED  6-BIT 

P1CUTRE  TU  SCANNER 

READ  IiM  NtCESSARY  VAkIAPLES 
101  - DISK  KILE  NUMBER 

ITK1  - FILE  NO  UK  INPUT  IK  1 APE 
II  -NO  OF  LINES  TU  rtt  I kANSE EKKt 0 

JJ  - NU  UF  SPOT/LJNE  TU  HE  TRANSFERRED 

NF DPT  - 1ST  OATA  POINT  TO  bE  TRANSFERRED 
NRITS  - NUMhfck  UF  OIIS/PIXEL 

PEAU( ICk,SOuO)  IDl, ITF 1 , 1 I , JJ,NKDPT,NBI  IS 
C ADVANCE  TU  CORRECT  FILE  NU  IF  NECESSARY 

IF(ITFl,Lt.l)GU  Tu  3 
011  1 K = 2,ITK1 
1 LALL  AOF(IOl) 

3 CONTINUE 

00  100  1=1,11 

C GET  ThE  OATA  FKUM  THE  DISK 

C 

huFFEk  1N( lul , 1H0F 1,H,JJ,IS1A1) 

JO  CALL  ST  AT  US  ( 101 ) 

Gu  TU  (<J0, 30, 30), I ST  AT 
C 

C CHECK  FUN  h-BIT  Ok  8-hIT  OA  I A 

30  CONTINUE 

IFINHI TS.Eu.BlCALL  SlXHITlJJ) 

CALL  UT  SCAN  ( 1 HIJF  1 , NF  OPT  , J J , I ERR  ) 

1F(  IfcNK.NE  .0)NRj  Tfc(  IT  rY,9000)  It'Ktf 
100  CONTINUE 
C 

C PICTUNE  SHUULO  BE  PROCESSED 

WRITE (ITTY,9010) 

NE  T URN 

SOOO~FORMAT(SI 10) 

9000  FORMAT ( IX, 'ERROR  IN  0UTPU1  TO  SCANNER* , IS) 

9010  FORMAT! IX, 'PICTURE  PROCESSED') 

END 
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htJriWliuT  1 iMt  SIXHIUJJ) 

LJVwijh  / I i.i  a F A / JftUF  1 1409b) , IblJF<>(4096) 

I>Al«  1CM,  1LP,  I TT  1/1, 6,  1/ 

S'irtkUuTlNC  PURPl.St:  CO-MVthl  b HIT  l)A|A  Tu  b bIT  DATA  dv 

KhMUtflNG  TWO  LbB  (SHIFT  WIOHT  £ ) 

C 

l.  JJ  - NO  OF  PlXtLS/LlNB 

I 'O  10  1 = 1 » J J 

10  IrtUFl  (I)siBUF  1(1). SHIFT. •*£ 

9t  f IJPN 
F.wO 
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c 

c 

c 

c 

c 

c 

c 

c 

r 

c 

r. 

c 

c 

c 

c 

c 


SUBROUTINE  GNCOPY 

COMMON  /IdaTA/IBUF  1 ( u 0 9 e> ) , I8UF2(4  0 9t>) 
data  ich, ilp, i it r/7, 6,  i / 

SUmHOuTINE  PURPOSE:  DU4P  A FICIUKh  FkU*  uivt  DEVICE 

10  ANOTHER 

(DISK  ID  DISK)  (UlS*  Tu  I APF  ) 
( TAPE  Id  CAPE)  ( I APE  TO  DloKj 


« E A D IN  ixECtSSAkT  VAFIAhLES: 

101  - INPUT  DEVICE  NJMBtk 

ITF1  - FILE  ND  OF  1NPUI  iF  TAPt 

ID ? - OUTPUT  DEVICE  NUMBEk 

1 TF  2 FILE  NO  OF  UU1PUT  IF  (APE 

II  - NO  Li F LINES  TO  HE  TRANSFERRED 

JJ  - ivO  OF  SPOT/lInE  TU  HE  TRANSFERRED 

PEAO( I Ik, 8000) lul , I TF 1 , 102, ITF2,  II , JJ 

kEnIND  THE  uEvlCtS 

K fc  A I N U 1 D 1 


C ADVANCE  TO  CUkkECl  FILE  nu  IF  NECESSARY 

IFdlU.Lt.  1 )G0  TU  3 
DO  1 K = 2 » 1 TF  1 
1 CALL  ADf(lui) 

3 CONTINUE 

IF ( I TF2.LE. 1 )G0  TO  5 
DO  4 K=2 i I TF2 

4 (all  ADFCID2) 

5 C ON f I NUc 

DO  100  1=1,11 
C GET  the  Input  DaTa 

hUFFEK  1N( lul , IbUF 1,B,JJ, ISTAT) 

20  CaLL  STATUS(IOl) 

GO  10(20, 30, 30) , ISTAT 

30  Cunt i nue 

C KkllE  Irtt  DATA  10  lHE  0UTPU1  DEVICE 

BUFFER  OUT ( 1D2,IHuF l,h,JJ, ISTAT) 

So  Call  st  a ruse  id 2) 

Gu  10  (50, H 0,60), I S T A 1 
• uO  CONTINUE 
100  CONTINUE 


C _ DATA  TRANSFER  SHOULD  PE  L0*PLt1E 
~ End  FILE  ID2 
HtulNO  I D 1 
«RI TE( ITT Y, 9000) 

PE  T UkN 

9000  FORMAT  (u,  'COPY  CuMPLElED') 

ROOO  FORMAT  (Hi  10) 

* End 
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SoHnUUTINE  GPAMAP 

CUwMOtM/I0Al  A/IPUF  1 ( A 096)  , lFiuF2(  A09h) 

Dl^tNSlOr  IN( 1 Ou)  , 10(100),  I OUT  C 1 0(;0) 

OaTA  ICk, ILP, ITY/7,6, 1/ 

PunPOSE — TO  Pt*AP  ThE  GnAY  LEVEL  OF  A .yINOO*. 

V*wlAriLFS  TO  «t  PE  AO  I N HY  GWAMAP 

I D 1 - DEVICE  NO  OF  INPUT  PlLTuPt 

1 1 K I - FILE  NO  OF  Input  if  i fipp 

IU2  - DEVICE  NO  OF  UUT  °UT  PiCIUhE 

ITF?  FILE  NO  OF  OUTPUT  IF  1 APf 
IMaX  - MAXIMUM  PICTURE  VaLI'F 

II  - NO  OF  LINES  lu  hE  TPAfJSFEKPED 

JJ  - NO  OF  SPUTS/LIME  TO  Ht  MANSF'tKNED 

St  I UP  IMEwNAL  DHtNSIuMNG  l 1 M I IS 
MAXHPS=1 no 
MAXPTSslOOO 

KEAu  IN  USUAL  VAklAhLtS 

KtAD(  I CP,  8000)  ID1 , ITF1  f 11)2,  I1F2. 1 *AX,  1 1,  JJ 
IMAX=1MAX+ 1 

NUTE  THAT  TFiAX  HEKE  PEFFPS  TO  THE  TOlAL  NOMHEP 
OF  GPAY  LEVELS 

PLNlNU  Tnb  DEVICES 
WEnlNO  101 

C AOv/ANCE  lU  CuPPLCT  FILE  NU  if  necfssapy 

IF l ITE 1 .LE. 1 )G0  TO  3 
00  1 K*2,ITF1 
1 CALL  AOE(lDl) 

5 continue 

IF  ( I TF2.Lt. J >U0  TO  5 
DO  u N=2,ITE2 
*4  CALL  ADF  ( 102 ) 

5 CONTINUE 

C HEAD  IN  eWFAK  POINTS  OF  MAPPING  FUNCTION 

00  21  IHP=1 ,K AXBPS 
WtAU( ICP,8000)  11,12,13 
Ir ( I 3.Nt .0)  GO  TO  41 
lN(IBP)sIl+l 
21  IU( 1BP)=I2 
K«ITE ( ITY , 31  ) 

31  FOWMATC  KEAD  IN  TOO  MANY  GPAY-LEVFL  mAPPINu  hWtAKPTS- 
1 STOPP 1 Nu ' ) 

STOP 


\\\  AD> 

"il  4J 
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Jl  IHIM>.Lt..'iAxP!S)  GO  TO  6 
PRllEC ITY,S1  ) 

SI  MMtiAH*  MU  I ENUUGH  MEMORY  FUK  GKAY  LEVEL  * APR  1 i^G-S  1 OPP  1 NG  • 
S flip 

lhP=It>P-2 
1*0  9 1=1,IMaX 
DU  7 J:  1,1  HP 

7 IF  ( 1 .GE  . IM  J ) . AND.  I .Lt . I N ( J ♦ 1 ) ) GO  TO  11 
i%  R I r E ( I T Y , 8 ) 

h FORMAT  ( * (ME  GRAY  LtVEL  SPtC  I F I C A T I ON  IS  WRONG-STOPPING') 
STOP 

11  CONTINUE 

X = FLOAT ( i 10  ( Jfl )-IU( J) )*( I -I M( J) )) /FLOAT 

* (IN(J+1)-IN(J)J 

IF (X.Gt  .0)  10UT ( I)=IFlX(Xt.5+IU( J) ) 

IFIX.LT.0)  I UUT  ( l)=IFlX(Xf .499+IUCJ)  ) 

° CONTINUE 

C PtRFORM  Tut  MAPPING 

’ DO  100  1=1,11 

PUFFER  INI ID1 , IbUF 1 ,B, JJ, 1ST AI  ) _ 

?0  CALL  STATUS (10 1J 

GO  TU(20,30,30) , ISIA1 
30  CONTINUE 

DU  ao  J=1,JJ 
I V = I BUF  1 ( J 3 + 1 
4 0 I_HUF2(J)  = IUUr(lV) 

C WRITE'  OUT  TmE  mapped  P1C1UPE 

puffer  OUT (IU2, I HUF2 , R > J J , 1ST  AT) 

SO  CALL  STATUS (I D2) 

GU  TU  (50,100,100),1STAT 
IOC  LUNllNUt 

l close  FiLts 

End  FILE  1 U2 
REWIND  ID1 


RETURN 

"OOo  F URM A T (Ml  10) 
END 


Sill  wfnil  1 Nt  J I { f ► 

l U 1 ••  / I I.  A I A/  1 r.Uf  1 l 'JUPo  ) , 1 HUK  ? I '40'}r>  ) 

• A 1 a J l n,  ILP,  111  v/7,  t>,  1 / 


r 


t.  SuhhUOTInl  PUPPuSt.  - L>'t  A It  KaNDUH  L 1 .Mf-.  JiTTKP 

l 'OhmUU  1 1 1\  t v/AkIAi-LKS 

L 101  - UtVlCfc  nu  UK  Mic_  JIvPUi  PICTURE. 

C 1 i r 1 - MLt  NO  OF  INPUT  IF  l«Pt 

C 1 0 c?  - OfcVICt  NONnrP  Ut  OUTPUT  piCFUpF 

1 1IKP  K 1 L f NU  UK  OUTPUT  ]t  TAPE 

C )i  - Ni  I UK  L I M b 1 U t>F  iPAwijFtKPtu 

C JJ  - JO  UK  SPul  S/UNfc  IU  Ht  TP ANbKtPKfcO 

C iMi  - .S  I ANuAki)  UKVIAHjw 

L If  u — iA  i_  A i'< 

t.  ibrfcu  - SfcK.U  KUP  i>KN  GtUt. PAiup 

C 

c i •'  t i al  i a i ouk  i 

L'l 1 10  lslrUO^r) 

1 0 I miF  1 ( I ) =0 

c 

C rrjau  IN  f Pfc.  V AP I A^Lf  b 

p t A u ( ( (.  * , fi  0 0 0 ) I 0 1 , 1 T K 1 , 1 IV  , 1 1 1 ? , 1 1 , J j 
K K Al)  ( ICk  » "000  J l SO  , IHU  , 1st  1. 1) 

L 

C (uuVtrtl  MtAw  K bfu  Utv  I fl  KLDAIIM;  POINT 

Su=PlUAT ( ISu)/100o. 

X it  AN  = FLUAId^N) 
r ptxjNO  T pit  ufcvICFs 

Ph  .'i  I NU  101 

L a 0 v A Ni;  t TO  lOKPtCT  M L F NO  ]F  OF.CfcbSAKY 

IK  ( I TK  1 .Lfc  . 1 )UU  To  3 
I'iJ  1 * = 2, 1 IP  1 
1 l alL  aOF ( 1 01 ) 

3 C'JN  T 1 ImUF 

IK( I lP2.Lt . 1 )U0  TO  b 
00  0 x = c»Ilft' 
a CALL  AUMJU*) 

S LU'.TIUUF 

I OPtPATt  On  1 Ht  PlLTUPt 

»/*)  100  1 = 1,11 


L 

<*t.l  ThF.  If  POT  Ll-Jt 


hUF  F fcP  I N ( ] D 1 , l Hot  l , n , J J , 1 S I A I J 
«*(i  LALL  SMTOSUDlJ 

GO  TO  (<?0,,>S^S),  I.SI«I 
C t»F  T fl  GAUSSIAiJ  rtV  ai'U  St  1 I Jill 

e 5 CALL  h NON* (SAMP,  ISf.tl'»S.)»  XV.F.AN) 

I J 1 T 1 = I *-  IX  (SAMP+O.S) 

l. 

C CLEAN  Tut  OUTPUT  LINE  10  2fchO 

On  JO  J=1,JJ 
Ju  1 ^UF <?  ( J J = o 

IE  ( 1 A*S(  I Jll  T ) .GE..)J)GU  TJ  lu 
f P-JT  I*  1HE  J II  TEW 

IE  ( I J]  T 1 .oE.GHn  10  So 
u=i»«siunTj  + i 
•11=0 

Ui  uO  J=1J,JJ 
J l = J 1 ♦ 1 

AO  J huF  <?  I J 1 ) = I r>UF  1 ( J J 
Go  TU  Vy 
SO  LOivT  I l»UF 
IJsIJITT-U 
J 1 =0 
DO  hO 
J1=J1*1 

o 0 I Htrf-2  ( J ) = i HUF  1 ( J 1 ) 

VO  CONTINUE. 

c 

c put  our  the  liuf 

SUE  p h.  N OUT  ( 1 1)<?  # 1<3UF2,H,JJ,  IsTaT  ) 
7o  call  statost  io<;) 

C*i|  to  ( 70, 1 00,  1 00)  , IS  T AT 
too  CONTIhlUt 

L "*o  F I Lt  I Ug 
Kt>lND  1 0 1 
k t 1 UWN 

»‘OuO  MJKMAT(rtllO) 

E 00 


BEST  AVAIL/; 


lO1  r 

lb.: 


206 


* 


b ; h k i » y I i , > i • N i J I 6 E 

l i / it>A  I A / 1 1 Of  l ( <H'9r  ) , 1 HL‘f  ? I 4U  V*  ) 

i 1 .U  u II  n,  ILi  , L N V //.'»,  1 / 


r 

SUrkUOI INF 

VAklAHl.rb  kf  in  In  U iV  U I S t 

r 

i U 1 

- 

Dr  V 1 C f DO  Of  iNkOl  kl(  TOKf. 

r 

lin 

- 

f ILL  'id  uf  INPUT  If  I Akt 

t. 

i'V 

- 

UFVlCf  DU  Uf  UUTkUl  PlCIUkt 

c 

l Tf  x 

F 1 L fc  N'J  uf  DU  Tkij  | II-  | Akt 

i 

I 1 

- 

NU  uf  Li'VtS  TO  I'll  Tii  Ai'iSf  fckf’t  D 

c 

JJ 

- 

•\'l  I'f  Ski'T  /L  l Mh  1 I'  kc.  Tf'ANSFFk 

L 

J i*i  i»  * 

- 

hi  A X 1 1U‘M  1'  i l I Ukf  V A L t ■ F. 

L 

IS!) 

- 

b T Ai  DaKIi  l)t  V 1 A 1 I [)!■ 

f 

I N N 

- 

bkfcc  If  IED  kt  An 

c 

I St  to 

- 

b 1 1 |J  f'UK  UFN  lit  Df  k A 1 Or 

L 

ISL 

s 

()  - ADO  A Dlf  f f.kfNl  MV  ID  EACH 

C 

r 

= 

1 ~ A 1)0  A DlffhkFNT  KV  T Ci  tACn 

W 

c 

i n i al i zt 

I 

r Of  1 li  I kuf  ? 

b k U T 
LllvF. 


1 0 


00  lu  1 = 1,4046 

1 nUf  H I I =0 
1 Dllt  ? I 1 ) =0 


C ntAl)  IN  I hfc  VA^lfthLtS 

«tAO(  ICK,**000  J H' l , I 7 F 1 , lOd,  1 1 F <>,  II  , JJ 
"t  Aij  l ] Ck  , H0U0  J I max  , 1 SO  , I Vi  , 1 St  F i/ , I SI 

c 

C L'KVfc^I  MF  An  and  SIC)  DtV  I (j  Hu  a | I m;  HuIUl 

Si  • =F  LuA  f ( I So ) / 1 000 . 

> i"F  A N = r LO«  T ( I MW  i 

L "UlM1  INI  I T A NL  OUTPUT  DfVlCfcS 

k t ‘I  1 H D 1 u 1 

C mDvANLN  Tu  COkkfCT  F 1 L f HU  If-  Nfct.tSoAHY 

1M  J TM  .Lt.  1 )G0  Tu  3 
DO  1 M=r>,  1 IF  1 
I CALL  Al'F(IlU) 

3 COMI..Uf, 

If ( ITF^.Lfc. 1 )UO  TU  5 
DU  4 ksH,  I Tf? 
a L ALl  AUl-  ( I(j<?) 

S Luull.-in 

C UkfkATt  UN  1 Ht  kllTukt 

DU  1 0 1>  1 = 1,11 

( if  A l)lJ  I MU  Mil  TO  Tnf  LIDf,  Gt  T I ut  kv 

if  ( ISL.fu.  1 Jl  ALL  kwOki-.  ( ba^k,  I si  LI) , Su , A *<  fc  A n ) 
MJFF  t«  lMlul,  IhuF  1 , t) , J J , 1ST  AT) 

<-0  r.ALL  STATUS!  1 u 1 J 

1-U  TO  (<?0,io,30),lSlAT 
3>)  Ciiw  f lr.ut 


In  klLTUkt- 
JD  klCH'kt. 


’Best  available  copy 
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hi  SO  J = l,J.i 

( IF  AOulMO  x * TU  E»Cm  SPOT,  l/t  I I Ht  "V 

1M  ISL.fct..O}CALL  xMiki-!(SA.**F  , 1 St 1 1) , SO  , XF.t  A*  ) 
C AW>  1 r.  I He  uUlSt 

1 wt  F J ) = I FI  X (FLO  A f l leun  ( J ) ) «-SAKp-M.  .S) 

) M lUoF?(  J ) .LI  .0)  1 HuF  ? ( J ) =0 
li-  ( i»UF«;(j)  .bl  . l*iAX)  I<>uh?'(Jj  = I*Ax 
So  CowT  I i.oL 
t 

C «MI  I F HUT  T hF  tvdlsr  PlClukt 

t.sL+  t 

nljFF  fcx  1>U  I l 10?  / irtUF?/i?/jJ  / 1ST  <*  I ) 
so  CALL  STaTuSIIi)?) 

SO  TO  (oO/7ii,  70),  1ST  AT 
7 0 COW  f I wlJF. 

CALL  xMjKMSAmP,  1 sfc  F 0 , SI  • , X t A !. ) 

CALL  kMiM  ( SAi*'P  / ISFtO/SO/Xr-'FAS) 

CALL  XWOM’’  foANP,  1SFU  ,S!),  xMEAiO 
100  COF.  Tl'lUt 
C 

i'  HtnIM;  Ot  VICES 

two  F ILF.  lUP 
x t A 1 W 1 1 101 

i r!F.  Ujl"l 

> * 0 0 0 F Ok.  'AT(xJio) 

fcwi> 
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wrjjgg 


> T i n f ul  oUfc  <1 

L 1 1 ■’  i ll.\/  1 Om  I A / 1 fcwK  1 ( A 0 Vh  ) , 1 ~>UF  £’  ( <«  ) 

ii*iTa  1L«,  H h, n Y//,h,  1/ 

(.  I IFF  ►'UwHoSF  - lYKoHAUK  |)Y  KMUUlol,  o ) I,  J I 1 / A 1 i >l>v 

l V A « 1 A it  l.  k.  3 K t A 0 I Cl  1 0 1 1 fc  i» 


c 

I I 'M 

- OF  V ILF 

Ml  OK  j’lHljl  *•' 

T f 1 uKt 

c 

1 IF  1 

- FILL  JO  CiF  Ilo’UI  IF 

T mHi 

c 

1 'if 

- OF  v ICt 

.1(1  uF  1 1'  1 1 H 'J  1 

H i C 1 (J  M t 

r 

1 IF? 

F 1 ..  L M.) 

ok  nor hut  if  i 

-H  t 

i. 

1 1 

- NO  OF  LJ 

Mb  1,  I H i j 1 H 

l L T OKt 

c 

JJ 

- Nil  OF 

Son  (S/LI  fit  To 

HI  1 n ,, i'.  o F 1 * r h 

t: 

vt- 1 f 

- Ml)  (jF 

Oils  lu  O t S H 1 K T 1 1 ( * 1 1 1 H I ) 

c 

*\  t A 0 1 f \> 

lot  N K.  C F S S 

A k Y VAKjAflLrS 

<1  ail i n il,  Hf. (,<))  1 •>  1 , 1 1 f 1 , 10«M  J TF  ?,  1 1 , JJ  ,*\r  M 
1 F l f i h ] 1 . t u . 0 ) Nhj  f t ( 11  i *^000)Nol I 
ll-  ( Nhl  1 . I>  1 . 7 ).•  k 1 Tfc  ( 1 T Y , WIKJOJN*  IT 


C *'  c «v  1 l Till-.  OtVlLtS 

‘'c  i‘  1 i»i>  I l)  1 

L «jv  ANl.F  In  LUKNi-CT  f-  I L K.  ml  IF  Mr  t,  t S S A K Y 

IF  Cl  TF1  .Ll  .1  ) |»U  ru  5 

ii. i 1 r\  = <^,  lIF  1 

1 ( aLL  aDF(irn) 

5 CuimI  iUUt 

if  1 1 1 1 a.LF  . Ubll  To  s 
U,i  a A- f , \ \ r £ 

'4  CALL  A l ) r l III?) 

S uiillMUh 

r. 

C OtbrAuF  THF.  F I L 1 ll*t 

I0lV  = ,'**:'IOl T 
Hfi  10u  1 = 1,11 

•OJFFFR  I l ( I 0 1 , IrtUF  i IS1  »T  I 

?!'  i:  ILL  STa  I U S ( 11)1  ) 

mi  in  (<;  0,4  o,3U), IS  if  i 
in  t Mo  I 1 out 


f 

c ^lir  lull  1 r«F  F>uF>E*J 

n jF  F F w UUl  ( J Ot?t  loT»>l  l 

nil  CALL  i>7  A J J3  l l Ur*  ) 

uo  Tu(6u,  1 00,  1 0 0 ) , J S I mT 
Mill  r.'tul  l'iut 
L 

C C l » ) t OUT  T h F.  F ] L fc  S 

t'4(>  FRt  10.? 

*£.«If'0  I 0 1 

n£  1 ii^ri 

h 0 (j 0 F‘ikv»k  T ( ft  1 1 0 ) 

iJiliWI  h ' J K i-'l  A I l 1 X f ’ S H 1 F T VAl.  iiL  •KJnIx  •% ' ■ I I = 1 t 1 *1  I 

t -M  D 


' ' )r  "IM  1. 1 'Nt  ►>  F 41.  3 

!.•!  * i*,trw  I Of  1 a / I r ( ‘MiVfO  , i n jk  < i '.OVb) 

■ilt  ILK,JLh/7,r/ 

L MjaMio  I l USE  I.*  l1'!  REDALS 

l 1 0 1 - Ut.VlCt  l J l«t-  ilxFUI  FlCIU^t 

t.  I'M  - FILE  .•>•0  OF  INPUI  JF  TAPE 

• lu?  - OLvlCt  '<0  OF  OUTPUT  ►'iClOKt; 

i He*’  F iLF  Mj  OF  • »U  I P 0 1 IF  lAPt 

L II  - M>  UF  Li^tS  III  I PaNSF t KKfc 0 

l Jj  * Ml  UF  SKUIS/UNt  Ht  1 KAWSiF  tWrtt’O 

C 1 3 1 Zb  - HUvtt*  LF  2 PfcUUCTlOF.S  I It  SJ7t 

(.  IMAX  - MAX  MUM  HILIUWt  VaL'JF 


f I m 1 I a i.  1 Z F I OOF  1 t»  MuF2 

I'll  10  I s 1 , 41.^6 
MUF  1 1 1 ) = 0 
1"  MJF2(I)sy 
r. 

C KtAL>  M lUt  VaklanLLS 

l<F.AD(  1M,  *000)  Ml  , I IF  1 , 102,  1TFP,  1 1 , JJ,  ISI/t 

IF ( ISIZfc.Lh.OjRfeTUkK 

C i-c  I Int  COHPfcLT  ^0*  A.xu  COLUMN  OFFSETS 

1-OVsIF  iX(FL0mT  ( 1 i )*  ( 1 1.  /FLOAT  (ISIZE))/p.)-i 

iC'JL=lFIMFU'Al  l JJ)M  1 1.  / F L t>  A T ( ISJ/fc)  )/?.)-! 

1 SZ  = ISI /t  * *«• 

t rfC'lMU  THE  INPUT  A :Mi>  IJUTPUT  OE  V 1CE3 

H't^MO  10  1 

l.  ui-'VANCt  Tu  LUMWfeCT  FILE  MJ  IF  F'ECEbSxrfi 

IF ( 1 1 F 1 .LF . 1 )bU  ru  3 
(HI  1 KSZ,  1 I F 1 
1 CALL  AUF(Hl) 

3 (.  UNI  I NUt 

! F ( I TF2.U  . 1 )i»U  TO  S 
I'U  4 K=P,lTFP 
H l.xLL  AOF  ( I OP  ) 

S LUMlUUF. 

0 

• UPf-.iATfc  Ou  I Ob  FlCiUPt 

( 

I..  l «3t«r  I 111  fl'P  PUHDFrt 

l Si 

. iJ  1 Op  1 = 1#  l"U*v 
L=L+  1 


10^  CONTINUE 

OU  100  I=1,II,ISIZE 

W112  H=1,JJ 

13?  I BUF  ? ( 1 3 ) =0 

do  io7  tiii=i,isize 

HUF  F EH  1N(ID1,I6UF1,B,JJ, ISTAT) 
20  CALL  STaTUS(IDI) 

GO  TO  (20, 30, 30) , 1ST  AT 
TITCDnTTNUE 
*=ICOL 

" DO  143  1 3=  1 , JJ  , I SI  ZE 

K =K  ♦ 1 

I3MaX=I3*1SIZE-1 
00  143  14=13, I3MAX 
RTTRTJr2CK)  = I8oF2(K)  + lHUFl  (14) 

107  CONTINUE 

"L=LM 

K = I COL 

DO  144  IS=1,  JJ,  1SIZE 
K = K ♦ 1 

— pnr  TB0F?CK)=IBUF2(K)/ISZ 

c 

C WRITE  OUT  TmE  NOISY  PICTURE 

BUFFER  0U1  ( I D2 , 1 BUF«? , H , J J , 1ST  AT  ) 
60  CALL  STATUS(ID2) 

GO  TO  (60,70,70), ISTAT 
“ TIT  CONTINUE  * ' 

100  CONTINUE 

c 

C INSERT  THE  oOTTUM  BORDER 

DO  104  1=1, JJ 

1 04  I BUF  2 ( 1 ) =0 

rrcr  lire  r*i,n 

BUFFER  0UT(ID2,1BUF2,h, JJ, ISTAT) 

106  CALL  STATUSCIO?) 

GO  TO  (106, 105, 105), ISTAT 

105  CONTINUE 
C 

C FFwTW  THE  DEVICES 

END  FILE  102 

'**£*  I NO'  1 D 1 

RETURN 

8000  FORMAT (81 10) 

ENO 
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subroutine  bubals 

common/  1 OAT  A/  IBuF  I * I HUF  2 ( U 0P6  ) 

TiATA  ICf,1LP/7,6/ 

r SUBROUTINE:  VARIABLES  USED  IN  HGRALS 

ID1  - OtVlCt  NU  UK  INPUI  PICT URE 

1 TK 1 - KILE  NO  OF  INPUT  IF  I APE 

102  - DEVICt  NO  UF  OUTPUT  PICTURE 

l TK  2 KILE  NO  UF  OUTPUT  IF  TAPE 
II  - NO  OF  LlNtS  TU  TRANSFERRED 

JJ  - NO  UF  SPOTS/LINE  TU  BE  TRANSFERRED 

IS1ZE  - PU*Er  OF  2 REOuCTIuNS  IN  SIZE 
IR.AX  - MAXIMUM  PICTURE  VALUE 


C INTIALIZE  IBUF1  * I BUK2 

00  10  I=l,«096 
I BUF 1 ( I ) =0 
10  lBUF2(I)=0 
C 

c read  in  the  variables 

~RF#OCICR,bOOO) IDI , ITF1,IU2,ITK2,I1,JJ,1SIZE 


PUFFER  lN(IDl,IBUFl,b,JJ, ISTAT) 
IQ  CALL  STATUS(IOl) 

CO  TO  (?0,30, 30), ISTAT 
3o  CONTINUE 
f 07  CONTINUE 
K=ICOL 

T>0  SO  J=1 , JJ, IS1ZE 
C SAMPLE  THE  Image 

k = k + 1 

IbUF?(K)=IbUFl (J) 

SO  CONI INUE 
L =L  ♦ 1 

C WRITE  OUT  TmE  NUISY  PICTURE 

“SOFTER  OUT  ( I D£  > lBUF*?,b,JJ,  ISTaT) 
60  CALL  STATUS (I D2) 

“CO  TO  (60, 70, 70) , ISTAT 
70  CONTINUE 

— nrcrroNTiNUE 

c 

c insert  the  bottom  buroer 

du  ioa  im,jj 

'100  1BUF2(I)=0 

l>0  10S  I=L,II 

fluFFTR  UUT(102,T6UF2,B,JJ, ISTaT) 
106  CALL  ST  A TUS ( 1 02 ) 

CO  TO  ( 106, 105, 105) , ISTAT  ' 

105  CONTINUE 

C ' 

C REWIND  THE  DEVICES 

END 'FILE  ID? 

KErtINO  101 

RETURN 

BOOO  F URM A T (8110) 

END 
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SECTION  C 


PROGRAM  FLOWCHARTS 
PARCON 


5 


DITS  MAIN 


217 


TODWIN 


0 


^ BEGIN  J 


3 


INITIALIZE 
IBUF1  & 
IBUF1  TO 
ZEROES 


READ  M,N,MM,NN 
1 1 ,JJ  , ID1 , ITF1 , 
ID2.ITF2 


REWIND  ID1 
REWIND  ID2 


IWORD= (NN+JJ )- 1 
I L I N E = ( MM  + 1 1 )-  1 

3E 

ADVANCE 
TO  CORRECT 
OUTPUT 
FILE 

* 


ADVANCE 
TO  CORRECT 
OUTPUT 
FILE 


READ  INPUT 
RECORDS  UNTIL 
DESIRED  RECORD 
IS  REACHED 


Sheet  2 


DUMPSC 


t 

NOTE: 


i 


DUMPTV 


This  subroutine  is  not  implemented  at  the  present  time. 


Sheet  8 
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JITTER 


Sheet  12 
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DO  60  J=IJ,JJ 


Sheet  13 


I BUF2 ( J ) 


IHU 


END  FILE 
ID2 

REWIND  ID1 
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HGHALS 


Sheet  17 
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7±E= 

/ WRITE 
/ OUTPUT 
/ RECORD 

<§> 


END  FILE 
ID2 

REWIND 

ID1 

X 

( EX'T  ) 
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REDALS 
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APPENDIX  B 

SUBJECTIVE  EFFECTS  OF 
SELECTED  DIGITAL  DEGRADATIONS 


APPENDIX  B 


SUBJECTIVE  EFFECTS  OF 
SELECTED  DIGITAL  DEGRADATIONS 

Section  IV-B  gives  a preliminary  subjective  image  test  set. 
In  asking  a trained  photo interprecator  to  evaluate  the  effects 
of  the  various  digital  degradations  selected,  the  prominent 
response  was  verbal.  These  verbal  comments,  recorded  in  Tables 
B-l  through  B-  23 , are  indicative  of  the  subjective  effects  noted. 


238 


* 


TABLE  B-l.  ADDITIVE  GAUSSIAN  SPOT  NOISE 
APPLIED  TO  A HIGH  RESOLUTION  B/W 
VERTICAL  DIGITIZED  PHOTOGRAPH 


Standard 

Deviation  = 

.3%  of 

Gray 

Level 

Range 

Good 

detail , good 

shadows 

Standard 

Deviation  = 

.6%  of 

Gray 

Level 

Ranqe 

Good 

deta i 1 , good 

shadows 

Standard 

Deviation  = 

1.2%  of 

Gray 

Level 

Ranqe 

Good 

detail,  good 

shadows 

Standa  rd 

Deviation  = 

2.4%  of 

_Gra£ 

Level 

Ranqe 

Some  loss  in  detail  of  KC - 135  tail  section,  probably 
caused  by  graininess. 


NOTE:  Loss  of  tonal  differences  between  B - 5 2 (painted 

A/C)  and  background  may  cause  the  interpreter 
some  difficulty. 

Shadow  detail  is  good  on  all  image  scenes. 

Standard  Deviation  = 4.88  of  Gray  Level  Range 

Graininess  causes  some  loss  detail  along  edges  of  targets. 
Tonal  difference  between  B - 5 2 and  portion  of  it's  background 
appears  to  be  enhanced  by  this  graininess. 
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TABLE  B-l.  ADDITIVE  GAUSSIAN  SPOT  NOISE 
APPLIED  TO  A HIGH  RESOLUTION  B/W 
VERTICAL  DIGITIZED  PHOTOGRAPH  (Continued) 


Standard  Deviation  = 9.77%  of  Gray  Level  Range 

Detail  of  tail  sections  on  both  aircraft  affected  by 
graininess.  There  is  a loss  in  target  background  tone  that 
is  caused  by  this  "salt  and  pepper-like"  appearance. 

Detail  of  left  tip  of  KC - 135  wing  completely  lost.  Men- 
suration would  be  difficult. 

Standard  Deviation  = 19.52%  of  Gray  Level  Range 

Detail  on  KC - 135  only  apparent  on  right  wing,  portion  of 
left  wing  and  fuselage.  Detail  of  B - 5 2 is  completely  lost 
under  magnification.  Loss  of  background  detail  caused  by  the 
"salt  and  pepper-like"  appearance.  Accurate  target  men- 
suration would  be  impossible. 

Standard  Deviation  = 39.04%  of  Gray  Level  Range 

Almost  complete  loss  of  detail  under  magnification. 


> 
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TABLE  B- 2 . ADDITIVE  GAUSSIAN  LINE  NOISE  APPLIED 
TO  A HIGH  RESOLUTION  B/W  VERTICAL 
DIGITIZED  PHOTOGRAPH 


Standard  Deviation  = .3%  of  Gray  Level  Range 

Edges  as  sharp  as  original,  good  shadow,  good  tonal 
differences. 

Stand a rd  D eviation  = .6%  of  Gray  Level  Range 

Tone  differences  not  as  great  on  KC-135,  primarily 
between  nose  of  A/C  and  fuselage. 

B - 5 2 fuselage  blending  in  with  hard  stand  parking  ramp  - 
poor  tonal  differences  causing  lack  of  target  detail. 

Overall  mesh  effect. 

St andard  Deviation  = 1.2%  of  Gray  Level  Range 

Shadow  detail  not  as  sharp  as  on  original  image  - 
difference  in  image  tone  on  KC-135  not  as  great  as  on 
ori gi nal  image . 

Standard  Deviation  =•  2.44%  of  Gray  Level  Range 

Many  "noise  lines"  across  target  causes  some  "break-up" 
of  image  and  loss  of  detail. 

Loss  of  target  detail  due  to  absence  of  tonal  differences 
between  B-52  and  hard  stand. 

Standard  D e v i ati o n 4 . 88%  of  Gray  Level  Range 

Noise  apparent  on  overall  image  scene,  horizontal 
"noise  lines"  break-up  target,  especially  left  edge  of 
KC-135  wing  - mensuration  of  wing-span  would  prove 
difficult  - same  applies  to  B-52,  left  edge  of  wing 
almost  lost  due  to  lack  of  target  and  background  tonal 
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TABLE  B-2.  ADDITIVE  GAUSSIAN  LINE  NOISE  APPLIED 
TO  A HIGH  RESOLUTION  B/W  VERTICAL 
DIGITIZED  PHOTOGRAPH  (Continued) 

differences  and  noise.  Poor  detail  in  tail  section  of 
B - 5 2 because  of  noise  and  lack  of  tonal  differences. 

Standard  Deviation  = 8.77%  of  Gray  Level  Range 

Excessive  noise  breaks-up  target  - causing  lack  of 
detail  - wide  "noi s^-bands " very  disturbing.  Almost 
complete  lack  of  detail  in  left  wing  edge  of  B-52  and 
tail  section. 

Standard  Deviation  = 19.55%  of  Gray  Level  Range 

"Wide  black  and  white  noise  bands"  breaks-up  target 
detail.  Target  mensuration  would  be  extremely  difficult 
if  not  impossible. 

Sta ndard  Deviation  = 39.  % of  Gray  Level  Range 

Noise  level  causes  greater  loss  of  image  detail  than 
image  scene  above. 

NOTE:  Adequate  shadow  detail  is  present  on  each  image 

scene . 
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TABLE  B-3.  GAUSSIAN  LINE  JITTER  APPLIED 
TO  A HIGH  RESOLUTION  B/W  VERTICAL 
DIGITIZED  PHOTOGRAPH 


Standard 

Deviation  = 

.05 

Pixels 

Standa  rd 

Deviation  = 

. 1 

Pixels 

Standard 

Deviation  = 

.2 

Pixels 

Edges 

not  as  sharp  as 

on  other  images 

Standard 

Deviation  = 

.39 

Pixels 

Excessive  "jitter"  causing  "break-up"  in  edges  of  target 
and  target  shadow.  Detail  of  KC - 135  tail-section  almost 
completely  lost. 

Detail  of  back  section  of  fuselage  and  tail  section  of 
B - 5 2 almost  completely  lost. 

This  "jitter"  is  quite  apparent  or.  vertical  painted 
lines  on  hard  stand. 

Standard  D e v i at  ion  = .78  Pixel s 

Excessive  jitter  breaking  up  target  and  target  shadow  - 
Detail  lost  on  nose  section  and  back  part  of  fuselage  and 
tail  of  B-52. 

Standard  Deviation  = 1 . 56  Pixels 

No  target  detail  apparent  - except  on  back  edge  of 
KC - 135  wings  and  back  edge  of  B-52 ' s left  wing.  Greater 
loss  of  detail  apparent  on  painted  B-52  vs  silver  (painted 
KC- 1 35  ) . 
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TABLE  B-3.  GAUSSIAN  LINE  JITTER  APPLIED 
TO  A HIGH  RESOLUTION  B/W  VERTICAL 
DIGITIZED  PHOTOGRAPH  (Continued) 


Standard  Deviation  _=  3.  13  P i xels 

Detail  of  B - 5 2 completely  lost.  Edges  of  KC - 135 
broken-up  to  the  point  where  all  target  detail  is  lost 
except  in  wing  sections. 

Standard  Deviation  = 6.25  P i xels 

100%  loss  of  target  detail  and  shadow  - shape  is  still 
apparent. 


NOTE:  Shape  of  target  remains  apparent  on  each  image  scene, 

particularly  when  no  magnification  is  used.  There  is 
poor  contrast  between  target  and  target  background 
making  accurate  interpretation  difficult. 
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TABLE  B-4.  VARYING  THE  SAMPLING  OF  A 
HIGH  RESOLUTION  B/W  VERTICAL  DIGITIZED 
PHOTOGRAPH  UNDER  HIGH  ALIASING 


Sampling  Reduction  = 2 

Borne  loss  of  detail  in  nose  section  (radar)  of  B - 5 2 . 
Sampling  Reduction  = 3 

Some  loss  of  detail  in  tail  section  of  KC - 135  and 
nose  of  B - 5 2 . 

Sampling  Reduction  = 4 

Loss  of  detail  in  nose  section,  tail  section  and 
engines  of  KC-135  - overall  target  image  lacks  sharpness. 
Much  of  B-52  (especially  right  wing  tip)  " bl ends- i nto " 
target  background.  Target  image  lacks  sharpness. 

Sampling  Reduction  = 5 

Edge  of  left  wing  KC-135  very  "choppy".  Engines  on 
that  wing  not  very  apparent  - very  little  detail  in  tail 
section  of  aircraft.  Target  shadow  and  background  tone 
differences  are  not  as  good  as  on  original  image  scene. 
Poor  detail  in  B-52.  Target  tends  to  blend  in  with 
background . 
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TABLE  B-5.  VARYING  THE  SAMPLING  OF  A HIGH  RESOLUTION 
B/W  VERTICAL  DIGITIZED  PHOTOGRAPH  UNDER 
REDUCED  ALIASING 

Sampl i ng  Reduc t i on  = 2 

Good  detail,  good  shadow,  good  tonal  differences 
between  target  and  background. 

Sampling  Reduction  = 3 

Overall  image  (especially  KC-135  appears  a bit  fuzzy) 
Sampling  Reduction  = 4 

Targets  do  not  appear  sharp,  loss  of  detail  in  nose 
sections  and  tail  sectionsof  A/C.  Wing-tips  of  B-52 
tend  to  "blend-in"  with  background. 

Sampling  Reduction  = 5 

Targets  lack  overall  sharpness,  much  detail  lost. 
Right  engine  of  KC-135  not  apparent  under  magnification. 
Right  wing  of  B-52  blending  in  with  background. 
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TABLE  B-6.  ADDITIVE  GAUSSIAN  SPOT  NOISE 
APPLIED  TO  A LOW  RESOLUTION 
B/W  VERTICAL  DIGITIZED  PHOTOGRAPH 

Standard  Deviation  = .3%  of  Gray  Level  Range 

Standard  Deviation  j= .6%  of  Gray  Level  Range 

Some  detail  is  lost  in  targets.  Second  aircraft  from 
the  left  appears  to  be  a "delta  wing"  aircraft,  rather 
than  a swept  wing. 

Stan d ard  Devi ati on  = 1.22%  of  Gray  Level  Range 

SJ: andard  Deviation  = 2.44%  of  Gray  Level  Range 

Very  "spotty"  - loss  of  edge  detail  because  of  spot- 
tiness. Target  background  (aircraft  parking  area)  tone 
affected  by  "salt  and  pepper-like"  interference  taking 
away  from  the  target  itself,  by  " bl endi ng- i n " with  the 
ta  rge  t . 

St  a ndard  Devi ati o n = 4.88%  of  Gray  Level  Range 

Aircraft  are  still  detectable,  however,  they  are  not 
identifiable.  Edges  of  hard  stand  lack  any  detail  or 
sharpness.  Area  between  third  and  fourth  hard  stand 
appears  to  contain  an  aircraft  (false  alarm). 

Standard  Deviation  = 9.77%  of  Gray  Level  Range 

Spot  noise  effects  entire  image  scene.  "False  alarm" 
detected  still  apparent.  Target  can  still  be  interpreted 
as  an  A/C  parking  area. 
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TABLE  B-6.  ADDITIVE  GAUSSIAN  SPOT  NOISE 
APPLIED  TO  A LOW  RESOLUTION 
B/W  VERTICAL  DIGITIZED  PHOTOGRAPH  (Continued) 


Standard  Deviation  = 19.51%  of  Gray  Level  Range 

Area  can  still  be  interpreted  as  an  A/C  parking  area. 
No  other  information  can  be  obtained  from  this  scene. 

Standard  Deviation  = 39.04%  of  Gray  Level  Range 
Target  area  not  identifiable. 
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TABLE  B-7.  ADDITIVE  GAUSSIAN  LINE  NOISE  APPLIED 
TO  A LOW  RESOLUTION  B/W  VERTICAL 
DIGITIZED  PHOTOGRAPH 


Standard  Deviation  = .3%  of  Gray  Level  Range 
Standard  Deviation  = .6%  of  Gray  Level  Range 

Standard  Deviation  = 1.22%  of  Gray  Level  Range 

Line  noise  effecting  (distorting)  shape  of  the  two 
A/C  on  the  third  and  fourth  hard  stand. 

Standard  Deviation  = 2.44%  of  Gray  Level  Range 

Line  noise  effecting  shape  of  all  parked  aircraft. 

Standard  Deviation  « 4.88%  of  Gray  Level  Range 
Aircraft  are  still  identifiable. 

Standard  Deviation  = 9.77%  of  Gray  Level  Range 

Aircraft  are  detectable;  however,  the  "wi de-bands " 
or  line  noise  is  "breaking-up"  the  target  and  making 
interpretation  difficult. 

Standard  Deviation  = 19.52%  of  Gray  Level  Range 

All  aircraft  are  detectable.  None  are  identifiable. 

St a nda rd  Deviation  = 39.04%  of  Gray  Level  Range 

It  is  still  apparent  that  the  target  is  an  aircraft 
parking  area.  Two  aircraft  are  detectable,  but  not 
identifiable. 

♦ 
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TABLE  B-8.  GAUSSIAN  LINE  JITTER  APPLIED 
TO  A LOW  RESOLUTION  B/W  VERTICAL 
DIGITIZED  PHOTOGRAPH 

Standard  Deviation  = .05  Pixels 

Standard  Deviation  = .1  Pixels 

Standard  Deviation  = .2  Pixels 

Edges  of  fourth  aircraft  from  the  left  are  broken. 

Standard  Deviation  = .39  Pixels 

Standard  Deviation  = .78  Pixels 

Aircraft  shape  distorted  by  "jitter".  Edges  of 
taxi-ways  and  hard  stands  are  also  distorted. 


Standard  Devia ti  on  = 1.56  Pixels 


Hard  stand  still 

identifiable. 

Aircraft  can 

be 

detected 

Standard  Deviation 

= 3.13  Pixels 

Hard  stand  still 

i den  ti f i abl e . 

Aircraft  can 

be 

detected 

Standard  Deviation 

= 6.25  Pixels 

Target  not  identifiable. 
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TABLE  B-9.  VARYING  THE  QUANTIZATION  OF  A 
LOW  RESOLUTION  B/W  VERTICAL 
DIGITIZED  PHOTOGRAPH 


Number  of  Gray  Levels  = 32 

Poor  contrast  between  target  and  target  background 
make  identification  difficult.  All  four  aircraft  are 
detectable,  however. 

Number  of  Gray  Levels  = 16 

Many  of  the  target  components  "blend  into"  their  back- 
ground. Identification  of  the  aircraft  not  possible. 

They  are  detectable. 

Number  of  Gray  Levels  = 8 

Aircraft  not  detectable.  Target  can  be  identified  as 
an  aircraft  parking  ramp. 

Number  of  Gray  Levels  = 4 

No  detail  in  image;  however,  the  overall  shape  of  the 
target  area  might  allow  an  interpreter  to  make  an  accurate 
i nterpretati on . 
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TABLE  B-10.  VARYING  THE  SAMPLING  OF  A LOW 
RESOLUTION  B/W  VERTICAL  DIGITIZED 
PHOTOGRAPH  UNDER  HIGH  ALIASING 


Sampling  Reduction  = 2 

Target  area  is  identifiable,  and  aircraft  are  detectable 
Sampling  Reduction  = 3 

Target  area  is  identifiable,  and  aircraft  are  detectable 
Sampling  Reduction  = 4 

Target  area  is  identifiable.  Loss  of  detail  does  not 
allow  for  the  detection  of  the  aircraft. 

Sampling  Reduction  = 5 

No  detail  in  image.  Shape  of  A/C  parking  area  is  still 
apparent . 


TABLE  B-ll.  VARYING  THE  SAMPLING  OF  A LOW 
RESOLUTION  B/W  VERTICAL  DIGITIZED 
PHOTOGRAPH  UNDER  REDUCED  ALIASING 


Sampling  Reduction  = 2 

Target  area  identifiable.  Aircraft  are  still 
detectable . 

Sampling  Redu ction  = 3 

Target  area  identifiable.  Two  aircraft  detectable. 
Sampling  Reduction  = 4 

Target  area  identifiable.  Aircraft  are  not  detectable. 
Sampling  Reduction  = 5 

Parking  area  still  maintains  shape.  No  target  detail 
apparent.  No  A/C  detectable. 
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TABLE  B-12.  ADDITIVE  GAUSSIAN  SPOT  NOISE  APPLIED 
TO  AN  INFRARED  IMAGE  OF  PARKED  AIRCRAFT 


Standard  Deviation  = 3%  of  Gray  Level  Range 
Edges  broken,  detail  remains. 

Standard  Deviation  = .6%  of  Gray  Level  Range 
Edges  broken,  detail  remains. 

Standard  Deviation  = 1.22%  of  Gray  Level  Range 
Edges  broken,  detail  remains. 

Standard  Deviation  = 2.44%  of  Gray  Level  Range 
Edges  broken,  detail  remains. 

Standard  Deviation  = 4.88%  of  Gray  Level  Range 

Edges  not  sharp,  some  loss  in  detail  especially  on 
boom  (located  rear  of  tail  section)  and  tail  section 
i tsel f . 

Standard  Deviation  = 9.77%  of  Gray  Level  Range 

"Salt  and  pepper"  effect  takes  away  from  target 
sharpness  and  detail.  The  right-rear  portion  of  the 
fuselage  blends  in  with  the  background  because  of  "spotty" 
effect.  Tip  of  right  wing  appears  to  be  rounded  rather 
than  straight  and  sharp. 

Standard  Deviation  = 19.52%  of  Gray  Level  Range 

Greater  "spotty"  or  "salt  and  pepper-like"  effects 
cause  loss  of  target  detail  and  edge  sharpness.  Boom 
completely  lost  under  magnification.  Right  wing-tip 
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TABLE  B-12.  ADDITIVE  GAUSSIAN  SPOT  NOISE  APPLIED 
TO  AN  INFRARED  IMAGE  OF  PARKED  AIRCRAFT  (Continued) 


appears  to  be  rounded.  Much  of  target  blends  into 
background  when  viewed  under  magn i f i ca t i on . 

Standard  Deviation  = 39.04%  of  Gray  Level  Range 

Loss  of  all  target  detail  under  magnification, 
appears  almost  like  a half-tone  newspaper  picture. 


the 


Image 
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TABLE  B-13.  ADDITIVE  GAUSSIAN  LINE  NOISE 
APPLIED  TO  AN  INFRARED  IMAGE  OF  PARKED  AIRCRAFT 


Standard  Deviation  = .3%  of  Gray  Level  Range 


Standard 

Deviation  = .6%  of 

Gray  Level  Range 

Standard 

Deviation  = 1 . 22% 

of 

Gray 

Le  ve  1 

Range 

Standard 

Deviation  = 2.44% 

of 

Gray 

Level 

Range 

Standard 

Deviation  = 4 . 88% 

of 

Gray 

Level 

Range 

Some  loss  of  detail  on  in-board  engine,  right  wing. 
Edges  of  target  very  "choppy"  due  to  wide-bands  or  lines. 

Standard  Deviation  = 9.77%  of  Gray  Level  Range 

Some  loss  of  detail  on  in-board  engine,  right  wing. 
Edges  of  target  very  "choppy"  due  to  wide  bands  or  lines. 

Standard  Deviation  = 19.52%  of  Gray  Level  Range 

Some  loss  of  detail  on  in-board  engine  right  wing. 

Edges  of  target  very  "choppy"  due  to  wide  bands  or  lines. 
More  loss  of  detail  due  to  wide  bands  or  lines.  Refueling 
boom  almost  completely  lost  under  magnification. 


Standard  Deviation  = 39.04%  of  Gray  Level  Range 

Rear  section  of  fuselage  lost  because  of  "banding" 
and  lines.  Boom  completely  lost  under  magnification. 


U- 
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TABLE  B-14.  GAUSSIAN  LINE  JITTER  APPLIED  TO 
AN  INFRARED  IMAGE  OF  PARKED  AIRCRAFT 


Standard 

Deviation  = 

.05  Pixels 

Standa  rd 

Devia  ti on  = 

.1  Pixels 

Standard 

Deviation  = 

.2  Pixels 

Standard 

Deviation  = 

.39  Pixels 

Standard 

Deviation  = 

.78  Pixels 

Target  edges  "broken"  by  jitter.  This  is  apparent  on 
all  target  components  except  the  rear  edge  of  the  right 
wing. 

Standard  Devi  a t i on  _= 1 .56  Pixels 

Lack  of  detail  caused  by  "broken"  edges. 

Standard  Deviation  = 3.13  Pixels 

Target  lacking  detail.  No  edge  sharpness  at  all. 
Complete  loss  of  target  components,  i.e.,  engines,  boom. 

Standard  Deviation  = 6.25  Pixels 

Complete  loss  of  target  detail.  Shape  could  be 
interpreted  as  a "thermal  shadow." 
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TABLE  B-15.  VARYING  THE  QUANTIZATION  OF  AN 
INFRARED  IMAGE  OF  PARKED  AIRCRAFT 


Number  of  Gray  Levels  = 32 
Number  of  Gray  Levels  = 16 

Light  gray  border  along  edge  of  target.  This  does  not 
effect  the  detail  too  badly. 

Number  of  Gray  Levels  = 8 

Gray  border  along  edges  of  target  distort  the  detail. 
Engines  on  right  wing  appear  to  be  separated  from  the 
wing.  Boom  located  in  tail  section  is  lost.  Shape  of 
rear  portion  of  fuselage  is  distorted. 

Number  of  Gray  Levels  = 4 

Step-like  edges  in  tail  section  of  A/C.  Target  lacks 
detail  in  general.  Front  portion  of  fuselage  blends  in 
with  background  ("smearing  effect"). 
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TABLE  B-16.  VARYING  THE  SAMPLING  OF  AN 
INFRARED  IMAGE  OF  PARKED  AIRCRAFT 
UNDER  HIGH  ALIASING 

Sampling  Reduction  = 2 

Poorer  target  - background  contrast  than  other  three 
images,  especially  in  tail  section. 

Sampling  Reduction  = 3 

Better  contrast.  Good  tone  differences  between 
target  and  target  background. 

Sampling  Reduction  = 4 

Good  tone  differences  between  target  and  target  back- 
ground. 

Samp  1 i n g Reduction  = 5 

Good  tone  differences  between  target  and  target  back- 
. ground. 
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TABLE  B- 1 7 . VARYING  THE  SAMPLING  OF  AN 
INFRARED  IMAGE  OF  PARKED  AIRCRAFT 
UNDER  REDUCED  ALIASING 


Sampling  Reduction  = 2 

Poorer  target  - background  contrast  than  other  three 
image  scenes,  especially  in  tail  section  of  A/C. 


Samp  ling 

Reduction  = 3 

Good 

target/background 

contrast; 

better 

detail 

Sampling 

Reduction  = 4 

Good 

ta  rget/background 

contrast; 

better 

detail 

Sampling 

Reduction  = 5 

Some  detail  lost  in  boom  section  of  aircraft.  Good 
target/background  contrast. 
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TABLE  B-18.  ADDITIVE  GAUSSIAN  SPOT  NOISE  APPLIED 
TO  A SYNTHETIC  APERATURE  RADAR  IMAGE 


Standard  Deviation  = .3%  of  Gray  Level  Range 

Standard  Deviation  = .6%  of  Gray  Level  Range 

Standard  Deviation  = 1.22%  of  Gray  Level  Range 

Standard  Deviation  = 2.44%  of  Gray  Level  Range 

Standard  Deviation  = 4.88%  of  Gray  Level  Range 
Image  appears  slightly  grainy. 

Standard  Deviation  = 9.77%  of  Gray  Level  Range 

Shape  of  aircraft  on  parking  ramp  distorted, 
graininess  increased  causing  some  loss  of  background 
detail  and  target  sharpness. 

Standard  Deviation  = 19.52%  of  Gray  Level  Range 

Detail  lost,  especially  in  hanger  area.  Aircraft 
shape  distorted.  Change  in  run-way  construction 
material  (on  over-run,  not  as  apparent  as  it  is  on 
original ) . 

Sta ndard  Deviation  = 39.04%  of  Gray  Level  Range 

Detail  almost  completely  lost  in  hanger  area  closest 
to  hard  stand.  Shape  of  aircraft  distorted  beyond  recog- 
nition. Complete  loss  of  background  detail,  road  not  as 
apparent  as  in  the  other  three  image  scenes  due  to  "salt 
and  pepper"  (spotty)  effects. 


Very  little 
(if  any) 
change  was 
noted  be- 
tween these 
four  image 
scenes  and 
the  original 
image. 
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TABLC  B-19.  ADDITIVE  GAUSSIAN  LINE  NOISE  APPLIED 
TO  A SYNTHETIC  APERATURE  RADAR  IMAGE 


Standard 

Deviation 

= .3%  of 

Gray  Level  Range 

Standard 

Deviation 

= .6%  of 

Gray  Level  Range 

Some 

di storti on 

in  shape 

of  aircraft 

noted . 

Standard 

Deviation 

= 1.22% 

of  Gray  Level 

Range 

Some 

di storti on 

in  shape 

of  aircraft 

noted . 

Standard 

Devi ati on 

= 2.44% 

of  Gray  Level 

Range 

Some 

d i s tort i on 

in  shape 

of  ai  rcraf  t 

noted . 

Standard 

Deviati on 

= 4.88% 

of  Gray  Level 

Range 

Lines  "breaking-up"  edges  of  hard  stand,  taxi-ways 
and  hangers.  Parked  A/C  distorted. 

Standard  Deviation  = 9.77%  of  Gray  Level  Range 

Lines  more  prominent.  Hangers  to  right  of  hard  stand 
may  be  interpreted  as  separate  building  rather  than  one 
long  building  because  of  the  line  noise.  Detail  in  parked 
aircraft  lost. 

Standard  Deviation  = 19.55%  of  Gray  Level  Range 

The  only  thing  that  an • i nterpreter  can  say  about  this 
image  scene  is  that  it  is  an  airfield  and  aircraft  parking 
area.  The  edges  of  the  hard  stand  are  still  recognizable. 
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TABLE  B-19.  ADDITIVE  GAUSSIAN  LINE  NOISE  APPLIED 
TO  A SYNTHETIC  APERATURE  RADAR  IMAGE  (Continued) 


Standard  Deviation  - 39.04%  of  Gray  Level  Range 

Line  noise  makes  any  detail  interpretation  impossible. 
Targets  are  "broken-up"  by  the  noise.  The  scene  can, 
however,  be  recognized  as  an  airfield  and  aircraft  parking 
area.  The  hangers  located  in  the  vicinity  of  the  parking 
area  (hard  stand)  are  no  longer  apparent. 


TABLE  B-20-  GAUSSIAN  LINE  JITTER  APPLIED  TO  A 
SYNTHETIC  APERATURE  RADAR  IMAGE 

Standard  Deviation  = .05  Pixels 
Standard  Deviation  = .1  Pixels 
Standard  Deviation  = .2  Pixels 

Standard  Deviation  - .39  Pixels 

Edges  of  run-way  and  taxi-ways  appear  "choppy" 
because  of  jitter. 

Imagery  is  still  useful  to  the  interpreter. 

Standard  Deviation  = .78  Pixels 

Target  shapes  are  distorted  due  to  jitter.  Straight 
lines  appear  extremely  "choppy." 

Standard  Deviation  = 1.56  Pixels 

The  target  scene  can  be  recognized  as  an  airfield; 
however,  all  detail  is  lost.  Parked  A/C  is  no  longer 
detectable  and  the  hangers  are  no  longer  recognizable. 

Standard  Deviation  = 3.13  Pixels 

The  runway  can  still  be  detected;  however,  I believe 
that  is  true  only  because  I am  familiar  with  the  image 
scene . 

Standard  Deviation  = 6.25  Pixels 
Same  as  above. 
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TABLE  B-21.  VARYING  THE  QUANTIZATION  OF  A 
SYNTHETIC  APERATURE  RADAR  IMAGE 


Number 

of 

Gray 

Levels  = 32 

Number 

of 

Gray 

Levels  = 16 

Number 

of 

Gray 

Levels  = 8 

"Gray-like,  blotchy  appearance"  takes  away  from 
target  detail.  This  is  especially  true  in  hanger  area 
to  the  right  of  the  A/C  parking  area.  In  addition,  much 
background  detail  is  lost  by  this  "bl o tch i nes s . " 

Number  of  Gray  Levels  = 4 
Same  as  above. 
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TABLE  B- 22 . VARYING  THE  SAMPLING  OF  A SYNTHETIC 
APERATURE  RADAR  IMAGE  UNDER  HIGH  ALIASING 

Sampling  Reduction  = 2 
Sampling  Reduction  = 3 

Parked  A/C  no  longer  detectable.  Some  loss  of  detail 
in  overall  scene. 

Sampling  Reduction  = 4 

Parked  A/C  not  detectable.  Great  loss  of  detail  in 
overall  scene. 

Sampling  Reduction  = 5 

Parked  A/C  and  hangers  not  detectable.  Target  scene 
can  be  recognized  as  an  airfield.  The  digitization  caused 
a "false-alarm”  situation.  The  aircraft  might  be  inter- 
preted as  a parked  A/C. 
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TABLE  B-23.  VARYING  THE  SAMPLING  OF  A SYNTHETIC 
APERATURE  RADAR  IMAGE  UNDER  REDUCED  ALIASING 


W 


Samp! i ng  Redu c t i on  = 2 

Parked  aircraft  not 
targets  and  background 
image  scene. 

Sampling  Reduction  = 3 

Parked  aircraft  not  detectable.  Poor  contrast  between 
targets  and  background  causing  loss  of  detail  in  overall 
image  scene. 

Sampling  Reduction  = 4 

Loss  of  contrast  and  overall  sharpness  causes  loss  of 
detail  in  any  targets  that  may  be  present  ’n  the  image 
scene . 

Samp  ling  Reduction  = 5 

The  target  scene  can  be  recognized  as  an  airfield; 
however,  all  detail  is  lost. 


detectable.  Poor  contrast  between 
causing  loss  of  detail  in  overall 
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APPENDIX  C 

INTERPRETATION  OF  MOIRE'  PATTERN  FORMATION 
AS  A DIGITAL  SAMPLING  PROCESS 


2*9 


APPENDIX  C 


INTERPRETATION  OF  MOIRE'  PATTERN  AS  A 
DIGITAL  SAMPLING  PROCESS 

Moire'  patterns  are  those  shimmering  patterns  seen  in 
everyday  life  when  looking  through  layers  of  screen  mesh,  folds 
of  silk,  or  even  parallel  bridge  railings.  Since  these  Moire' 
patterns  are  exceptionally  sensitive  to  displacements,  their  use 
has  been  abstracted  for  use  in  science  to  detect  optical  flat- 
ness, crystalline  inhomogeneities,  stress,  faint  signals,  etc. 

At  its  conceptual  simplist,  the  Moire'  pattern  is  formed 
by  viewing  one  periodic  structure  through  another  periodic 
structure.  Such  an  example  was  shown  previously  in  Figure  5-4. 

The  formation  of  the  Moire'  pattern  is  best  understood  in 
terms  of  the  spatial  frequency  representation  of  a digital 
sampling  process.  Figure  C-l  is  the  spatial  frequency  represen- 
tation of  the  test  image,  Figure  5-1.  The  test  image  is  actually 
composed  of  square  wavs  and  thus  possesses  an  infinite  number 
of  harmonics.  For  simplicity  and  with  no  serious  results  we 
ignore  the  ha rmon ics  in  this  discussion.  Now  remember  that 
the  test  image  was  seen  by  viewing  through  its  slightly  rotated 
self.  Figure  C-2,  giving  Figure  5-4.  Thus  the  transparent 
rotated  copy  of  the  test  image  forms  a spatial  sampler  of  the 
nature  discussed  in  Section  1 1 1 - A . The  resulting  sampling 
action  Is  one-dimensional.  This  implies  that  the  spatial 
frequency  spectrum  will  be  replicated  throughout  the  spatial 
frequency  plane  with  the  origin  moved  to  each  of  the  points 
as  shown  in  Figure  C-3.  The  results  of  the  replication  are 
seen  in  Figure  C-4,  the  spatial  frequency  representation  of 
Figure  5-4. 

From  Figure  C-5,  the  MTF  response  of  the  human  visual  sys- 
tem, it  is  apparent  that  certain  spatial  frequencies  are  more 
prominent  to  the  human  than  others.  By  moving  the  image  of 
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Figure  5-4  closer  to  or  further  from  the  eye  it  is  possible 
to  superimpose  contours  of  MTF  response  onto  Figure  C-5»  as 
shown  in  Figure  C-6.  Thus  the  spatial  frequency  components 
of  Fi gure  5-4  whi ch  are  most  apparent  are  D,  E,  and  F.  E is 
the  average  brightness  while  D and  F form  an  aliased  low 
frequency  nearly  vertical  spatial  sine  wave.  This  sine  wave 
is  the  quite  evident  Moire'  pattern  observed  in  Figure  5-4. 
Thus  Moire'  patterns  are  seen  to  be  the  aliased  phenomena  of 
a digital  image  sampling  process  and  are  subject  to  the  laws 
and  properties  thereof.  For  instance  the  reason  that  the 
Moire'  pattern  is  horizontal  while  the  original  patterns  are 
vertical  is  apparent  in  the  spectral  plane  of  Figure  C-6. 

Very  complex  and  rich  forms  of  digital  image  testing 
procedures  can  thus  be  developed  using  the  specialized  alias- 
ing properties  of  digital  image  systems  and  the  MTF  properties 
of  the  eye.  Section  V considers  a few  of  the  possibilities. 
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APPENDIX  D 

MATHEMATICS  OF  MOIRE1  PLOTTING 

A mathematical  understanding  of  the  means  by  which  Moire1 
patterns  can  be  used  to  plot  function  is  developed  here. 

Let  the  lines  of  the  test  image,  Figure  5-9,  consisting 
of  horizontal  lines  be  indexed  by  the  sequenced 

h = -2,  -1,  0,  1,  2,  

and  the  lines  of  Figure  be  indexed  by  the  sequence 

k = ....  -2,  -1,  0,  1,  2 

as  shown  there.  It  will  be  seen  in  Figure  that  the  Moire' 

patterns  will  be  observed  at  the  locus  of  the  points  of  inter- 
section of  the  two  patterns  of  Figures  5-9  and  5-10  where  h-k 
is  another  integer  p 

p = -2,  1,  0,  1,  2 

That  is,  the  Moire'  patterns  are  solutions  of  the  subtractive 
indicial  equation 

h - k = p ( D - 1 ) 

Now  the  equation  for  the  lines  of  Figure  is 

y = bh  + tf  (bh)  (D-2) 

Where 
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y 

b 

h 

t 

f(.) 


the  vertical  position  of  any  point  on 
the  line 

the  normal  spacing  of  the  lines 

is  defined  above 

a constant,  the  scale  factor 

the  function  by  which  the  line  spacing 
deviates  from  the  normal  spacing,  b 


The  equation  for  the  lines  of  the  figure  consisting  of 
horizontal  lines  when  canted  at  angle  a is 


x = y cot  a - bk  cos  a/  sin  a 


(D-3) 


where 


x = the  horizontal  position  of  any  points  on 
the  line  and  all  other  values  are  defined 
above 

Equation  (d - 2 ) can  not  in  general  be  solved  for  h.  Instead, 
solving  Equation  (D-3)  for  k gives 

k = ^ C_Ota  - X ( 4) 

b tan  0 ' 

Inserti ng  Equation  (D-4)  into  Equation  (D-l)  gives 

hb  . _JL c.P_sg_- ^ — n_a_  + pb  (d_5) 

Now  rotating  the  coordinate  system  clockwise  throuqh  an  angle  a 
gives  the  old  coordinates  x,  y in  terms  of  the  new  coordinates 
x1  , y* 


»> 

a 

>> 

cosa  + x ' sina 

( D - 6 ) 

x = y 1 

s i na  + X ‘ COSa 

(D  - 7 ) 
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Combi ni ng  Equations  (D-6)  and  (D-7)  into  Equation  (D-5)  gives 


hb  = ■ ^ — + pb 

COSa  r 


( D-8 ) 


Substituting  Equations  (D-8)  and  (D-6)  into  Equation  (D-2)  gives 
_ 1 


sina 


^ — - COSa)  + pb  + tf|— ^ — — + b p| 
y COS  a J K ycosa  H f 


Now  as  a -*•  0 


x ' * a [pb  + tf (y +bp )J 

Upon  appropriately  shifting  the  origin  and  dropping  primes,  the 
Moire'  pattern  obeys  the  equation 

x‘  £ f(y) 

Thus  it  is  seen  that  the  Moire'  pattern  plots  the  function  f(-) 
at  an  amplification  inversely  proportional  to  a,  the  angle  of 
cant,  of  the  original  patterns  of  Figures  5-9  and  5-10. 
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APPENDIX  E 


MEASUREMENT  OF  MTF  BY  THE  OBSERVER 

Section  V-B  described  an  approach  by  which  the  MTF  of  a 
digital  system  could  be  determined  by  an  observer  without  test 
equipment.  The  mathematical  justification  is  given  here. 

Let  AjN(x,y)  and  AgUT(x,y)  be  the  amplitude  of  the  sine 
wave  at  location  (x,y)  in  the  original  and  reproduced  test 
images . 

Now  from  Equation  ( 1 ) x is  also  the  spatial  frequency  of 
the  sine  wave  at  the  point  (x,y).  Thus 

Ajn  (x,y)  = y 

Now  from  the  definition  of  MTF  it  is  seen  that 

^OUT  ^ x ) = MTF(x)  A j n ( x ,y ) 

= MTF ( x )y 

For  any  given  conditions,  there  exists  a minimum  amplitude 
d(x)  of  any  spatial  sine  wave  of  frequency  which  any  given  human 
can  perceive.  Thus  there  exists  a b(x).  Figure  5-7,  such  that 

A0UT*x,Mx^  = d(x) 
and  there  exists  an  a(x)  such  that 

AIN(x,a(x))  = d ( x ) 


It  follows 


A0UT (x  »b (x  ) ) = d(*)  = AJN(x,a(x)) 
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or 

MTF ( x ) b ( x ) = d ( x ) = AjN(x,a(x)) 
or 

MTF(x)b(x)  = a(x) 

or 

MTF<*>  * 

which  is  the  results  obtained  at  Step  4 of  the  procedure  of 
Section  V-B. 
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